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s  The  discovery  that  fracture  surfaces  exhibited  limited  areas  with 

unusual  dimple  features,  interpreted  to  have  been  caused  by  high 

temperatures,  led  to  research  on  the  temperature  rise  of  moving  crack 

tips.  Experimentally,  this  entailed  mechanical  measurements  using 

-3  -1  3  -1 

strain  rates  between  10  sec  and  10  sec  ,-*the  determination  of 

y  x...  ” 

stress  intensity  rates,  straining  in  vacuum,  measurement  of  fracture 
topographies  by  stereophotogrammetry ,  and  the  recording  of  light 

emission  during  fracture.  Three  theoretical  models  were  developed  to 

/ 

explain  the  above  hypothesis;  the  first  explored  the  extent  of 
thermal-mechanical  coupling  by  finite  element  analysis,  the  second 
developed  the  geometry  of  adiabatic  shear  for  plane  stress  specimens, 
and  the  third  consisted  of  a  numerical  analysis  of  crack  tip 
temperatures  based  on  a  new  development  of  microstructural  plasticity 
for  fast  moving  cracks.  Results  from  all  phases  of  this  research 
provided  support  for  the  contention  that  the  molting  point  can  be 
exceeded  at  certain  local  areas  along  the  crack  flank  for  the  high 
strength,  low  thermal  conductivity  titinium  alloys  investigated. 
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INTRODUCTION  AND  OVERVIEW 

The  objective  of  this  project  was  to  develop  an  understanding  of 
the  influence  of  microstructure  on  the  fracture  of  ductile  titanium 
alloys  and  specifically  on  their  fracture  toughness.  The 
microstructures  of  most  titanium  alloys  are  complex  and  their 
relationship  to  fracture  toughness  is  rather  involved  as  reported 
earlier  by  Hirth  and  Froes  [1], 

The  project  was  started  with  a  study  of  void  initiation  in  P.  C. 
titanium  by  correlating  tensile  data  with  detailed  observations  of 
fracture  surfaces  examined  by  scanning  electron  microscopy  (.SEM)  .  Grain 
size  was  chosen  as  structural  variable,  and  a  range  from  13  ym  to  1310 
ym  was  covered.  In  the  past,  the  3-D  character  of  fracture  surfaces  had 
been  neglected  and  most  measurements  were  taken  from  the  2-D  projection. 
Our  work  was  done  with  SEM  stereo-pairs  and  the  features  of  fracture 
surfaces  were  measured  with  a  light  spot  stereoscope.  Concentrating  on 
void  initiation  during  deformation  leading  to  ductile  fracture  by  void 
coalescence,  it  was  found  that  up  to  grain  sizes  of  about  50  ym  the 
primary  void  initiation  site  was  at  grain  boundary  triple  points.  At 
larger  grain  sizes,  void  initiation  sites  were  seen  at  slip  band 
intersections,  interaction  of  dislocation  avalanches  with  grain 
boundaries,  grain  boundary  ledges,  twin  intersections,  and  triple 
points.  By  reconstructing  the  topography  of  matching  fracture  surfaces 
it  was  possible  to  determine  the  sequence  of  void  initiation  and  the 
subsequent  growth  of  voids  towards  their  coalescence  for  specific  cross 
sect  ions . 

Turning  to  titanium  alloys  the  research  focused  on  crack 
propagation  selecting  Ti- 10V-2Fe-3Al  for  plane  strain  fracture  and 
Ti-SMn  for  the  plane  stress  mode.  Our  early  SEM  stereo  observations  of 
fracture  surfaces  revealed  in  both  alloys  limited  areas  of  droplet-like 
features  at  the  dimple  walls.  The  two  most  obvious  interpretations  led 
us  to  believe  that  these  features  wore  due  to  an  artifact  from  oxidation 
or  that  the  melting  point  hail  been  exceeded  at  the  moment  of  final 
separation  and  the  droplet- 1 i ke  dimple  wall  remnants  represented 
solidified  molten  material.  The  experimental  approach  involved 
iractures  in  vacuum  and  argon,  exposure  to  short  heat  I  lashes,  and 
attempts  to  record  radiation  possibly  coming  from  the  fracture  location. 
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Theoretically,  calculations  were  made  initially  to  obtain  a  ball  park 
figure  to  see  whether  indeed  a  temperature  rise  exceeding  the  melting 
point  at  a  few  local  areas  during  the  final  separation  was  a  reasonable 
conclusion  and,  if  so,  should  be  pursued  further.  A  summary  of  these 
early  deliberations  is  provided  in  Appendix  I  [2].  In  that  paper  the 
transformation  of  plastic  deformation  during  straining  preceding 
fracture  into  heat  was  estimated  as  a  first  contribution,  but  was  found 
negligible.  Following,  an  equation  developed  by  Rice  and  Levy  [3]  was 
adapted  to  our  situation  with  the  result  that,  theoretically,  a 
temperature  of  more  than  1,000°C  at  the  tip  could  be  generated  during 
crack  propagation.  And  finally,  dislocation  theory  was  applied  to  the 
rupture  of  micro  ligaments  which  always  develop  between  voids  in  ductile 
metals,  this  being  the  last  step  in  the  fracture  process;  an  additional 
300°C  -  50u°t  can  be  generated  locally  by  this  mechanism.  Adding  the 
throe  temperature  contributions  reassured  us  that  our  interpretation  of 
the  experimental  results  as  a  melting/solidification  phenomenon  was 
reasonable.  Details  of  microstructural  effects  at  the  crack  tip  in 
general  are  provided  in  Appendix  II  [4],  In  addition  to  two  Masters 
Theses  [5,6],  two  Ph.D.  Theses  have  documented  by  various  research 
techniques  \  7 , 8)  our  hypotlies  is  . 

One  can  anticipate  that  the  development  of  heat  at  the  crack  tip 
depends  on  the  rate  of  crack  propagation.  This  was  explored  with  the 

metastable  beta  Ti - 10V-2Fe-3A 1  alloy  which  was  subjected  to  strain  rates 
from  quasistatic  (10  ^s  to  near,  ballistic  (10^s  in  tensile 

loading.  Yield  stress,  ultimate  tensile  stress,  reduction  in  area,  and 
strain  to  fracture  increase  with  increasing  strain  rate,  z.  With  a  four 
order  of  magnitude  increase  of  z  the  fracture  toughness  was  seen  to 
increase  by  82°0.  The  tortuosity  of  the  tensile  fracture  surfaces  was  of 
interest  in  this  connection  and  changes  due  to  the  increasing  strain 
rates  were  measured  with  a  semi-automatic  stereographic  mapping  system 
(built  by  Bryant  [9]).  In  order  to  assess  the  magnitude  of  the  thermal 
mechanical  coupling  at  the  crack  tip,  a  finite  element  analysis  was 

carried  out  which  is  summarized  in  this  report  (pp.35-bh).  Ip  should  be 
add-’d  that  the  boundary  conditions  for  the.  finite  element  method  were 
obtained  from  experimental  measurement  on  Ti-10-2-3.  This  work 
confirmed  the  high  temperature  rise  at  the  crack  tip. 

The.  alpha/beta  alloy  Ti-HMn  was  investigated  in  sheet  form  with 

t  nipii.  I,  i  s  on  the  ve  i  i  l  i  t  icn  of  the  "m  i  croroughoued  or  droplet-like 
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features  on  the.  fracture  surfaces  by  heat  and  the  possible,  involvement 
of  adiabatic  shear  including  the  corresponding  glide  geometries.  The 
first  task  involved  the  exposure  of  regular  dimples  to  flashes  (0.5  sec) 
of  a  neutral  hydrogen/oxygen  torch  flame.  Re-examination  of  these 
treated  surfaces  showed  indeed  a  resemblance  of  the  "mi croroughened" 
features  observed  after  fracture.  Since  surface  oxides  and  nitrides 
rapidly  form  on  titanium  above  650°C,  it  was  suggested  that  oxidation 
might  be  responsible  for  the  observed  "microroughening. "  To  investigate 

this  possibility,  tensile  tests  were  performed  in  a  chamber  which  had 

_  ^ 

been  evacuated  to  10  Pa  and  backfilled  with  argon  four  times. 
Although  oxygen  and  nitrogen  were  not  completely  absent  in  these  tests, 
the  short  duration  of  the  high  temperatures  and  the  rarified  nature  of 
the  concentration  of  these  two  gases  severely  limits  the  extent  of 
oxidation  reactions.  Examination  of  the  resulting  fracture  surfaces 
revealed  the  same  features  as  the  air-fractures  samples.  This  result 
was  considered  a  clear  indication  that,  although  some  oxidation  of  newly 
formed,  heated  titanium  surfaces  is  inevitable,  the  surface  structures 
found  on  our  samples  do  not  form  primarily  by  an  oxidation  mechanism. 
The  geometrical  aspects  of  adiabatic  shear  are.  summarized  in  this  report 
on  pp .  67  to  98,  and  a  brief  account  of  visible  light  coming  from  the 
fracture  area  is  included.  All  results  are  providing  evidence  for  a 
high  temperature  rise  ahead  of  the  crack  tip. 

Our  earlier  analysis  of  t.lie  temperature,  rise  at  the  crack  tip  was 
based  on  the  Rice-Levy  equation  which  does  not  include  the  effect  of  a 
heat  increase  simultaneous  with  the  moving  crack.  Contrary  to  the 
earlier  models  by  otho.r  authors  who  based  their  calculations  on 
continuum  concepts,  our  model  (executed  by  Dr.  K.  Jagannadham,  North 
Carolina  State  University,  Raleigh,  N.C.,  under  subcontract)  is  making 
use  of  the  plastic  work  term  which  in  our  opinion  must  change  in  the 
region  ahead  of  the  crack  tip  during  crack  propagat ion .  Appropriate 
parameters  used  iri  the  model  are  dislocation  density  in  the  process 
none,  f  r  iet  iona  1  ‘•t  ress  for  dislocations,  dislocation  velocity,  plastic 
strain,  plastic  work  rate,  and  loading  rate.  Thus  the  workhardening 
cl.  l  r  i..  t  ■■  r  i  >  t  i  c  in  tie-  "p- las  Lie  /.one"  in  included  in  the  model.  The 
resells  of  extensive  computer  calculations  ire  relating  a  verv 
substantial  temperature  rise  in  tile  process  none  to  the  time  scale  of 
c  rack  propag.il  i  on  . 

In  summary,  the  proicct  research  explored  the  effect  of 
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microstructure  on  the  ductile  fracture  of  C.P.Ti,  Ti-SMn,  and 
T i - 1 0 V - 2 Fe - 3 A 1 .  It  was  found  that  the  temperature  rise  ahead  of  the 
crack  tip  can  be  more  than  1,000°C,  more  likely  1,300°C;  when  adding  to 
this  value  the  heat  developed  during  the  rupture  of  micro  1 igaments 
between  voids,  the  melting  point  of  these  alloys  was  exceeded  in 
localized  areas  of  the  crack  flank.  This  temperature  regime  for  the 
temperature  rise  at  the  crack  tip  could  be  calculated  since,  for  the 
first  time  in  a  high  temperature  metal,  a  determination  of  the  final 
temperature  at  fracture  was  possible  with  better  than  average  precision. 
It  should  be  remembered  that  titanium  alloys  have  very  low  heat 
conductivity  and  high  strength  which  makes  it  possible  for  adiabatic 
she.ai  to  occur  at  normal  loading  rates.  With  increasing  loading  rates 
the  fracture  toughness  was  found  to  increase  by  more  than  SCTo. 
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SECTION  II 

VOID  INITIATION  PRECEDING  DUCTILE  FRACTURE 
IN  POLYCRYSTALLINE  TITANIUM 

Based  on  a  Masters  Thesis  (Materials  Science-)  by  Marjorie  A.  Erikson*, 
UVa.  1983. 


ABSTRACT 


The  objective  of  this  research  was  a  determination  of  primary  void 
initiation  sites  preceding  void  growth  in  the  process  of  deformation 
leading  to  ductile  fracture  in  polycrystal  line  C.  P.  titanium. 
Stereophotogrammetry  was  employed  to  evaluate  SEM  fracture  surfaces. 
Two  distinct  regimes  of  void  initiation  were  found  to  operate. 

(i)  Void  initiation  occurred  primarily  at  grain  boundary  triple 
points  in  specimens  with  grain  sizes  up  to  50  pm; 

(.iil  In  material  with  larger  grain  sizes  additional  void 

initiation  sites  were  found  to  be  at  slip  band 
intersections,  interaction  of  dislocation  groups  with  grain 
boundaries,  twin  intersections,  and  grain  boundary  ledges. 

Quantitative  topographical  information  was  obtained  from 
measurements  with  a  light  spot  stereoscope.  Reconstructing  the 
topography  of  matching  fracture  surfaces  with  SEM  accuracy  and 
resolution,  the  time  sequence  of  void  initiation  along  a  line  on  the 
fracture  surface  was  obtained. 
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VOID  INITIATION  PRECEDING  DUCTILE  FRACTURE 
IN  POLYCRYSTALLINE  TITANIUM 

1.  INTRODUCTION 

While  void  initiation  has  been  covered  in  fair  detail  for  particle 
containing  alloys  [1,2],  experimental  evidence  in  pure  or  relatively 
pure  polycrystalline  ductile  metals  has  received  only  scant  attention, 
except  for  fracture  at  higher  temperatures  [3].  Under  creep  conditions 
grain  boundary  sliding  often  leads  to  void  formation  at  grain  boundary 
triple  points  (.GBTP)  and  theoretical  models  for  the  high  temperature 
regime  have  been  published  [4,5].  One  cannot  expect  that  these 
mechanisms  are  operative  in  high  melting  point  metals  at  room 
temperature,  and  plastic  deformation  has  been  invoked  for  void 
initiation  by  assuming  the  development  of  pile-ups  when  dislocations 
were  unable  to  penetrate  grain  boundaries  [6,7,8],  Also,  grain  boundary 
ledges  have  been  thought  to  initiate  microcracks  [9];  other  mechanisms 
are  based  on  crossing  slip  bands  [10],  or  intersecting  twins  [11],  and 
bend  planes  [7,12] . 

The  present  investigation  has  as  its  objective  to  determine  the 
locations  of  primary  void  initiation  sites.  Fracture  experiments  were 
made  with  C.P.  titanium  under  tensile  loading,  and  fracture  surfaces 
were  examined  with  a  scanning  electron  microscope  (SEM).  The  crack  path 
in  most  metals  is  a  tortuous  one,  and  this  is  particularly  true  for 
titanium.  It  was  found  necessary  to  employ  stereo-photogrammetry  in 
order  to  obtain  true  measurements  of  fracture  surface  features. 

2.  SAMPLE  PREPARATION  AND  EXPERIMENTAL  PROCEDURE 

Tensile  samples  were  prepared  from  a  0.198  cm  thick  sheet  of  C.P. 
titanium.  The  material  was  99.838°t>  pure  with  an  impurity  content  that 
included  0,  N0 ,  C,  H9  and  Fe,(.  After  polishing  the  gauge  sections  of 
each  sample  with  grit  size  decreasing  to  600  grit  and  heat  treating,  the 
gauge  lengths  and  cross-section  areas  of  each  sample  were  carefully 
measured.  Once  fracture  occurred  pieces  containing  the  fracture 
surfaces  were  carefully  cut  approximately  0.5  cm  from  the  fracture  edge 
.lnd  mounted  onto  an  aluminum  specimen  stub  for  subsequent  SEM  study. 

A  range  of  heat  treatments  was  done  to  vary  the  grain  size  of 
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samples  to  allow  a  study  of  grain  size  versus  tensile  properties  and 
fracture  surface  characteristics.  Because  of  sample  purity  only  the 
a-phase  of  Ti  is  present  at  room  temperature,  so  stress  relief, 
recrystallization  and  grain  control  were  the  primary  goals  of  heat 
treating. 

Ten  samples  were  prepared  for  study.  Heat  treatments  were  designed 
to  cause  a  variation  in  grain  size  from  13  pm  to  1300  pm.  Table  1  lists 
the  various  samples  with  the  heat  treatment  they  were  given  and  their 
resulting  grain  size. 

To  prevent  oxygen  from  diffusing  into  the  samples  during  heat 
treating,  each  sample  was  enclosed  in  a  quartz  tube  evacuated  to  5x10  ^ 
torr . 

For  grain  size  determination,  a  piece  from  the  shoulder  of  each 
sample  was  cut,  mounted  onto  an  SEM  stub,  encased  in  epoxy,  and  polished 
and  etched  in  a  solution  of  10  ml  HF,  5  ml  HNO^  and  85  ml  distilled 
water  for  1|  minutes. 

The  tensile  specimens  were  pulled  to  fracture  on  an  Instron 
Universal  Testing  Machine.  The  crosshead  speed  was  either  0.05  cm/sec 
or  0.005  cm/sec.  Load  vs.  time  was  recorded  and  from  this  tensile 
properties  were  determined. 

Fractographs  of  the  entire  surface  were  taken  of  opposing  surfaces 
at  low  magnifications,  using  stereo  techniques,  as  references  for  future 
higher  magnification  micrographs.  These  low  magnification  fractographs 
were  also  used  to  measure  final  cross-sectional  area  at  the  fracture  to 
determine  the  reduction  in  area  of  each  sample. 

A  number  of  pictures  were  taken  at  higher  magnifications  (3G0x  to 
1000.x),  using  stereo  techniques  to  enable  a  detailed  three-dimensional 
study  of  the  fracture  surface.  Thus,  fracture  surface  features  could  be 
determined  accurately  by  making  measurements  of  relative  height  as 
well  as  relative  horizontal  distances. 

Studies  of  microstructures  were  related  to  studios  of  fracture 
surface  features  to  pinpoint  primary  void  initiation  sites  in  a-Ti. 
Five  to  ten  micrographs  of  each  sample,  polished  and  etched  for  grain 
size  measurements,  were  enlarged  and  used  to  determine  grain  size  using 


TABLE  1 


GRAIN  SIZES  OF  FRACTURE  SPECIMENS 


SAMPLE  HEAT  TREATMENT  Average  Grain  Size 


1 

i  hr.  at  700°C 
furnace  cool 

13  ym 

1-55 

ym 

2 

i  hr.  at  700°C 
air  cool 

21  ym 

2-72 

ym 

3 

2  hrs .  at  700°C 
air  cool 

46  ym 

7-188 

ym 

4A 

2  hrs.  at  1000°C 
air  cool 

178  ym 

30-1210 

ym 

4B 

2  hrs.  at  1000°C 
air  cool 

200  ym 

27-1215 

ym 

5A 

5  hrs.  at  1175°C 
furnace  cool 

427  ym 

40-1550 

ym 

5B 

5  hrs .  at  1 150°C 
furnace  cool 

420  ym 

35-1490 

ym 

6 

8 i  hr.  at  1150°C 
furnace  cool 

792  ym 

70-1880 

ym 

7A 

20  hrs .  at  1150°C 
furnace  cool 

1247ym 

82-2150 

ym 

7B 

20  hrs.  at  1150°C 
furnace  cool 

1 3 lOym 

79-2120 

ym 

the  grain  boundary- 1 ine  intercept  method.  These  pictures  were  also  used 

to  measure  intertriple  point  distances  by  direct  measurements. 

2 

Measurements  were  made  over  areas  of  ~1000  ym  for  small  grained  samples 
2  2  2 

and  1000“  -  2000“  ym“  for  larger  grained  samples. 

Measurements  of  dimple  diameter  and  interdimple  distances  were 
taken  directly  from  stereo  micrographs  which  allowed  evaluations  in 
3-dimensions.  Dimple  diameters  were  taken  as  the  line  through  the 
center  of  the  dimple  perpendicular  to  the  long  axis  of  the  dimple. 
Interdimple  distances  were  measured  from  the  center  of  one  dimple  to  the 
center  of  an  adjacent  dimple  taking  relative  height  variations  into 

account  bv  making  measurements  using  a  Hilgor-Vntts  stereo  viewer. 

3.  RESULTS 

Significant  mechanical  properties  of  typical  tensile  samples  are 
listed  in  Table  2 . 

An  evaluation  of  measurements  obtained  from  SEN  fractographs 
indicated  immediately  that  in  small  grain  specimens  interdimple 

distances  corresponded  primarily  to  intertriple  point  distances.  This 
can  be  clearly  seen  in  Figure  1;  the  measurements  were  taken  from 
specimens  1-3  with  average  grain  sizes  of  13  ym,  21  ym,  and  46  ym, 
respectively.  The  fit  of  the  second  peak  of  interdimple  spacings  with 
intertriple  point  spacing  peaks  is  amazingly  good.  Notably,  there  is  on 
each  of  these  graphs  a  higher  peak  for  interdimple  spacings  which  lies 
between  1  ym  and  3  ym .  The  significance  of  this  maximum  peak  will  be 
discussed  in  the  following  section.  Figure  2  gives  an  impression  of 
dimple  size  distributions  in  a  fracture  surface  of  specimens  having 

grain  sizes  below  50  ym . 

The  character  of  fractographs  changed  drastically  for  specimens 
with  grain  sizes  of  200  ym  and  larger  (Figure  3).  Instead  of  only  a 
broad  distribution  of  dimple-  sizes  over  the  whole  fracture  surface, 

elongated  ridges  are  also  seen  over  a  relatively  large  area.  Dimple 
size  distributions  tor  tour  grain  sizes  are  given  in  Figure  4.  For 
these  samples,  iuLerdimple  distances  and  intertriple  point  distances  are 
unrelated  although  the  maximum  peak  of  the  first  group  is  retained  at  1 
ym  to  3  yin.  The  ch  mge  in  the  fracture  surface  characteristics  between 
specimens  with  "small"  and  "large"  grains  indicates  the  presence  of  two 
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lgnre  _.  Sh.'i  fractograph  of  specimen  «2  having  an  average  grain  size  of 
21  pm.  Range  of  grain  sizes  from  2  to  72  pm.  The  bottom 
levels  of  dimples  A,  C  and  B  are  at  10  pm,  5  pm  and  0  base, 
respectively,  as  determined  by  stereo-photogrammet ry .  The 
ridge  on  the  right  of  dimple  11  (arrow)  is  almost  40  pm  higher 
than  the  bottom  level  of  II. 
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distinct  mechanisms  of  void  initiation. 

The  measurements  displayed  in  Figure  1  and  4  were  made  from  stereo¬ 
pairs.  A  simple  addition  to  a  flying  spot  stereo-microscope  made  it 
possible  to  draw  topographical  maps  of  fracture  surfaces.  An  example  is 
given  for  a  rather  simple  fractograph,  Figure  5,  for  which  a  stereo-pair 
is  shown.  The  corresponding  topographical  map  is  given  in  Figure  6. 
While  stereo-viewing  conveys  a  distinct  three-dimensional  impression  of 
the  general  layout  of  a  fracture  surface  and  also  of  individual  dimple 
geometries,  it  would  not  have  been  possible  to  conceive  the  large 
difference  between  the  lowest  points  of  the  dimples.  The  largest 

height  difference  measured  in  this  micrograph  amounts  to  40.8  pm. 
Dimple  has  an  average  depth  of  24  pm  while  dimples  42  and  43  have 
average  depths  of  13  pm  and  10  pm,  respectively.  The  stereo-photo- 
grammetry  revealed  not  only  that  the  dimple  depths  are  different  by 
larger  amounts  than  what  the.  eye  had  perceived  but  that  the  base  surface 
is  inclined,  i.e.,  it  increases  from  top  right  of  Figure.  5  to  bottom 
left.  The  even  greater  complexity  of  the  topography  in  large-grained 
specimens  is  now  appreciated  and  stereo  pairs  are  essential  for  their 
eva  1  tint  i  on  . 

For  <i  complete  evaluation  of  fracture  surfaces,  however,  one  has  to 
go  one  step  further.  From  the  irregular  features  seen  in  all 

stereo-micrographs,  one  must  draw  the  conclusion  that  the  growth  of 
voids  depends  on  local  stress  conditions  at  a  mi crostructura  1  scale  and 
consequently  is  anisotropic.  Therefore,  one  can  expect  that  matching 
fracture'  surfaces  will  look  differently,  and  the  drawing  of  contour  maps 
has  shown  this  convincingly.  Figure  7a  depicts  a  cut  through  contour 
maps  of  matching  fracture  surfaces.  In  order  to  classify  the  difference 
in  the  features  of  these  matching  surfaces,  the  vertical  scale  has  been 
enlarged  by  a  factor  of  five.  The  true  profile  for  the  "lower"  surface 
i  shown  in  Figure  7b. 

PISbWSSivN 

The  i  mmed  i  it  e  problem  is  to  rationalize,  the  effects  of  grain 
i  •  1  iry  -.tru  lures  on  void,  initiation.  The  dimple  measurements  show 
i kit  there  ire  two  relationships  between  i nterd imp le  distances  and 
i nt  ert  r : p 1 e  punt  distances:  one  group  with  grain  sizes  below  30  pm  has 
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an  almost  perfect  correspondence  between  these  data  points,  while  for 
another  group  with  grain  sizes  above  200  ym  this  relationship  is  not  the 
dominating  one. 

Crack  initiation  at  triple  points  is  common  in  creep  as  first 
explored  by  Zener  [13];  here  grain  boundary  glide  was  involved  and 

thereby  a  triple  point  becomes  the  obvious  location  for  a  stress 
concentration.  Movement  between  grains  along  their  boundary  is  not  to 
be  expected  in  titanium  at  room  temperature.  However,  incompatibility 
stresses  between  grains  [14}  should  be  high  at  triple  points.  In-s itu 
straining  in  a  TEM  showed  that  the  first  glide  dislocations  are  often 

generated  at  triple  points  (see  Figure  8  [15]),  aud  slip  lines  have  been 

seen  to  originate  at  triple  points  [16].  Recently  Kurzydlowski  et  al. 
concluded  on  the  basis  of  in-situ  TEM  work  in  stainless  steel  that  the 
emission  of  dislocations  from  a  triple  point  occurs  at  a  stress  of 

G/320,  with  G  the  shear  modulus  [17],  These  authors  claim  that  the 
stress  concentration  factor  of  a  triple  point  is  not  far  from  2.8 
[17,18].  All  of  this  indicates  that  after  extensive  deformation,  a 
triple  point  represents  an  exceptionally  high  local  stress  concentration 
in  a  polycrvstalline  material.  In  addition,  it  is  known  that  impurities 
reduce  the  energy  in  grain  boundaries  [19]  and  this  most  likely  applied 
to  our  case.  The  experimental  result  that  grain  boundary  triple  points 
have  been  the  primary  sites  for  void  initiation  is  thus  understandable. 
The  question  now  to  be  treated  is  concerned  with  the  observation  that 
for  specimens  with  grains  larger  than  200  pm  this  does  not  hold. 

Since  the  reduction  in  area  for  all  specimens  was  60°o  or  higher, 
the  plastic  deformation  in  the  neck  was  substantial.  Titanium  deforms 
at  room  temperature  by  prism  glide,  {1010}  <1120>,  and  by  twinning  on 
six  planes  [20];  the  most  frequently  observed  twins  have  {1012}  and 
{1122}  as  composition  planes.  The  influence  of  twinning  on  void 
initiation  has  to  be  included  in  our  deliberations  since  it  is  well 

known  that  twinning  plays  an  important  role  in  larger  grained  titanium 

while  deformation  twins  have  not  been  seen  in  small  grained  titanium 

[21],  Mechanical  twinning  will  occur  in  Ti  for  {1012}  ^ 1 0 1  1 > ,  for 

example,  by  stresses  in  tension  parallel  to  the  c-axis  and 

correspondingly  for  compression  stresses  perpendicular  to  it.  That  both 
types  of  stresses  can  act  at  grain  boundaries  cannot  be  doubted  since' 
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after  even  modest  strains,  grain  boundaries  will  have  developed  ledges, 
the  stress  fields  of  which  have  been  calculated  by  Das  and  Marcinkowski 
[9],  Amateau  et  al.  reported  the  occurrence  of  twins  in  a-Ti  (0^  = 
0.25°o  as  compared  to  our  0n  =  0.102°o)  after  a  strain  of  4.8°0  for 
specimens  with  grain  sizes  comparable  to  our  specimen  #4  [22].  These 
authors  identified  "thick"  twins  as  {1012}  and  {1122}  types  and  "thin" 
twins  of  {112x}  types  with  x  =  1,  3,  or  4  (second  order  twins)  in 
accordance  with  the  findings  of  Rosi  et  al.  [23].  The  second  order 
tw'ins  were  seen  in  the  {1122}  twins  and  the  crack  nucleation  occurred 
between  second  order  twins  and  matrix  based  on  observations  by  light 
microscopy  at  180°K.  Figure  9  shows  a  fracture  surface  in  support  of 
this  mechanism.  In  stereo-viewing  one  can  see  a  fine  ridge  structure 
with  a  spacing  of  a  few  micrometers  lying  within  larger  bands  10  pm  to 
30  pm  wide.  The  latter  fit  the  description  of  the  "thick"  {1122}  twins 
and  the  finer  lamellae  the  secondary  twins.  A  grain  boundary  is 
indicated  by  arrows.  Intergranular  fracture  obviously  has  produced  the 
relatively  flat  area  marked  A.  Arrows  B  are  pointing  to  two  deep 
microcracks.  While  Figure  9  represents  an  extreme  case  of  crack 
initiation  by  twinning  processes,  Figure  10  is  indicative  of  slip  and 
twinning  which  are  both  present  in  areas  about  a  few  hundred  micrometers 
across  where  void  initiation  by  glide  and  initiation  of  microcracks  by 
twinning  operated  in  close  proximity  in  grains  between  200  and  1,300  pm. 

In  Figure  11  we  have  plotted  the  maxima  of  interdimple  spacings 
against  grain  size  and  one  can  see  that  the  interdimpde  spacing  follows 
the  grain  size  for  the  grains  and  then  levels  off  regardless  of  grain 
size.  All  measurements  of  interdimple  spacings  have  been  made  from 
only  those  areas  of  fracture  surfaces  which  showed  full  ductility.  As 
seen  from  Figure  4,  the  broad  distribution  of  dimples  is  not  uniform  but 
shows  in  addition  to  the  maximum  of  about  50  pm  a  number  of  smaller 
peaks.  The  same  is  true  for  the  intertriple  point  spacing  with  its 
maximum  near  the  average  grain  size.  It  should  be  noted  that  for  the 
four  grain  sizes  (Figure  4)  a  small  peak  of  intertriple  point  spacings 
coincides  with  a  peak  of  interdimple'  spacings.  Again  the  indication  is 
that,  triple  points  remain  important  void  initiation  sites.  However, 
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other  void  initiation  processes  are  now  more  competitive.  Besides 
initiation  sites  at  grain  boundary  ledges  and  secondary  twin-matrix 
intersections,  other  processes  must  have  been  active.  The  "bend  plane" 
mechanism  is  normally  cited  as  an  important  crack  nucleation  model  for 
hexagonal  crystals  but  is  not  applicable  for  highly  strained 
polycrystalline  Ti  since  that  mechanism  is  based  on  bicrystal  studies  in 
zinc  at  low  temperatures  where  cracks  appeared  parallel  to  the  basal 
plane  after  elongations  to  fracture  of  only  0 . 5°0  to  2 . 8°<,  [12].  All  of 
our  samples  showed  multiple  glide  in  the  necked  region  which  means  that 
extensive  glide  was  produced  by  a  number  of  glide  systems  crossing  each 
other.  Microcrack  initiation  at  slip  intersections  in  close  packed 
crystals  has  been  considered  as  one  of  the  most  important  void 
initiation  mechanisms  [10],  and  it  is  proposed  that  this  mechanism  also 
participates  in  polycrystalline  Ti.  The  significant  peak  for 
interdimple  spacirgs  in  Figure  4  is  consistently  at  or  near  50  pm  in 
larger  grain  size  material  and  gives  credence  to  a  glide  crossing 
mechanism  for  prime  void  initiation  sites  for  large-grained  Ti. 

Close  inspection  of  SEM  fractographs  shows  that  the  largest  number 
of  voids  grows  only  to  a  diameter  of  about  1  pm  as  documented  by  the 
highest  peaks  for  interdimplc  spacings  in  Figures  1  and  3.  Earlier  work 
on  the  initiation  of  microcracks  in  heavily  workhardened  metals  and 
alloys  [25,26,27,28,29]  pointed  to  the  importance  of  dislocation  cell 
walls  for  microvoid  initiation.  This  is  exactly  the  situation  in  the 
micro- 1 igaments  between  the  large  voids  discussed  so  far.  Here  e>1  on 
account  of  the  continuous  reduction  in  thickness  of  the  micro- 1 igaments 
to  a  fewr  micrometers.  A  detailed  description  of  microvoid  initiation  and 
growth  has  been  given  elsewhere  [30].  The  microdimples  under  discussion 
are  seen  primarily  on  the  rims  of  larger  dimples  and  indicate  areas 
where  localized  rupture  has  taken  place. 

A  determination  of  the  profile  of  a  dimpled  fracture  surface  can 
provide  information  as  to  the  precise  sequence  of  void  initiation, 
provided  the  resolution  of  the  technique  is  high  enough  to  discern  the 
relationships  between  dimples  for  the  two  opposite  fracture  surface 
profiles.  The  profiles  in  Figure  7  have  been  evaluated  for  a  time 
sequence  by  first  "closing"  the  fracture,  i.e.,  moving  the.  upper  and 
lower  profiles  over  each  other  so  that  no  white  areas  remain;  then 
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moving  the  profiles  apart,  one  can  monitor  the  initiation  sequence  (see 
Figure  12  for  the  appearance  of  initiation  of  voids  [1-3]).  Measuring 
the  wnole  sequence  of  16  initiation  sites  yields  an  average  distance 
between  sites  of  nearly  10  pm.  With  an  average  grain  size  of  13  ym  in 
this  specimen  the  claim  of  triple  points  as  primary  initiation  sites 
(see  Figure  1)  is  verified  within  the  error  limits.  This  determination 
was  made  under  the  assumption  that  void  growth  is  proportional  to  time, 
and  that  the  profile  would  be  going  through  the  center  of  the  dimples 
which  is,  of  course,  not  fulfilled.  A  more  careful  analysis  would 
involve  profiling  at  close  intervals  which,  if  done  for  a  larger  part  of 
the  fracture  surface,  would  contribute  to  the  quantitative  topographical 
characterization  of  the  time  sequence  of  void  initiation.  This 
measurement  is  adding  the  dimension  of  time  to  the  study  of  voids  in 
fracture  surfaces  at  microscopic  resolution  for  the  first  time. 

5 .  CONCLUSIONS 

(a)  In  specimens  up  to  grain  sizes  of  50  ym,  grain  boundary  triple 
points  are  primary  void  initiation  sites. 

(b)  When  grain  sizes  range  between  200  ym  and  1,300  ym  additional 
void  initiation  sites  are  provided  by  twin  intersections, 
glide  band  intersections,  dislocation  interactions  with  grain 
boundaries,  and  grain  boundary  ledges. 

ic)  Microvoids  arc  seen  at  the  rims  of  dimples  and  between 
dimples;  they  are  due  to  initiation  at  dislocation  cell  walls 
and  occur  during  rupture  at  the  final  separation. 

(d)  The  time  sequence  of  void  initiation  can  be  derived  from 
fracture  surface  profiles  which  were  constructed  from 
quantitative  topographical  information. 

( e )  It  has  been  demonstrated  that  quantitative  stereo-photo- 
graminetry  for  the  characterization  of  fracture  surfaces  is 
mandatory . 
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SECTION’  III 

THE  EFFECT  OF  THERMAL-MECHANICAL  COUPLING 
AT  THE  HEAD  OF  A  RAPIDLY  LOADED  CRACK  TIP 


Based  on  a  Ph.D.  Thesis  (Materials  Science')  by  J.  Daniel  Bryant,  UVA, 
1987* 

ABSTRACT 

An  analytical  method  has  been  developed  to  quantify  the  effect  of 
de format iona 1  heating  in  the  near  crack  tip  region  of  a  material 
subjected  to  rapid  loading.  The  finite  element  method  (FEM)  has  been 
employed  to  model  the  changes  in  the  stress  and  strain  fields  near  the 
crack  tip  due  to  the  effects  of  strain  induced  temperature  rises  for  a 
material  undergoing  elastic-plastic  behavior.  The  effects  of  these 
temperature  rises  are  manifested  in  both  altered  material  properties, 
i.e.  thermal  softening,  and  in  thermal  expansion.  The  dynamic 
interaction  between  the  evolving  stress,  strain,  and  temperature  fields 
incurred  during  rapid  loading  is  approximated  by  dividing  the  loading 
history  into  discreet  steps.  By  using  the  nodal  displacements  and 
elemental  temperatures  developed  in  the  previous  load  step  as  initial 
conditions  for  the  subsequent  step  the  model  emulates  the  true 
thermal -mechanical  coupling  experienced  during  shock  loading. 

Comparison  of  the  stress  and  strain  fields  of  the 
thermal -mechanica 1 ly  coupled  and  uncoupled  models  has  shown  the  effects 
of  temperature  rise  to  be  present  well  in  advance  of  the  crack  tip; 
significant  effects  include  diminished  maximum  stress  and  increased 
maximum  strain  in  the  coupled  model.  By  incorporating  the  results 
obtained  from  instrumented  test  specimens  used  in  the  determination  of 
loading  parameters,  insights  have  been  gained  into  the  effects  of 
nonadiabatic  conditions  and  the  app 1 i cab i 1 ity  of  the  method  to  the  case 
of  the  propagating  crack. 
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THE  EFFECT  OF  THERMAL-MECHANICAL  COUPLING  AT  THE 
HEAD  OF  A  RAPIDLY  LOADED  CRACK  TIP 

1.  INTRODUCTION 

Deformational  heating  has  long  been  acknowledged  as  having  a 
significant  influence  on  the  process  of  dynamic  fracture.  Early  works 
by  Zener  and  Hollomon  describing  the  effect  of  strain  rate  on  plastic 
flow  in  steel  elucidate  the  opposing  roles  played  by  the  well  known 
phenomenon  of  strain  rate  hardening  and  that  of  thermal  softening  as  a 
result  of  adiabatic  heating.  More  recent  works  have  concentrated  on  the 
effects  of  the  instability  reached  when  adiabatic  softening  dominates 
the  fracture  process.  The  formation  of  adiabatic  shear  bands  in  metals 
discussed  by  Rogers  in  his  review  and  by  Stelly  et  al  and  others  in  the 
field  is  associated  with  the  instability  dictated  by  the  slope  of  the 
adiabatic  stress -strain  curve.  Evidence  of  rapid  temperature  rises 
during  adiabatic  heating  has  presented  itself  in  the  form  of 


transformation  bands  and  incipient  melting  of  the  fracture  surface. 
Figure  1  shows  such  a  heat  affected  area  on  the  fracture  surface  of  a 
beta  titanium  alloy  subjected  to  high  strain  rates. 

A  detailed  numerical  model  of  the  propagation  of  adiabatic  shear 
bands  was  recently  accomplished  by  Kuriyama  and  Meyers  [1].  In  their 
analysis,  the  region  ahead  of  a  crack  is  modeled  for  a  material  obeying 
an  adiabatic  stress -strain  relation  and  subjected  to  pure  shear  (Mode 
II)  stresses.  In  concentrating  on  the,  near  crack  tip  region,  the  model 
is  thus  analogous  to  the  approach  of  fracture  mechanics.  The  case  of  a 
crack  tip  experiencing  a  Mode  I  type  opening  stress  was  addressed  by 
Hoff,  Rubin  and  Hahn  [2].  A  compact  tension  specimen  subjected  to 
dynamic  loading  was  modeled  using  finite  element  analysis  (FEM) ;  the 
results  describe  the  changes  in  the  stress  and  strain  fields  produced  by 
strain  rate  sensitive  and  insensitive  materials.  The  model,  however, 
did  not  account  for  the  effects  of  adiabatic  heating  in  the  rapidly 
strained  region  at  the  crack  tip. 

In  the  present  study,  the  stress  and  strain  fields  at  the  crack  tip 
in  a  compact  tension  specimen  are  calculated  for  the  case  of  deformation 
dominated  by  adiabatic  heating  etiecLs.  To  address  the  problems  of 
strain  rate  sensitivity,  an  experimental  testing  program  has  been  run  in 
parallel  and  thus  acts  to  supply  the  parameters  needed  to  describe  both 
the  material  response  to  accelerated  deformation  rates  and  the  loading 
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history  for  the  specimens  being  modeled.  Two  manifestations  of 
adiabatic  heating  on  the  stress-strain  distribution  are  considered: 
thermal  softening  and  thermal  expansion.  The  effects  of  the  dynamic 
interaction  of  the  stress,  strain  and  temperature  fields  are  determined 
for  an  elasto-plastic  temperature  dependent  material  using  the  finite 
element  method.  In  this  approach,  the  thermal-mechanical  coupling 
relating  these  three  fields  is  determined  by  the  modulation  of  the 
elemental  temperatures  at  discreet  intervals,  these  temperatures  being 
determined  by  heat  generated  through  plastic  deformation.  By  studying  a 
rapidly  loaded  stationary  crack  undergoing  adiabatic  heating,  the  stress 
and  temperature  fields  derived  can  be  used  to  predict  the  response  of 
more  complicated  systems,  including  the  nonadiabatic  case  and  that  of 
the  propagating  crack. 

2.  THEORETICAL  AND  EXPERIMENTAL  PROCEDURE 
(i)  Experimental  Method 

Specimens  for  tensile  and  KT  tests  were  machined  from  a 

Ic 

billet  of  Ti-10V-2Fe-3Al  of  the  composition  given  below. 

Table  1:  Chemical  Composition  (weight  percent) 

C  N  Fe  V  0  A1  Y  Ti 

.03  .013  1.8  9.9  .095  3.4  10  ppm  balance 


Specimens  were  homogenized  above  the  beta  transus  temperature  at  825°C. 
This  was  followed  by  a  furnace  cool  to  780°C  and  a  quench  to  room 
temperature  water.  A  subsequent  aging  treatment  was  performed  at  500°C 
for  one  hour. 

Fracture  toughness  tests  were  made  through  testing  of  compact 
tension  and  three  point  bend  specimens.  Testing  was  performed  on  a 
closed  loop  hydraulic  testing  unit.  As  the  loading  rates  used  in  these 
experiments  were  far  higher  than  could  be  recorded  using  pen  and  chart, 
a  high  speed  data  acquisition  system  was  developed.  The  load  signal  was 
recorded  using  an  A'lT  6300  computer  at  intervals  determined  by  the 
on-board  clock.  To  determine  the  exact  moment  of  crack  propagation  (and 
hence  the  load  at  the  onset  ol  crack  motion)  a  strain  gage  was  cemented 
directly  ahead  of  the  pre-crack  and  the  signal  monitored  using  the  same 
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could  be  determined 


equipment.  In  this  way,  the  fracture  toughness 

#  • 

as  well  as  the  loading  rate  parameter  K  where  K  =  K  /t  and  t  is  the 

1  -L  i. 0  C  C 

critical  time  to  reach  crack  growth.  Loading  rate  was  varied  from  1.5 
Mpa/ms  1  to  1.28xl04  MPa/ms  ^ .  The  data  recorded  provided  the  required 
loading  history  used  in  the  finite  element  model  developed.  The  same 
data  aquisition  system  was  used  to  measure  crack  propagation  speed; 
using  a  crack  propagation  gage  (an  array  of  parallel  conductors)  the 
resistance  drop  was  measured  as  the  advancing  crack  severed  the  gage. 
The  instrumented  compact  tension  specimens,  as  well  as  instrumented 
tensile  specimens  used  in  high  strain  rate  tests,  are  shown  in  Figure  2. 

(ii)  Modeling  Method 

The  finite  element  method  (FEM) ,  a  system  of  mathematical 
modeling  which,  for  static  stress  problems,  determines  the  displacements 
of  the  nodes  contained  within  the  element  mesh  in  response  to  forces  and 
displacements  impressed  upon  that  mesh  as  boundary  conditions.  The  mesh 
is  constructed  by  the  user  to  emulate  the  geometry  of  the  system  to  be 
analyzed.  In  this  analysis,  the  stress  and  strain  fields  near  the  crack 
tip  in  a  compact  tension  (CT)  specimen  were  modeled  with  the  aid  of  the 
ANSYS  finite  element  codes  [3].  This  mesh  configuration  was  chosen  for 
two  reasons.  First,  the  stress  field  near  the  crack  tip  of  such  a 
specimen  has  been  extensively  studied;  this  allows  for  comparison  of  the 
current  model  results  with  previous  work.  Secondly,  by  performing 
mechanical  tests  upon  specimens  of  this  same  configuration,  the 
parameters  of  the  model,  namely  the  loading  history  and  the  time 
indicating  the  onset  of  crack  growth,  can  be  determined  directly  from 
experimentally  measured  quantities. 

In  the  present  work,  a  four  noded  isoparametric  stress  solid  was 
chosen  as  the  element  type.  This  two-dimensional  model  was  made  to 
emulate  a  three-dimensional  model  of  unit  thickness  undergoing  plane 
strain.  A  diagram  of  the  element  mesh  is  given  in  Figure  3.  The  load 
is  applied  at  the  node  marked  in  the  upper  left  portion  of  the  diagram; 
this  corresponds  to  the  center  point  of  the  load  pin  in  the  actual  CT 
specimens.  As  shown,  the  direction  of  crack  propagation  is  to  the 
right.  As  is  standard  in  finite  element  analysis,  the  total  specimen 
being  modeled  is  divided  along  its  lines  of  symmetry  and  a  single 
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The  element  mesh  used  to  model  the  rapid  loading  of  a 
compact  tension  specimen.  The  load  is  applied  at  the 
node  marked,  and  the  crack  is  to  propagate  to  the  right. 
Only  the  upper  half  of  the  CT  specimen  is  modeled,  due  to 
svmmet  ry . 
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section  analysed.  The  line  of  symmetry  in  this  problem  is  along  the 
line  Y=0 ,  which  includes  the  crack  tip  itself.  To  preserve  the 
equilibrium  of  the  mesh  in  response  to  the  loading  force,  displacement 
constraints  are  placed  upon  the  series  of  nodes  along  the  X  axis  to  the 
right  of  the  crack  tip.  The  mesh  generated  consists  of  144  nodes 
forming  125  elements.  In  the  region  close  to  the  crack  tip,  the  density 
of  the  elements  increases,  as  it  is  here  that  the  greatest  stress 
gradients  will  occur.  An  enlargement  of  the  near  crack  tip  region  is 
shown  in  Figure  4,  with  the  nodal  numbers  marked.  Concentric 

semicircles  surround  the  crack  tip  itself  (node  1)  providing  an 
elemental  spacing  of  one  micron  directly  next  to  the  crack  tip,  which  is 
represented  in  this  model  as  sharp  and  unblunted. 

The  compliance  of  the  element  mesh  in  response  to  an  external  force 
or  imposed  displacement  is  translated  into  nodal  displacements  via  the 
stiffness  matrix.  The  stiffness  matrix  incorporates  the  material 
properties;  for  the  elasto-plastic  model  developed  here  a  bi-lincar 
stress-strain  curve  W'as  used  to  describe  the  alloy's  mechanical 
properties,  as  shown  in  Figure  5.  The  temperature  dependence  of  the 
elastic  and  plastic  portions  of  the  lower  stress-strain  curves  in  Figure 
5  is  approximate  and  based  upon  empirical  data  taken  from  another 

titanium  alloy  for  which  high  temperature  data  was  available,  namely 
Ti-oAl-4Y. 

A  single  nodal  force,  situated  at  node  number  82,  acts  in  the 
vertical  direction,  thereby  simulating  the  action  of  t!he  load  pin  in  the 
actual  specimen.  The  maximum  load  applied  at  this  node  was  determined 
from  one  of  the  K  tests  run  at  a  loading  rate  near  the  center  of  the 
range  studied.  In  both  the  coupled  and  the  uncoupled  model,  this  load 
is  applied  in  discreet  load  steps;  the  maximum  load  is  applied  over  a 
series  of  10  steps.  There  are  two  reasons  for  this.  Firstly,  by  its 
very  nature,  a  problem  incorporating  plasticity  is  nonlinear;  the 
solution,  therefore,  can  only  be  arrived  at  through  a  process  of 
iteration.  By  dividing  the  total  load  into  a  series  of  steps,  the 

problem  ol  nonconve l'gence  may  bn  avoided.  The  second  reason  lies  in  the 

method  used  in  the  thermal -median ical ly  coupled  model  to  produce  a 
dynamic  interaction  between  the  elemental  temperature  and  the  stress  and 
striin  fields.  Incrementing  the  load  to  the  maximum  load  over  ten  steps 
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The  near  crack  tip  region  of  the  element  mesh.  The 
Tensity  of  the  elements  increases,  as  the  stress 
gradients  near  the  crack  tip  will  be  higher. 
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allows  for  the  periodic  adjustment  of  the  elemental  temperatures,  as 
determined  by  the  plastic  strain  energy  expended  in  the  previous  load 
stop . 

The  number  of  iterations  required  to  reach  a  particular  margin  of 
error  within  a  load  step  is  controlled  in  the  program  through  a 
convergence  criterion.  For  a  stress  solid  element  type,  the  convergence 
criterion  is  defined  by  the  ratio: 

C=Ae  ft 
P  e 

where  z  is  the  elastic  strain  and  At  is  the  change  in  the  plastic 

e  p 

strain  with  the  previous  iteration.  Obviously,  the  smaller  the 

convergence  criterion  set,  the  greater  precision  with  which  the  model 
will  simulate  a  material  described  by  a  given  elasto-plastic 

stress-strain  curve.  The  disadvantage  lies  in  the  greater  computational 
time  and  expense  for  this  added  precision.  A  convergence  criterion  of 
C=.05  was  chosen  in  both  models.  Depending  on  the  load  step,  between  5 
and  20  iterations  were  required  to  match  this  criterion. 

As  stated  in  the  introduction,  the  deformation  induced  heating  in 
an  element  near  the  crack  tip  acts  to  change  the  defined  temperature  of 
the  element  involved,  and  this  in  turn  affects  the  relationship  of  the 
stress  and  strain.  At  this  point  the  method  used  to  determine  the 
temperature  rise  and  the  effects  of  that  temperature  rise  will  be 
d i s  cussed . 

While  a  formal,  analytical  expression  relates  the  stress  and  strain 
in  an  elastic  material  to  the  stored  energy  per  unit  volume,  this 
expression  is  not  exact  for  the  case  of  an  elasto-plastic  material 
because  of  the  path,  dependence  of  the  integral.  In  a  plastically 
deforming  material,  the  stored  energy  per  unit  volume  is  approximated  by 
the  expression: 

Ah  =  ("ode 

E 

w ! ;  »*  r  *  *  z  md  £  -  are  the  vie  Id  strain  and  fracture  strain  respectively, 
y  1  ! 

The  issoc  i  at  e,l  temperature  rise  within  this  element  is  given  by: 
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where  C  is  the  specific  heat,  p  is  the  density  and  R  is  the  proportion 
of  energy  converted  directly  to  heat.  From  previous  works,  R  is 
generally  accepted  to  have  a  value  between  .90  and  .95.  This  equation 
cannot  be  used  analytically  to  describe  the  temperature  rise  for  a 
material  undergoing  plastic  deformation  given  the  path-dependent  nature 
of  plastic  deformation.  Because  of  this,  only  an  approximate  relation 
can  be  made  between  the  stress  and  plastic  strain  acting  in  an  element 
and  the  energy  expended  per  unit  volume.  The  actual  amount  of  energy 
expended  is  of  course  dependent  upon  the  size  of  the  element;  a  larger 
amount  of  energy  will  be  expended  in  deforming  a  larger  element.  This 
problem  of  determining  elemental  volume  factually,  elemental  area  with  a 
given  unit  thickness)  can  be  avoided  by  directly  determining  the 
temperature  rise  within  each  element.  As  the  density  and  specific  heat 
are  both  on  a  per  unit  basis,  the  volume  of  a  given  element  is  not 
required  for  establishing  the  elemental  temperature.  The  plot  shown  in 
Figure  6  shows  graphically  the  method  used  to  determine  the  elemental 
temperature  rise.  An  approximation  is  made  that  the  area  beneath  the 
stress-plastic  strain  curve  is  described  by  a  rectangle  with  an  area 
equal  to  the  product  of  the  stress  and  plastic  strain.  The.  area  beneath 
this  curve,  as  well  as  the  similar  curves  for  stress  and  plastic  strain 
acting  in  the  Y  and  shear  directions,  is  equated  to  the  energy  per  unit 
volume  expended  in  the  action  of  plastic  flow.  The  associated 


temperature  rise  within  an  element  is  given  by  the  equation: 


!  ( o  e  )  +  (o  e  )  +  2‘  a  £  ) 

A1  =  - 1  x  px  v  py  xv  pxy 

pc 

PI 


where  the  absolute  value  functions  have  been  employed  to  account  for  the 


positive,  scalar  nature  of  the  heat  generation. 


The  elemental  temperatures  lor  those  elements  inside  of  the 
e  1  as to-p 1  as t i c  boundary  are  defined  using  the  above  equation  at  10 
points  during  the  loading  historv,  i . e . ,  the  final  iteration  of  the  10 
load  steps.  The  subsequent  load  step  is  affected  by  this  modulation  of 
elemental  temperatures  in  two  ways.  First,  the  appropriate 
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stress-strain  curve  will  be  selected  for  an  element  at  that  defined 
temperature.  Secondly,  the  thermal  expansion  of  that  element  will  be 
incorporated,  and  this  will  be  reflected  in  the  stress  and  strain  values 
for  that  element  and  those  neighboring  it. 

(a)  The  Uncoupled  Model.  In  the  uncoupled  model,  the  elemen¬ 
tal  temperatures  are  maintained  at  a  constant,  ambient  temperature  of 
300°C  throughout  the  loading  history.  This  temperature  was  set  in 
the  pre-processing  stage  of  the  analysis,  and  as  a  result  the  material 
properties  relating  the  imposed  loads  to  the  nodal  displacements  are 
prescribed  by  the  upper  most  curve  in  Figure  5. 

(b)  The  Coupled  Model.  In  the  coupled  model,  the  requirement 
imposed  by  the  modulation  of  the  elemental  temperatures  at  the  end  of 
each  load  step  necessitates  a  different  approach  than  that  used  for  the 
uncoupled  model.  The  first  load  step  applied  is  identical  to  that  of 
the  uncoupled  model,  and  in  the  subsequent  load  steps  the  load  is 
applied  to  element  mesh  using  the  nodal  displacements  and  stresses  from 
the  previous  step  as  initial  conditions.  In  the  coupled  model  the 
elemental  temperatures  are  adjusted  for  each  element  subjected  to 
plastic  strain  in  the  previous  load  step.  This  process  is  repeated  at 
the  end  of  each  load  step.  With  the  iterations  required  to  reach  a 
convergence  criterion  of  .05,  a  total  number  of  192  iterations  with  a 
running  time  of  approximately  eight  hours  was  required. 

3.  RESUbTS  AMD  INTERPRETATION 

1 

The  objective  of  performing  this  analysis  was  to  gain  insight  into 
the  role  which  deformation  induced  heating  plays  in  dynamic  fracture. 
Specifically,  there  are  two  subjects  which  are  to  be  addressed.  The 
effect  of  thermal -mechanical  coupling  at  the  crack  tip  on  the  stress  and 
strain  fields  is  to  be  discussed.  In  addition,  an  associated 
temperature  field  is  predicted  under  the  assumption  of  adiabatic 
deformut iona 1  heating.  In  the  model,  adiabatic  boundaries  are  imposed 
between  elements;  no  transfer  or  heat  is  permitted  away  from  the  crack 
tip  area.  The  adiabatic  temperature  field  predicted  by  the  model  can 
tleui  be  inulyxed  from  the  standpoint  of  transient  heat  generation  and 
heat  t  rails  f " r  . 

The  nodal  displacements  determined  from  the  analysis  are  used  to 
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predict  the  elements  of  the  stress  and  strain  matrices  at  each  node.  In 
addition,  the  stress  and  strain  values  of  the  four  nodes  defining  the 
perimeter  of  an  element  can  be  averaged  to  yield  the  elemental,  or 
centroidal,  stresses  and  strains.  The  stress  and  strain  data,  as  well 
as  temperature  profiles,  are  graphically  presented  in  the  form  of 
connected  contour  plots  determined  from  centroidal  values. 

(i)  Analysis  of  Stress  and  Strain  Fields 

In  the  first  set  of  plots.  Figures  7  and  8,  the  plastic  strain 
in  the  Y.  or  vertical,  direction  is  compared  for  coupled  and  uncoupled 
models.  The  higher  values  of  plastic  strain  in  the  thermal-mechanical  ly 
coupled  model  are  reflected  in  the  difference  in  maximum  plastic  strain: 
iy.5°0  for  the  uncoupled  model  as  compared  to  26. 6°,,  for  the  coupled 
model.  The  higher  plastic  strain  levels  reached  in  the  coupled  model 
are  a  consequence  of  the  thermal  softening  in  the  near  crack  tip  region. 
The  perimeter  of  the  plastically  deformed  region,  or  the  elasto-plastic 
boundary,  can  be  interpolated  from  the  stress-data  printout.  For  both 
models,  plasticity  is  seen  only  to  occur  within  the  fifth  semi-circular 
regime  of  elements;  this  defines  an  approximate  plastic  zone  radius  of 
42  micrometers.  This  is  in  good  general  agreement  with  the  plastic  zone 
radius  defined  by  Pugdale  for  the  plane  strain  case,  given  by: 


1  1 

(k* 

Ic 

6tt  ' 

\°2 

Using  experimentally  determined  values  for  the  yield  stress  and  critical 
stress  intensity,  this  is  calculated  at  56  micrometers. 

The  shape  of  the  contour  plots  of  von  Mises'  stress  for  the  coupled 
and  uncoupled  models  are  seen  to  be  roughly  the  same  (Figures  9  and  10). 
A  lower  stress  in  the  near  crack  tip  region,  as  well  as  a  lower  maximum 
stress,  can  be  seen  in  the  coupled  model.  This  is  a  result  of  not  only 
the  thermal  softening,  which  lowers  the  stress  bearing  capacity  of  the 
element,  but  also  of  the  thermal  expansion.  In  the  elements  under 
positive  stress,  the  thermal  expansion  caused  by  the  increased  elemental 
temperature  acts  to  reduce  the  stress  by  providing  additional  compliance, 
to  the  imposed  forces  at  the  nodes. 
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igure  7.  The  contour  plot  showing  the  plastic  strain  in  the  f,  or 

vertical,  direction  for  the  uncoupled  model.  The  data 
column  on  the  right  lists  the  maximum  strain  (MX} ,  the 
minimum  strain  (MN) ,  the  number  of  contours  (NCON)  and 
the  strain  increment  between  contours  (VINC). 


,CTS -PLASTIC 


ur<>  ') .  The  contour  plot  of  the  van  Mises  stress  near  the  crack 

tip  for  the  uncoupled  t'EM  model. 
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The  most  interesting  result  of  the  analysis  can  be  seen  in  the 
stress  fields  in  the  Y,  or  vertical  direction,  as  it  most  clearly 
depicts  effects  of  Mode  I  type  loading.  In  Figure  11,  the  stress  field 
for  the  uncoupled  model  is  shown.  This  can  be  compared  to  the  stress 
field  for  the  thermal-mechanically  coupled  model  in  Figure  12,  and 
several  important  differences  can  be  noted.  The  intensity  of  the 
gradient  in  the  stress  directly  ahead  of  the  crack  tip  is  greatly 
reduced  in  the  model  which  incorporates  thermal  softening.  The  maximum 
value  of  stress  in  the  Y  direction  (in  the  units  used  in  the  model  of 
Kgf/mm)  in  the  uncoupled  model  is  seen  to  be  553  kgf/mm,  while  in  the 
coupled  model  the  maximum  stress  is  reduced  35  percent  to  a  value  of  360 
kf g/mm . 

(ii)  Transient  Thermal  Analysis 

Using  the  thermal-mechanically  coupled  finite  element  analysis 
done  for  a  compact  tension  specimen,  the  temperature  fields  generated 
through  de format iona 1  heating  can  now  be  subjected  to  an  analysis  of  the 
thermal  flow  as  a  function  of  position  and  time.  The  model  described  in 
the  previous  section  is  specifically  for  the  case  of  a  rapidly  loaded 
compact  tension  specimen  where  the  heat  of  plastic  flow  is  allowed  to 
perturb  the  material  properties  in  that  element  alone.  No  attempt  has 
been  made  to  consider  nonadiabatic  conditions;  the  heat  generated  in 
elements  undergoing  plastic  deformation  was  not  allowed  to  dissipate  to 
the  neighboring  elements.  This  was  an  unavoidable  consequence  of  the 
modeling  method  used.  Rather  than  running  a  thermal  finite  element 
analysis  using  simply  the  temperature  field  predictions  from  the  final 
iteration,  a  one  dimensional  thermal  analysis  was  performed  based  upon 
heat  generation  rate.  It  is  the  author's  belief  that  this  method  more 
clearly  elucidates  the  role,  of  heat  conduction  in  the  near  crack  tip 
reg ion  . 

(iii)  Predicted  Temperature  Field 

The  temperature  field  predicted  through  the  thermal- 
mechanically  coupled  stress  model  is  shown  in  Figure  13.  In  the  right 
hand  column  next  to  the  diagram  is  the  information  required  to  interpret 
the  diagram.  The  colored  contours  are  at  levels  of  lil°K,  and  the 
ambient  temperature  of  the  specimen  was  defined  at  300°K  (this  is  of 


The  Mode  I  type  stress  plotted  tor  the  uncoupled  model  in 
Figure  11  can  be  compared  with  the  equivalent  plot  for 
the  thorma 1 -median ica 1 ly  coupled  FKM  model.  Important 
points  to  note  are  the  shape  of  the  stress  contours,  and 
the  maximum  stress  near  the  crack  tip,  which  has  been 

_  O 

reduced  to  3r>0  Kgf-mm  “  on  account  of  the  thermal 
softening  and  thermal  expansion  due  to  deformation 
induced  heating. 


The  temperature  field  predicted  in  the.  thermal  - 
mechanically  coupled  finite  element  model.  Temperature 
values  are  in  degrees  Kelvin,  with  a  maximum  predicted 
temperature  increase  of  852°C. 


course  the.  minimum  temperature  value,  as  shown  in  the  right  hand  data 
column).  The  maximum  temperature  of  879°C  (1152°K)  occurs  near  the 
boundary  between  the  first  and  second  regime  of  semicircular  element 
arrays;  this  places  the  position  of  maximum  temperature  rise  at  a 
position  approximately  1.5  microns  from  the  crack  tip.  The  split  in  the 
stiape  of  the  temperature  field,  which  can  be  seen  near  the  vertical 
position,  is  a  consequence  of  the  shape  of  the  field  describing  the 
plastic  shear  strain  near  the  crack  tip.  The  plastic  shear  strain 
changes  sign  near  this  region.  This  results  in  a  smaller  heating 
contribution  due  to  plastic  flow  under  this  component  of  stress. 

It  should  be  noted  that  this  model  does  not  incorporate  the  very 
large  strains  present  in  an  actual  crack  tip  within  the  process  zone. 
The  process  zone,  which  may  have  a  radius  of  less  than  a  micrometer,  has 
been  described  as  a  region  within  the  plastic  zone  characterized  by  very 
large  strains,  which  may  be  more  than  500  percent.  The  finite  element 
analysis  here  was  restricted  to  plastic  strains  of  less  than  35  percent. 

It  is  logical,  therefore,  that  the  analysis  of  the  heating  effects 

within  the  process  zone  developed  in  a  previous  work  are  applicable  to 
this  case  and  may  be  superimposed  upon  those  predicted  by  the  FF.M . 

(iv)  One-Dimensional  Trans ient  Therma 1  Analysis 

The  temperature  field  predicted  by  the  therma 1 -mechanical ly  coupled 
stress  analysis  provides  the  spatial  distribution  of  temperature  on 
which  the  thermal  analysis  may  be  performed.  The  experimentally 
determined  loading  history  and  critical  time  (the  time  from  initial 

loaning  to  initial  crack  growth)  gives  the  times  over  which  the  thermal 
model  will  be  applied.  The  analysis  done  here  is  based  on  a  treatment 
originally  proposed  by  Cars  law  and  Jaeger. 

Heat  generation  near  the  crack  Lip  can  best  be  modeled  in  one 

dimension  as  a  semi  -  inf  inite  plate.  The  near  crack  tip  region  in  which 
heat  is  being  generated  is  mode1  led  as  a  band,  of  width  I,,  located  at  the. 
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plastic  strain.  The  equation  describing  the  temperature  field  as  a 
function  of  time  and  position  for  any  point  X  within  the  heat  generating 
zone  is  given: 


where  A  defines  a  constant  heat  generation  rate,  K  the  coefficient  of 
o 

thermal  conduciton,  and  k  the  thermal  diffusivity.  For  any  point  X 
outside  of  the  heat  generating  zone,  the  temperature  is  described  by: 

i 

!  X  >  L 

Obviously,  the  condition  cf  uniform  heat  generation  rate  is  approximate, 

and  only  valid  for  the  case  of  a  stationary  crack  tip  at  times  less  than 

t  .  At  times  greater  than  t  ,  the  crack  has  propagated  away  from  this 
c  c 

region  and  those  heat  generation  sites  active  during  loading  are  now 
quiescent.  If  the  value  of  L,  the  width  of  the  heat  generating  zone,  is 
estimated  from  Figure  13,  the  equations  can  be  applied  to  the  case  of 
the  stationary  crack  tip  undergoing  rapid  loading.  Using  the  physical 
properties  of  Ti-10-2-3  and  the  relation  of  the  thermal  diffusivity  to 
the  other  parameters  given  by: 

K 

r  = - 

PCP 

the  first  equation  may  be  rewritten  in  the  form: 

1  r 

AT  = -  1  -  2  i  2er  f  c 

PCP 

The  temperatures  predicted  under  adiabatic,  restrictions  are  simply  the 
values  dot  iued  by  the  first  coetficient  of  the  equation.  F.xprossed  as  a 
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The.  effect  of  thermal  conduction  on  the  temperature  field 
predicted  by  the,  finite  element  analysis  for 
T  i  - 1 OV- 2Fe - 3A 1  can  be  judged  vis  ing  the  curves  shown 
above.  The.  ratio  of  the  actual  temperature'  to  the 
adiabatic  prediction  is  read  off  of  the  ordinate.  The 
effective  loading  time,  is  read  off  of  the  abscissa.  !., 
the  width  of  the  heat  generating  zone,  was  taken  it. 
values  of  i tul  10  micrometers,  from  Figure  13. 


The  impact  of  thermal  conduction  away  from  the  heat  generating  region  is 
controlled  by  the  arguments  of  i2erfc  functions.  As  the  right  hand  side 
of  the  equation  approaches  unity  (the  second  and  third  terms  on  the  RHS 
approach  zero),  the  temperature  field  approaches  the  adiabatic  case 
predicted  by  the  finite  element  model.  In  Figure  14,  this  function  is 
plotted  for  several  values  of  L  and  X  estimated  from  the  spatial 
distributions  in  Figure  13. 

Several  conclusions  can  be  made  about  the  shape  of  these  curves. 
First,  at  very  short  times,  on  the  order  of  1  to  10  microseconds,  there 
is  very  little  effect  of  thermal  diffusion  in  the  near  crack  tip  region. 
At  larger  times,  the  effect  of  thermal  conduction  becomes  great, 
lowering  the  temperature  to  nearly  the  ambient  level  in  a  time  of  a 
millisecond.  This  is  a  consequence  of  the  very  narrow  width  of  the  heat 
generating  region. 

When  the  critical  times  taken  from  the  test  results  using  compact 
tension  specimens  are  applied  to  the  equation  above,  it  is  apparent  that 
thermal  conduction  plays  a  heavy  role  in  the  determination  of  the 
temperatures  reached  near  the  crack  tip.  The  smallest  experimentally 
determined  critical  time  t  of  2.2  milliseconds  places  the  temperature 
within  the  heat  generating  zone  at  a  value  much  closer  to  the  ambient 
temperature  than  that  predicted  by  the  adiabatic  model.  The  thermal 
transfer  model,  together  with  the  experimental  data,  leads  one  to 
conclude  that  the  effect  of  deformational  heating  on  a  rapidly  loaded 
stationary  crack  tip  is  slight. 

This  result  is  not  surprising  considering  the  relatively  large 
amount  of  time  required  to  load  a  specimen  to  critical  stress  intensity 
for  crack  growth.  If  the  model,  though  constructed  for  the  case  of  a 
stationary  crack,  is  applied  to  the  case  of  a  propagating  crack,  quite  a 
different  conclusion  is  reached.  If  one  may  safely  ignore  the  effects 
of  inertia,  as  a  recent  finite  element  work  has  shown,  the  finite 
element  model  developed  here  may  be  applied  to  the  case  of  a  moving 
crack  tip.  The  parameter  describing  the  time  scale  of  loading  may  now 
lie  deduced  from  the  experimentally  measured  crack  propagation  speeds  and 
the  .tress  calculations  discussed  in  the  previous  section. 

Fruit  Figure  11,  the  'lode  1  type  stress  drops  off  i  rom  its  peak 
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value  near  the  crack  tip,  to  a  value  roughly  half  of  this  at  a  distance 
20  micrometers  ahead  of  the  crack  tip.  For  a  point  ahead  of  the 
propagating  crack,  the  increasing  stress  at  that  point  as  the  crack 
advances  closer  is  equivalent  to  a  loading  rate  directly  proportional  to 
the  speed  of  the  advancing  crack.  Using  the  distance  over  which  the 
imposed  stress  drops  to  half  of  its  peak  value  and  the  measured  crack 
velocity,  one  calculates  equivalent  loading  times  for  points  ahead  of  a 
propagating  crack.  For  a  crack  speed  of  10  ms  \  the  argument  within 
square  brackets  in  equation  6  takes  a  value  of  .95,  or  95  percent  of  the 
temperature  rise  value  predicted  by  the  adiabatic  model. 

4.  CONCLUSIONS 

A  finite  element  model  has  been  developed  and  used  to  quantify  the 
dynamic  interaction  between  the  stress,  strain  and  temperature  fields 
near  a  rapidly  loaded  crack  tip  undergoing  deformational  heating.  The 
model  has  predicted  the  following: 

1.  The  action  of  deformation  heating  acts  to  increase  the  maximum 
plastic  strain  within  the  plastic  zone  by  36  percent. 

2.  The  von  Mises  and  Mode  1  type  stress  are  diminished  between  30  and 
35  percent  in  the  therma 1-mechanical ly  coupled  model.  The  lowered 
load  bearing  capacity  caused  by  thermal  softening  is  assisted  by 
thermal  expansion  in  the  near  crack  tip  elements. 

3.  The  temperature  field  produced  through  deformational  heating  has 

been  determined  for  the  adiabatic  case.  Temperature  increases  of 

up  to  830°C  are  predicted  for  some  regions  within  the  plastic  zone 
(excluding  the  process  zone). 

4.  The  effect  of  non-adiabatic  conditions  has  been  gaged  using  the 

temperature  field  calculated  from  the  FEM  and  the  equations  for 
heat  transfer.  It  has  been  determined  that  only  under  exceedingly 
high  loading  rates  (10^  MPa/m/s)  will  adiabatic  type  heating  play  a 
predominant  role  on  the  stress  and  strain  fields  of  a  stationary 
nr  ink  . 

5.  In  the  case  of  the  propagating  crack,  it  was  found  that  the 

effective  loading  rate  it  which  adiabatic  heating  effects  become 
significant  occurs  when  crack  propagation  speeds  reach  the  realm  of 
In  m/s  and  over. 

6  5 
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Based  on  a  Pli.D.  Thesis  (Materials  Science)  by  David  D.  Makel,  UVA,  1987 

ABSTRACT 

During  the  examination  of  the  fracture  surfaces  of  primarily  plane 
stress  TiSMn  tensile  samples,  regions  were  found  which  have  a  finely 
roughened  texture  and  are  covered  with  small  spheroidized  volumes  of 
material.  The  high  strength  and  low  thermal  conductivity  of  this 
titanium  alloy  suggest  that  these  unusual  features  may  result  from 
adiabatic  type  heating  effects  during  the  final  separation. 

To  investigate  the  possibility  of  high  temperatures  being  generated 
during  straining  to  fracture,  tensile  samples  were  tested  over  a  range 
of  strain  rates  and  the  resulting  fracture  surfaces  were  thoroughly 
characterized  with  stereo  pairs  of  scanning  electron  micrographs.  After 
characterization,  the  surfaces  were  nickel  plated,  imbedded  in  resin  and 
sectioned  to  reveal  the  microstructures  underlying  the  various  surface 
structures.  These  sections  reveal  highly  localized  zones  of  shear 
directly  below  the  unusual  surface  features. 

Localized  shear  and  high  temperatures  are  two  of  the  basic 
components  needed  for  adiabatic  shear,  a  deformation  mechanism  generally 
associated  with  considerable  compressive  stresses  and,  therefore,  not 
usually  considered  as  playing  a  role  in  tensile  testing.  To  determine 
whether  or  not  the  unusual  surface  features  are  caused  by  local  melting 
due  to  adiabatic  shear,  the  effects  of  high  surface  temperatures  were 
investigated,  as  well  as  the  possible  role  of  oxidation.  A  geometrical 
model  of  adiabatic  shear  is  presented  in  this  report. 

Analyses  of  the  heating  effects  which  occur  during  straining  and 
separation  indicate  that  temperatures  as  high  as  the  melting  temperature 
are  reached  at  separation,  and  that  elevated  temperatures  play  an 
important  role  in  the  material  instability  and  strain  localization  which 
lead  to  fracture. 
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A  GEOMETRICAL  MODEL  OF  THE  STRAIN  RATE 
DEPENDENT  FRACTURE  OF  Ti-8Mn 


1 .  INTRODUCTION 

The  primary  thrust  of  this  study  has  focused  on  the  effects  of 
strain  rate  on  the  strength,  toughness  and  fracture  mode  of  high 
strength,  low  heat  conductivity  materials,  in  this  case,  titanium  8  w°0 
manganese.  A  synopsis  of  the  results  of  this  research  can  be  broken 
down  into  three  major  categories: 

(1)  Observation  of  apparent  high  temperature  fracture  surface 
features  ; 

(2)  Identification  and  verification  of  adiabatic  shear  as  the 
primary  mechanism  responsible  for  the  temperature  rise;  and 

(3)  The  application  of  the  concept  of  adiabatic  shear  to  the 
overall  fracture  sequence. 

A  brief  description  of  the  experimental  methods  used  in  our  studies 
will  be  given  followed  by  the  significant  results  and  a  discussion  of 
the  implications  and  conclusions  which  can  be  drawn. 

2.  INITIAL  TESTING  AND  OBSERVATIONS 

(il  Tensile  Testing 

The  primary  alloy  used  in  this  research,  titanium  8  w°„ 

manganese,  was  chosen  because  it  is  easily  heat  treated  to  a  condition 
of  high  strength  and  moderate  toughness,  it  has  a  very  low  thermal 
conductivity  thereby  enhancing  any  adiabatic-tvpe  effects,  and  it  is  a 
binary  alloy  with  a  well-defined  phase  diagram  and  thoroughly 
investigated  microstructures.  The  use  of  a  standard  commercial  alloy 
eliminated  the  necessity  of  exhaustive  microstructural  classifications 
prior  to  the  mechanical  testing  and  post - f racture  examinations. 

A  single  rolled  sheet  of  Ti  8  w°0  Mu  was  used  for  all  of  the  samples 
so  that  the  results  of  different  samples  could  be  confidently  compared 
without  worrying  about  variations  in  chemistry  and  initial 
microstructure.  Sheet  samples  were  list'd  to  emphasize  the  shear  nature 
of  the  fracture  since  shear  separation  is  dominant  in  plane  stress.  The 
sample  configuration  used  is  shown  in  Figure  1.  The  surfaces  of  the 
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72 


samples  were  wet  ground  down  to  grit  size  600  to  reduce  the  stress 
concentration  caused  by  surface  scratches. 

Two  heat  treatments  were  used  for  the  samples;  solutionizing  for 
one  hour  at  800°C,  water  quench,  three  day  aging  at  600°C,  air  cool 
(wq-a),  and  solutionizing  one  hour  at  800°C,  20  minute  furnace  cool  to 
600°C,  three  day  aging  at  600°C,  air  cool  (fc-a).  Both  of  the  heat 
treatments  result  in  acicular  +  grain  boundary  o  in  a  5  matrix,  but  the 
furnace  cooled  samples  have  considerably  larger  a  platelets,  as  shown  in 
Figures  2  and  3.  Although  the  strength  and  elongation  is  comparable  for 
both  treatments,  the  furnace  cooled  samples  were  used  for  the 
low- to-moderate  strain  rate  tests  because  the  coarser  platelets 
facilitate  the  determination  of  local  strain  values,  as  will  be  seen 
later  on . 

Tensile,  tests  were  performed  on  two  different  machines,  low  rate 
(8.75x10  ‘s  1)  tests  on  an  Instron  twin  screw  frame  and  the 
low-to-medium  rate  (6.26x10  3s  1  to  8.54x10  “s  1)  tests  on  an  MTS 
servo-hydraulic  system  in  stroke  control.  Even  at  the  moderate  rates 
used  in  these  tests,  data  acquisition  was  a  considerable  problem  due  to 
the  slow  chart  pen  response.  To  record  the  load-elongation  data,  a 
personal  computer  with  an  a/d  board  with  a  data  storage  rate  of  10,000 
hz  was  used  in  con  inaction  with  a  storage  oscilloscope  with  a  writing 
speed  of  10  ysec  per  dot.  The  combination  of  the  oscilloscope  traces 
and  the  digital  data  files  proved  an  effective  means  of  recording  the 
tests,  the  results  of  which  are  summarized  in  Figure  4.‘ 

Figure  4  shows  the  slightly  positive  strain  rate  dependence  of  the 
material.  This  rate  dependence  is  understandably  small,  since  the 
dramatic  increases  in  strength  demonstrated  at  high  rates  generally  do 
not  appear  until  strain  rates  above  approximately  5xl03s  1  . 

The  load-elongation  curves  recorded  in  the  low-to-medium  strain 
rate  tests  showed  a  marked  (Increase  in  the  strain  hardening  rate  with 


r 


i 


increasing  strain  rate.  At  the  highest  rates  tested  the  true 
stress  true  strain  curves  actually  had  a  negative  hardening  slope 
indicating  that  work  softening  occurs.  This  result  suggested  that  (1) 
the  mechanisms  which  caused  hardening  at  lower  rates  are  not  as 


in  decreasing  the  effective  dislocation  link  length,  and 
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there  foro  the  instantaneous  yield  strength  of  the  material  and  (21 
thermal  softening  of  the  whole  gage  length  might  play  a  significant  role 
in  lowering  tin'  yield  strength  at  increasing  rates.  Dislocation  drag 
mechanisms  were  not  considered  because  these  effects  are  generally  not 
detected  at  rates  below  approximately  l.\103s  *. 

Temperature  measurements  o  1  Lb"  gage  length  were  desirable  to  test 
the  second  of  the  these  possibilities.  Thin  film  res i stance  -  typo 
surface  temperature  gages  were  applied  to  a  selected  group  of  samples 
which  were  strained  to  fracture  at  > .  7x  1  ()°s  1  .  A  typical  temperature 
vs.  t  i me  curve  is  shown  in  Figure  a  showing  a  maximum  surface 
temperature  rise  of  about  '1K°(1.  (loin;  ider  i  ng  that  the  test  took  only  'ifi 
msec,  to  Ir. in turn,  it  is  clear  that  there  is  considerable  thermal  lag  in 
the  measin  in..;  system.  This  test,  does,  however,  confirm  the  adiabatic 
nature  of  the  overall  temperature  rise  in  the  gage,  in  complete 

agreement  with  the  strain  rate  window  for  "cpias i si  at  i  c  ad i abat ic i t y , "  as 
de  t  e  rm i ned  by  u 1  sen  e  t  a  I  .  |  I  |  and  bn  1 ve  r  | 2  |  . 

Sine'  the  >.  ib.iilat  ions  and  expel  iment.al  results  indicate  that  a 
fully  adiabatic  model  can  be  effectively  used  to  describe  the 

temperature  in<  reuse  in  the  overall  gage,  previously  measured  yield 
strength  vs.  temperature  information  [  'i  |  was  used  to  determine  an 
upp :  a  >:■:  in  it  e  differential  strain  softening  rate  due  to  thermal  so  I  t  e  1 1  i  ng  . 
When  the  theoret  ii.al  .o  f  ten  i  ng,  rate  was  compared  with  the  measured 
sc. |  t  el.  ill,;  I, lies  it  ••  i  ev.it  ed  st  IM  III  l.lle,  It  was  loUIld  tll.lt  theilll.il 

.  ■ . ;  t  '  ■  n  i  i 1  ;  •  •  •>!  I.)  .  ,n  ly  .l'.c.onnl  f  or  the  dec  reuse  in  worf  ii  mien  iug  i  it  e  up 
t  ■  I  .  :  :  111;  rife,  .  I  i  o;  1  ()  1  ,  ‘  ,  above  wll  i  «.ll  •  t  lie  *  mechanisms  mu  .  t  U  Iso 

,  :  1  ;  : . •  •  :  e  .  f  the  •  i  ! .  •  " : .  t  • .  i  !  so  .  1 ; ow  an  interesting 
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By  reconstructing  the  samples  after  fracture  and  measuring  them  on 
a  profile  projector,  total  elongation  measurements  were  made.  These 
were  complimented  by  pre-  and  post  - f racture  measurements  of  gage  marks 
shown  in  Figure  1.  By  measuring  the  gage  marks  outside  of  the  slightly 
necked  regions  of  the  gage  lengths  uniform  elongation  measurements  were 
also  taken.  Elongation  data  demonstrate  a  decrease  in  elongation  with 
increasing  strain  rate. 

(ii)  Microscopic  Examinations 

Analysis  of  the  fracture  surfaces  in  the  scanning  electron 
microscope  consisted  of  two  steps  (.1)  complete  characterization  of  the 
various  features  found  on  the  surfaces  using  stereo  microscopy 
techniques  and  (2)  nickel  plating  of  the  surfaces  followed  by 
sectioning,  grinding,  polishing  and  etching  to  reveal  the 
microstructures  beneath  the  various  surface  features.  The  combination 
of  surface  and  subsurface  micrographs  provides  significant  insight  into 
the  final  fracture  sequence. 

A  low  magnification  overview  of  the  fracture  surfaces  shows  the 
typical  features  one  might  expect;  central  crack  areas  covered  with  mode 
I  and  III  ductile  dimples  and  lamella  resulting  from  separation  along 
grain  boundaries,  and  shear  lip  regions  covered  with  regular  mode  I,  II 
and  III  ductile  dimples  lying  in  directions  of  high  shear  relative  to 
the  tensile  axis.  Although  the  relative  amounts  of  central  and  shear 
lip  areas  showed  no  dependence  on  strain  rate,  the  water -quenched  and 
aged  samples  display  consistently  more  central  crack  area  (on  the  order 
of  50°ol  as  compared  to  the  furnace  cooled  and  aged  (on  the  order  of 
10Y)  . 

At  higher  magnifications  other  unusual  and  unexpected  features 
become  apparent.  The  most  dramatic  of  these  features  are  spheroidized 
surface  debris,  which  are  found  primarily  along  the  dimple  walls  at  the 
borders  of  the  shear  lip  areas.  An  example  of  this  debris  is  displayed 
in  Figure  t> ,  which  shows  matching  dimples  from  both  of  the  fracture 
sui"!  ii.es.  As  (.an  be  seen  in  Figures  oc  and  od,  sphere  id  i  zed  debris  is 


found  scattered  across  the  neighboring  suit  ace. 


The  surface  areas  where  the  spheroidized  debris  is  found  show 
another  interesting  feature  at  high  magnifications,  a  fine  roughening 
(."microroughening")  .  This  "microroughening , "  shown  in  Figure  7,  is 
always  found  in  the  areas  containing  the  spheroidized  debris.  Because 
it  has  a  grey  tint  in  the  SF.M  at  low  magnif  ications ,  "microroughening” 
facilitates  the  location  of  the  regions  containing  the  unusual  features 
which,  for  reasons  which  will  later  become  apparent,  are  called  Inner 
Surface  Shear  Zones  (.ISSZs). 

Dimple-free  bands  of  material  are  found  adjacent  to  some  but  not 
all  of  the  ISSZs.  These  areas,  an  example  of  which  is  shown  along  the 
right  side  of  Figure  8,  are  unusual  not  only  because  of  their  almost 
complete  lack  of  dimples,  but  also  because  they  contain  striations  which 
are  oriented  generally  in  the  directions  parallel  to  the  shear 

directions,  as  indicated  by  the  shear  dimples  around  the  spheroidized 
debris  covered  areas.  When  matching  fracture  surfaces  are  compared,  the 
dimple-free  bands,  which  will  be  called  Open  Surface  Shear  Zones 
(OSSZs),  almost  completely  surround  the  shear  lip  areas,  with 

approximately  one  half  of  the  OSSZs  on  each  of  the  matching  fracture 

surfaces . 

To  get  a  better  idea  of  the  relative  positions  of  the  various 
features  just  described,  overall  maps  of  the  surface  structures  were 

constructed,  an  example  of  which  is  given  in  Figure  9  This  diagram 
also  demonstrates  that  the  matching  surfaces  of  the  samples  can  be 
divided  into  a  small  number  of  major  areas,  each  of  which  display  a 

predominant  shear  direction,  as  indicated  by  the  arrows.  The  fact  that 
these  areas  display  common  shear  directions,  as  determined  by  the 
matching  shear  dimples  and  the  striations  on  the  OSSZs,  is  a  clear 
indication  that  the  final  separation  occurs  as  a  small  number  of  major 
shear  events,  three  in  the  case  of  Figure  9.  It  should  also  be  noted 
that  the  ma jo  -  shear  directions  displayed  on  the  fracture  surfaces  are 
in  every  case  compatible  with  the  overall  motion  of  the  separating 
sample  segments,  i.e  adjacent  major  shear  segments  never  display 

opposing  shear  directions.  In  Figure  9  this  compatibility  is  maintained 
across  the  surface  by  the  interior  region  of  twist  which  requires 

different  shear  directions  on  either  suit'. 

The  next  step  in  Liie  microscopic  investigation  was  to  coat  the 
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characterized  fracture  surfaces  with  a  heavy  layer  of  nickel  to  preserve 
the  fine  surface  structures  and  then  section  the  samples  parallel  to  the 
tensile  axis.  After  grinding,  polishing  and  etching,  the  subsurface 
microstructures  were  examined  and  compared  with  the  corresponding  stereo 
surface  micrographs.  In  the  case  of  the  typical  ductile  dimples,  the 
sections  demonstrated  that  the  crack  path  is  influenced  by  the  acicular 
a  precipitates  and  also  by  grain  boundary  a  with  the  a-(5  interfaces 
serving  as  preferred  sites  for  separation,  as  demonstrated  in  Figure  10. 
This  result  is  in  agreement  with  earlier  studies  of  Ti  Sw°„  Mn  and  other 
a-5  titanium  alloys. 

When  the  material  under  the  unusual  surface  features  (ISSZs  and 
OSSZs)  was  investigated,  an  unexpected  phenomenon  was  observed;  the 
microstructures  clearly  showed  extremely  high  shear  strain  localized 
into  thin  layers,  as  depicted  in  Figure  10a.  These  shear  bands  range  in 
thickness  from  1  to  3  mm  under  the  surface  and  are  occasionally  found  to 
intersect  each  other  under  surface  irregularities,  as  shown  in  Figure 
10b,  leaving  regions  of  rotated  material  between  the  two  bands.  These 
bands  of  concentrated  shear  provide  a  critical  clue  to  the  explanation 
of  the  unusual  features. 

3 .  VERIFICATION  OF  ADIABATIC  SHEAR 

Tlie  presence  of  sphero id i zed  surface  debris  posed  an  intriguing 
question  as  to  the  nature  of  its  formation.  Spheroidizat ion  due  to  the 
reduction  of  surface  area  of  molten  material  and  high  temperature 
oxidation  were  selected  as  the  two  most  probable  causes  of  these 
unexpected  surface  features.  In  both  cases,  very  high  surface 
temperatures  are  required,  but  in  the  case  of  oxidation,  the  temperature 
rise  is  on  the  order  of  600°C,  as  compared  to  approximately  1400°C  for 
melting.  For  this  reason  it  became  important  to  determine  whether  or 
not  the  material  had  to  reach  the  melting  temperature  to  form  the 
sphero id i zed  surface  structures. 

Fo  determine  t  lie  ro]o  of  oxidation,  samples  were  strained  to 
fracture  in  a  chamber  which  was  evacuated  and  backfilled  with  argon  four 
t  irs's  ,  effectively  eliminating  ilmost  all  of  the  available  oxygen  and 
nitrogen.  In  “  ve  ry  case,  the  simples  tested  in  this  inert  atmosphere 
display  the  same  features  as  those1  fractured  in  air,  clearly  indicating 
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that  oxidation  docs  not  play  the  primary  role  in  spheroidization  of  the 
surface  debris . 

It  was  also  important  to  determine  if  the  "microroughening"  found 
on  all  of  the  ISSZs  could  be  considered  a  temperature  effect. 
Unfortunately,  the  situation  of  a  high  temperature  surface  layer  lying 
over  a  much  lower  temperature  substrate  is  very  difficult  to  reproduce 
under  controlled  conditions.  As  a  rough  approximation  of  this 
condition,  areas  of  ductile  dimples  which  did  not  display 
"microroughening"  after  fracture  were  characterized  in  the  SEM,  flash 
heated  to  a  bright  white  color  with  a  neutral  hydrogen -oxygen  flame,  and 
cooled  in  flowing  argon.  With  this  technique  it  was  possible  to 
reproduce  the  fine  surface  roughening  found  in  the  ISSZs,  as  displayed 
in  Figure  11.  The  success  of  the  flash  heating  in  causing 
"microroughening"  strengthened  the  concept  that  the  formation  of  the 
unusual  features  is  closely  related  to  very  high  local  temperatures. 

When  the  conditions  of  concentrated  shear  and  high  temperatures  are 
present.  adiabatic  shear  is  a  logical  explanation.  Although 
concentrated  shear  was  clearly  demonstrated  in  the  subsurface  sections, 
the  presence  of  high  temperatures  had  not  been  d i  rect  1  v  observed  in  any 
of  the  experiments,  and  this  was  a  critical  step  to  confirm  the  presence 
of  adiabatic  shear  bands. 

After  ruling  out  oxidation,  melting  of  the  newly  forming  surface 
areas  was  the  most  likely  cause  of  the  observed  spheroidization,  and  the 
local  surface  temperatures  required  for  this  to  occur  must  reach  at 
least  1400°C.  Temperatures  this  high  not  only  promote  luminescence  by 
chemical  reactions  at  the  fresh  surface  but  also  cause  the  emission  of 
visible,  light  by  black  body  radiation.  It  was  reasoned  that  if  black 
body  radiation  does  occur  during  final  separation  this  radiation  should 
be  visible  to  the  naked  eye  and  also  recordable  by  means  of  fast 
photographic  film. 

Direct  observation  of  light  emission  during  separation  is 
complicated  by  two  factors:  (1 )  the  flash  is  very  short  in  duration  and 
is,  therefore,  easily  missed  and  (2)  separation  is  generally  accompanied 
by  a  sharp  report  which  causes  l  lie  eyes  to  naturally  blink  and  miss  the 
sni  a !  I  ,  short  f  lash.  Even  with  i  lies*'  comp]  icat  ing  factors,  certain 
samples  were  seen  to  >’iTiit  a  i  lent*  light  flash  at  separation  when  the 


F  i  gure  1  1  . 


F  i  nt;  surface  roughening  ( "m  i  croroughen  ing" )  produced  by 
rapidly  heating  initially  smooth  dimples  with  a  neutral 
propane-oxygen  torch  for  .approximately  1.5  seconds. 
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tests  were  conducted  in  a  totally  blacked  out  room. 

To  record  the  observed  light  flash,  ASA  1000  film  (Kodak  Recording 
Film  2475)  was  used  with  a  35mm  camera  and  a  Vivatar  macro  lens.  Prior 
to  testing,  all  of  the  lights  in  the  room,  including  the  indicator 
lights  on  the  tensile  testing  system  and  the  CRT  of  the  computer  were 
blacked  out  or  turned  off  to  make  the  room  as  dark  as  possible.  This 
made  taking  the  photographs  a  simple  matter  of  opening  the  shutter  prior 
to  the  test  and  closing  the  shutter  after  separation.  The  initial 
images  recorded  by  this  method  clearly  show  bright  particles  apparently 
ejected  from  a  localized  region  in  space.  To  more  directly  link  these 
bright  objects  to  the  newly  forming  fracture  surfaces,  a  variety  of 
backlighting  conditions  were  tested  until  one  was  found  which  gave  the 
location  of  the  light  relative  to  the  samples.  Figure  12  is  an  example 
of  the  dim  profile  of  the  sample  along  with  the  bright  particles  ejected 
during  separation. 

In  the  case  of  the  test  shown  in  Figure  12,  three  bright  areas 
could  be  seen  at  the  newly  formed  surfaces,  and  when  the  surfaces  were 
examined  with  the  SEM  it  was  found  that  these  three  regions  contained 
the  on  iv  areas  of  considerable  spheroidization.  This  is  a  very 
important  result,  since  it  clearly  links  the  formation  of  the  unusual 
.-.pile  ro  uli/.ed  dimple  walls  and  surface  debris  with  the  emission  of  light. 
An  extensive  literature  search  into  light  emission  of  separating  metals 
turned  up  a  short  report  of  this  effect  in  the  tensile  fracture  of 
another  titanium  alloy  [4],  In  this  study  it  was  found  that  light 
emission  could  be  enhanced  by  enriching  the  oxygen  content  of  the 
surrounding  atmosphere  but,  more  importantly,  it  also  showed  that  the 
flash  persisted  under  vacuum,  in  complete,  agreement  with  our  inert 
atmosphere  resu Its. 

The  light  emission  results  are  confirmation  of  the  tremendous  local 
temperatures  which  are  generated  at  separation,  and,  when  combined  with 
the  observed  localized  subsurface  shear,  corroborate  the  concept  of 
adiabatic  shear  bands  playing  a  vital  role  in  the  formation  of  the 
unusual  surface  features.  This  constitutes  the  first  experimental 
evidence  of  adiabatic  shear  in  tensile  tests  and  also  discounts  the 
concept.  th  it  considerable  compressive  stresses  are  necessity  tor 


adiabatic  shear  initiation,  since  these  forces  are  obviously  not  present 
in  these  tests. 


4.  APPLICATION  OF  ADIABATIC  SHEAR  TO  THE  FRACTURE  SEQUENCE 

The  relative  positions  of  the  various  fracture  surface  features,  as 
indicated  in  Figure  9,  when  combined  with  the  concept  of  adiabatic  shear 
at  final  separation  can  be  used  to  generate  a  model  of  the  sequence  of 
events  which,  occur  during  fracture.  Calculations  of  the  maximum 
possible  temperature  rise  during  uniform  elongation  clearly  point  out 
that  prior  to  localization,  only  small  Con  the  order  of  60°C) 
temperature  increases  are  possible.  During  localization,  which,  in  our 
case,  occurs  in  the  form  of  primarily  plane  stress  localization,  this 
value  can  be  increased  approximately  by  an  order  of  magnitude.  This 
still  leaves  a  tremendous  increase  unaccounted  for  if  the  melting 
temperature  is  to  be  reached. 

Local  strain  values  on  the  order  of  12  are  found  under  the  unusual 
features  and  strains  this  high  can  be  used  to  compute  temperature 
increases  far  exceeding  the  melting  point  if  constant  material 
properties  and  fully  adiabatic  conditions  are  assumed.  At  these 
elevated  temperatures,  however ,  the  strength  of  the  material  decreases 
tremendously,  thereby  decreasing  the  possible  temperature  rise.  As  an 
approximation  of  the  possible  temperature  increases  which  can  be 
expected  under  fully  adiabatic  conditions,  yield  strength  versus 
temperature  data  for  Ti-Sw°<,  Mn  sheet  was  approximated  by  a  linear  and  a 
bilinear  curve,  and  this  was  then  used  to  approximate  the  realistic 
temperature  rise  limit  with  strain.  The  results  of  these  calculations 
are  shown  in  Figure  13.  These  results,  however,  do  not  include  any 
consideration  of  the  strain  rate  hardening  effects  which  are  doubtless 
important  during  final  separation  and  serve  to  increase  the  maximum 
aeh ivvab le  temperature . 

The  presence  of  spheroidized  features  on  the  fracture  surfaces  of 
samples  strained  at  quasi-static  rates  is  an  indication  that  the  applied 
strain  rate  of  the  tost  system  does  not  necessarily  provide  an  accurate 
me.isur**  of  the  rate  of  s epa ration.  Elastic  energy  stored  in  the  testing 
system  sail  be  released  into  the  final  deforming  volume  at  a  rate 
controlled  only  by  the  elastic  wave  propagation  velocity.  When  this 


91 


Temp.  Rise  (°C) 


Characteristic  (fully  adiabatic) 
Temperature  Rise 


Temperature  Rise  with  Heat 
Conduction  Considered 


ol 


to  * 10  c.o  go  ‘20  iio  iUO  )Qo  2oo  2lo  Iho  2uo  26 o  3oo  J2o  3*<o  JfcO 


Crack  Tip  Veiocity  (m/sec) 


Figure  13.  Characteristic  (fully  adiabatic)  and  predicted  temperature 
increases  with  crack  propagation  at  velocity  v  for  a  10 
micron  process  zone  in  Ti8Mn.  Approximations  used  for 
predicted  curve  according  to  Rice  and  Levi  [A). 
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stored  elastic  energy  is  considered,  the  final  deformation  rate  can  be 
seen  to  be  nearly  independent  of  the  applied  strain  rate,  even  at  the 
highest  possible  applied  rates.  This  clearly  agrees  with  our 
experiments  and  points  out  the  importance  of  the  elastic  compliance  of 
the  loading  system,  be  it  a  load  frame  or  an  engineering  structure,  in 
the  initiation  of  adiabatic  shear  in  this  type  of  tensile  fracture. 

5.  ADIABATIC  SHEAR  AND  THE  SEQUENCE  OF  TENSILE  SEPARATION 

If  the  sequence  of  events  which  occur  during  separation  culminate 
with  adiabatic  shear  in  high  shear  directions,  the  following  model  can 
be  used  to  consolidate  all  of  the  proposed  mechanisms  after  the  final 
instability. 

1)  After  necking  begins,  the  central  region  of  the  sample  begins 
to  develop  voids,  further  increasing  the  effective  stress  on  the  local 
cross  section. 

2)  "Sample  scale"  shear  bands  begin  to  develop,  greatly 
localizing  the  continuing  deformation  into  high  shear  directions  and 
increasing  the  local  strain  rate,  as  depicted  in  Figure  14a. 

3)  As  the  central  voids  grow  and  coalesce,  the  central  crack 
extends  towards  the  outer  surfaces  and  the  effective  cross  section 
continuously  decreases. 

4)  When  the  effective  stress  and  the  resulting  effective  strain 
rate  reach  high  enough  values,  adiabatic  shear  is  initiated.  At  this 
point,  as  shown  in  Figure  14b,  the  strain  rate  and  the  strain  increase 
to  very  high  values,  causing  the  temperatures  to  increase  dramatically. 

5)  Final  separation  occurs  along  the  thermally  and  geometrically 
softened  high  shear  regions,  forming  the  ISSZs  and  the  OSSZs  ,  as  shown 
in  F igure  14d . 

The  fact  that  Open  Surface  and  Inner  Surface  Shear  Zones  are  found 
at  the  interior  and  exterior  edges  of  the  shear  lip  areas  can  be 
explained  by  considering  two  couf  igu; at  ions  during  final  separation.  In 
the  first  situation,  depicted  in  Figure  13,  the  void  sheet  has  grown  to 
the  outer  edges  of  the  shear  lip  and  the  final  shear  separation  occurs 
through  the  outer  edges  of  the  sample  and  also  through  the  neighboring 
voids  which  have  not  yet  grown  to  coalescence.  Figure  1e  shows  a  more 
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Open  Surface  Sheer  Zone 
Heat  Affected  Zones 


'igure  14.  Model  of  formation  of  unusual  surface  features,  a)  Shear 
is  localized  into  zone  lying  in  high  shear  direction,  b) 
Adiabatic  shear  is  initiated  at  the  intersection  of  the 
concentrated  shear  with  an  open  surface,  c)  Propagation 
of  adiabatic  shear  creates  a  narrow  band  of  softened 
material,  d)  Separation  along  the  softened  region  leaves 
Heat  Affected  Zones  in  localized  areas  adjacent  to  Open 
Surface  Shear  Zone. 
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Figure  15.  Final  separation  along  adiabatic  shear  band  located  a:  ■ 

....  r  ion  ■  --..p a -•->*-  v’g  shear  lip  and  an  external 
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guro  lb.  Final  separation  along  adiabatic  shear  band  located  at  the 
intersection  of  a  separating  shear  lip  and  an  internal 
crack  which  has  been  blunted  at  a  grain  boundary. 


unusual  situation,  that  of  a  void  sheet  which  has  not  grown  from  the. 
central  crack  directly  outward  and  has,  instead,  lias  grown  from  the 
interior  of  the  incipient  shear  lip  towards  the  interior  and  exterior 
open  surfaces.  This  situation  is  unusual  if  one  considers  the 
traditional  moving  crack  tip  model  of  separation,  but  it  is  consistent 
with  the  "double  cup  and  cone"  fracture  process  reported  by  earlier 
examinations  of  the  various  separation  modes  in  tension  [6]. 

A  final  mechanism  which  contributes  to  the  high  temperatures  which 
cause  the  spheroidized  surface  structures  found  in  this  study,  is  the 
heating  effect  of  the  numerous  dislocations  needed  to  form  the  dimple 
walls  during  void  coalescence.  Dimple  walls,  which  are  universal  in 
ductile  fracture,  have  been  demonstrated  to  result  from  cooperative  slip 
in  the  final  ligaments  between  growing  voids,  and  this  rupture  of 
micro- 1 igaments  requires  a  large  number  of  dislocations  concentrated 
into  a  very  small  volume.  The  resulting  high  local  strains  and  strain 
rates  provide  the  necessary  requirements  for  considerable  adiabatic 
temperture  increase.  Calculations  [7]  indicate  that  the  temperature 
rise  from  dislocation  motion  during  rupture  of  micro- ligaments  amounts 
to  several  hundred  degrees  Celsius  leading  to  local  melting. 

6 .  CONCLUSIONS 

11  There  is  a  continuous  decrease  in  the.  workhardening  rate  of 
our  TiSMn  sheet  tensile  samples  as  the  strain  rate  is  increased  from 
8. 75x10  4  sec  1  to  8.54x10°  sec  l. 

21  Adiabatic  heating  of  the  gage  length  during  uniform  straining 
causes  a  measurable,  increase  in  the.  gage  temperature,  and  in  the  TiSMn 
samples  tested,  this  increase  has  a  maximum  value  of  approximately  60°C. 
Fully  adiabatic  thermal  softening  from  this  heating  is  not  sufficient  to 
account  for  the  decrease  in  workhardening  rate  above  strain  rates  of 
approximately  6x10  1  sec  1 . 

3)  Stereo  Scanning  electron  micrographs  of  the  fracture  surfaces 
show  local  regions  which  have  a  finely  roughened  ( m i croroughened ) 
surface  text  ure  and  contain  spheroidizod  dimple,  walls  and  surface 
debris.  These  unusual  features  always  form  at  the  boundaries  of  shear 
lip  zones  adjacent  to  dimple  free  hands  (('pen  Surface  Shear  Zones  I . 


4)  Sections  cut  through  the  fracture  surface  show  that  these 
unusual  features  lie  above  localized  bands  of  highly  concentrated  shear. 
Shear  strain  values  greater  than  12  have  been  measured  in  these  bands, 
and  the  average  width  of  the  bands  is  approximately  2  microns. 

5)  The  microroughened  surfaces  and  spheroidized  surfaces  debris 
are  formed  in  tests  performed  in  an  argon  atmosphere,  ruling  out 
oxidation  as  the  primary  mechanism  responsible  for  their  formation. 

6)  Distinct  flashes  of  light  are  emitted  at  separation,  and  the 
sources  of  emission  can  be  localized  to  the  surface  areas  containing  the 
spheroidized  surface  features. 

7)  Initiation  and  propagation  of  adiabatic  shear  bands  occurs 
during  final  separation  in  localized  regions  of  some  of  the  shear  lip 
areas . 

8)  Adiabatic  temperature  increases  during  uniform  extension, 
localization,  crack  propagation,  adiabatic  shear  and  final  separation 
can  combine  to  cause  the  final  separating  material  to  exceed  the  melting 
temperature  in  localized  areas. 
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A  NUMERICAL  MODEL  FOR  CALCULATING  THE  TEMPERATURE 
RISE  AT  CRACK  TIPS  IN  A  TITANIUM  ALLOY 

K.  Jagannadham*  and  H.  G.  F.  Wilsdorf 
ABSTRACT 

Previous  attempts  to  calculate  the.  temperature  rise  at  a  moving 
crack  tip  used  the  continuum  approach.  The  current  model  is  based  on 
the  plastic  deformation  behavior  in  the  process  zone  which  is  a  smaller 
volume  element  within  the  classical  "plastic  zone."  The  essential 
elements  of  the  now  model  are  loading  rate,  plastic  work  rate, 
dislocation  density,  dislocation  velocity,  and  plastic  strain.  The 
effects  of  the  temperature,  rise  at  the  tip  simultaneous  with  crack 
propagation  have  been  determined.  The  model  predicts  correctly  the 
temperatures  which  developed  during  fracture  as  determined  by 
experiment.  In  addition,  new  insights  of  the  time-temperature 
relationship  during  fast  crack  growth  were  obtained. 
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A  NUMERICAL  MODEL  FOR  CALCULATING  THE  TEMFERATURE 
RISE  AT  CRACK  TIPS  IN  A  TITANIUM  ALLOY 


In  first  estimates  of  the  temperature  rise  ahead  of  the  crack  tip 
by  Krafft  and  Irwin  [1,2],  the  plastic  work  term  within  the  plastic  zone 
was  averaged  and  the  uniform  heat  generated  was  used  to  determine  the 
temperature  rise.  A  refinement  in  the  analysis  by  Rice  and  Levy  [3] 
consisted  of  solving  the  heat  equation  with  the  integrated  plastic  work 
term  over  the  entire  plastic  zone.  The  stresses  and  strains  within  the 
plastic  zone  obtained  from  the  small  scale  yielding  solutions  of  the 
Dugdale  model  were  employed  [4],  However,  the  analysis  is  still 
incomplete  in  the  sense  that  the  temperature  at  the  crack  tip  is 
obtained  using  static  crack  approximations  without  including  the  effect 
of  a  simultaneous  temperature  rise.  It  is  our  opinion  that  the  plastic 
work  term  changes  in  the  region  ahead  of  the  tip  in  the  plastic  zone  so 
that  special  emphasis  should  be  given  to  the  changes  that  occur  within 
tlu>  process  zone. 

We  have  reconsidered  the  problem  of  determining  the  temperature 
rise  in  the  region  ahead  of  the  crack  tip.  Our  calculations,  instead  of 
using  the  continuum  approach  make  use  of  the  dislocation  density  and  the 
associated  plastic  work  term.  The  present  formulations  take  into 
account  the  work  hardening  nature  of  the  plastic  zone  through  the 
equations  connecting  the  flow  stress  and  the  dislocation  density.  Thus 
the  process  zone  where  the  dislocation  density  far  exceeds  its  magnitude 
elsewhere  is  clearly  distinguished.  The  strain  rate  dependence  of  the 
frictional  stress,  the  temperature  dependence  of  the  flow  stress,  the 
frictional  stress  and  that  of  the  thermal  conductivity  are  also  included 
in  the  present  analysis  to  fit  the  experimental  condi ti u.s  The 
expressions  for  the  tensile  component  of  the  stress  field  luted 
with  the  Broberg  dynamic  elastic  crack  configuration  [  "  j  are  used  md 
further  resolved  for  the  shear  stress  component  along  the  slip  planes  to 
account  for  the  orientation  by  a  factor  equal  to  0.3.  The  stress  1  ield 
expressions  tor  the  plastic  crick  coni  iguration  can  he  modi!  jed  bv 
shilling  the  crack  tip  stress  singularity  into  the  middle  nl  the  plastic 
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zone.  The  dimensions  of  the  plastic  zone  are  determined  from 


t  =  i  +  aGb/p"  , 

o  i 

where  t  is  the  shear  component  of  the  stress  field,  t  the  frictional 

*  o 

stress  against  movement  of  dislocations  in  the  lattice,  a  a  numerical 
factor  equal  to  0.5,  G  the  shear  modulus,  b  the  Burgers  vector  and  p  the  , 

dislocation  density.  When  the  stress  associated  with  the  crack  is  lower 
than  the  lattice  frictional  stress,  no  plastic  deformation  takes  place. 

The  plastic  work  rate,  f  defined  as  the  work  done  by  the  crack  tip 
stress  field  in  generating  the  plastic  strain,  is  calc  lated  from 

f  =  oe  =  Opb  v 

where  e  is  the  plastic  strain  rate  and  v  the  dislocation  velocity  in  the 
plastic  zone.  The  velocity  of  the  dislocations  in  the  region  of  the 
plastic  zone  also  varies  with  the  distance  from  the  tip  in  the  form, 

v  =  (t  -  T  1  C  G/G 
o  t 

where  C  is  the  transverse  wave  velocity,  G  is  an  adjustable  parameter 
depending  on  the  velocity  of  the  crack  itself.  The  parameter  G  can  be 
varied  to  adjust  the  dislocation  velocity  and  thus  determine  its  effect 
on  the  temperature.  The  experimentally  correlated  expression  for  the 
velocity  of  the  crack  as  a  function  of  the  loading  rate  parameter  p  in 
the  form 


V  =  (p/968) 

1 

is  list'd  where  p  is  expressed  in  megapascals  (meter)  ‘/sec  and  V  in 
meters/sec.  A  finite  difference  formulation  of  the  one  dimensional 
heat  efju.it  ion  in  the  form, 

6T/6t  =  Ik'd  C)  6 2  T.'  6.x 2  +  f/d-C 

is  used  where  k  is  lie  thermal  cot. duct  i  v  i  ty ,  il  the  density  and  C  tin; 
specific  heat .  The  region  ihe.id  of  t  lie  crack  Lip  in  the  plastic  zone  is 
d.v;d<>d  into  several  intervals  and  the  pi  1st  ic  work  rate 


1  02 


f  is  determined 


as  a  function  of  time  .it  e.acli  point.  Further  substitution  of  the 
plastic  work  rate  into  the  finite  difference  scheme  (Korhonen  and 
h'leemola  (oj)  enables  the  temperature  at  each  point  to  be  determined  as 
a  function  of  time. 

Our  analysis  is  applied  to  the  Ti-  10V-2Fe-3Al  alloy  for  which  the 

material  parameters  may  be  prescribed  by  the  following  values:  the 

shear  modulus,  G  =  41.9  GPa ,  Poisson's  ratio  v  =  0.32,  the  thermal 

3 

conductivity  k  =  10.90  w/m°C,  density  d  =  4.b5  Mgm/m  ,  specific  heat  C  - 

_  o 

4.9  x  10  J/K’g°C,  and  the  Burgers  vector  b  =  2.86  x  10  cm.  The 
frictional  stress  which  opposes  dislocation  generation  and  movement  is 
taken  to  be  one  tenth  of  the  yield  stress  at  £  -  0.002.  The  dislocation 
drag  force  used  in  calculation  of  dissipated  plastic  energy,  is  assumed 
to  be  arising  from  phonon  scattering.  This  term  which  is  estimated  to 
be  the  major  contribution  to  the  viscous  forces  acting  is  given  by 

3kT  v 

t  b  =  Bv  =  -  - j 

drag  lOC^b2  (l-v2/c2)2 

which  already  includes  the  temperature  dependence  and  the  relativistic 
effects.  The  temperature  and  strain  rate  effects  on  the  frictional 
stress  opposing  d i s locat ion  movement  in  the  lattice  are  given  by 

t  =  (1300  +  21. 2S  Inc)f(T) 
o 

where  f(T)  =  (1  -  0.4  T/350)  for  0  <  T  <  350°C,  and  f(T)  =  (0.75  - 
0 . 3T/  200 ')  for  350  •  T  •'  12C0°C. 

The  temperature  dependence  of  the  shear  modulus  is  taken  in  the  form 
G ( T )  =  G ( 0 )  (1  -  0.5T/1200)  for  0  <  T  <  1200°C 

whereas  that  of  the  thermal  conductivity  and  specific  heat  is  given  by 

k(T)  =  k  (01  (1  +  0.6T/12001  , 

PfT)  =  C  (01  i  1  +  0.6T/1200)  . 

Grier  >;  row  th  takes  place  initially  elastically  until  the.  frictional 
stress  is  exceeded.  Then  crack  growth  will  start  and  the.  crack  will 
pinr.eed  into  its  own  plastic  /.one.  Tin*  width  of  the  plastic  none 
venerated  the  crank  tip  is  a  function  of  the  velocitv  of  the  crack, 
loading  rate  and  the  temperature  dependent  frictional  stress.  In  order 
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to  ascertain  the  influence  of  the  crack  velocity  and  loading  rate,  we 
have  determined  the  various  parameters:  temperature  at  the  crack  tip, 
dislocation  density  at  the  tip,  plastic  work  dissipated  at  the  tip, 
dislocation  velocity  at  the  tip,  and  the  plastic  strain  at  the  tip. 

These  are  shown  in  Figures  1  to  4  for  two  crack  velocities.  In  this 

comparison  it  becomes  clear  that  the  softening  arising  from  the 
dissipation  of  energy  at  the  tip  has  a  major  influence  on  the  various 
parameters.  At  slower  crack  velocity,  the  crack  tip  moves  into  the 
plastic  zone  at  a  speed  slow  enough  so  that  it  continues  to  increase  the 
temperature  rapidly  .  However,  when  the  crack  velocity  is  higher,  the 
crack  tip  moves  into  a  region  where  the  temperature  is  lower  which  has 
the  effect  to  slow  down  the  temperature  increase  in  relation  to  the 
origin  of  the  Lime  scale.  Further,  faster  crack  growth  leads  to  higher 
strain  rates  and  therefore,  the  frictional  stress  is  higher  also.  As  a 
result,  the  region  of  plastic  zone  formation  is  reduced.  The  effects  on 
various  parameters  are  illustrated  in  Figures  1  to  4.  We  have  chosen 
two  crack  velocities  =  0.4  m/sec  and  V  =  0.8  m/sec  at  the  loading 

rates  ot  p  =  b80  MPa/m  sec  *  and  p  =  890  MPa/iii"  sec  ^  respectively.  The 

dislocation  velocity  ahead  of  the  crack  tip  is  chosen  such  that  8  =  10 

for  V  =  0.4  m  sec  ^  and  8  =  20  for  V  =  0.8  m  sec  ^ . 

c  c 

It  is  seen  from  Figures  1  and  2  that  the  temperature  ahead  of  the 
crack  tip  for  =  0.4  m  sec  1  is  higher  than  for  V  =  0.8  m  sec 

This  leads  to  a  lowering  of  the  frictional  stress  and  the  siiear  modulus 
whereas  the  thermal  conductivity  and  specific  heat  increase.  As  a 
result,  the  dislocation  density  (Figures  1  and  21,  plastic  work 
dissipated  (Figures  1  and  2),  dislocation  velocity  (Figure  3)  and  the 
plastic  strain  (Figure  41  are  higher  at  slower  crack  growth  rates. 

These  results  are  important  in  understanding  the  formation  of 
molten  regions  at  the  dimple  walls  of  fractured  surfaces.  High  crack 
velocities  are  indeed  responsible  for  exceeding  the  melting  point  at 
specific  regions.  These  experimental  results  do  not  follow  from  the 
previous  analytical  calculations  which  predicted  only  the  changes  of  a 
temperature  rise  that  occurs  in  the  plastic  zone.  Those  previous 
calculations  do  not  give  the  idiabatic  temperature  changes  that  take 
place  around  voids  and  the  rupture  of  microligaments.  It  is  clear  from 
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our  calculations  that  the  melting  point  can  be  exceeded  only  at 
ligaments  involved  in  the  coalescence  of  voids  at  the  crack  tip  where 
the  strain  rate  is  exceedingly  high.  If  the  temperature  rise  in  the 
crack  tip  plastic  zone  alone  were  to  explain  the  formation  of  molten 
zones,  the  result  should  be  that  such  molten  zones  are  formed  everywhere 
along  the  edge  of  the  crack.  When  the  crack  tip  moves  very  fast,  the 
size  of  the  "plastic  zone"  is  very  small  and  as  a  result  the  width  of 
the  molten  areas  formed  before  the  final  failure  occurs  should  also  be 
smaller.  This  general  result  is  born  out  by  the  distribution  of  molten 
areas  on  the  fracture  surface.  Rapid  glide  processes  within  the 
ligaments  associated  with  the  crack  tip  contribute  the  additional 
temperature  increase  in  the  order  of  300°C  to  500°C  [7,S].  The 
temperature  increase  calculated  from  the  model  alone  is  well  below  the 
melting  point  since  these  calculations  have  not  included  the  rupture  of 
microl igamonts  between  voids  which  are  reduced  in  thickness  at  this 
stage  to  about  1  pm.  Earlier  work  has  shown  that  the  deformation  to 
final  rupture  follows  strictly  the  rules  of  crystal  plasticity,  i.e., 
shear  is  as  predicted  along  crystallographic  glide  pianos  and  glide 
directions  [9].  Modifying  a  model  for  the  calculation  of  heat 
production  through  dislocation  avalanches  by  Freudenthal  and  Weiner  [10] 
the  following  equation  was  used: 

TvbN  f  VK  \j 

T  *  link-  yXv  j  , 

again  v  is  the  dislocation  velocity,  N  is  the  number  of  dislocations,  X 
is  the  length  of  travel,  and  r.  =  k/d*C.  Estimating  the  shear  stress  for 
this  alioy  at  a  temperature  of  1200°C  to  be  190  MPa  and  assuming  that 

dislocation  avalanches  arc  moving  on  glide,  planes  20  nm  apart  which 

,  12-2  , 
corresponds  to  a  maximum  dislocation  density  of  10  cm  ,  one  finds  a 

temperature,  increase  of  540°C.  This  brings  the  localized  heating  above 

the  me  1 1  j  ng  po  i n t  . 


REFERENCES 


1.  J.  M.  Krafft,  Correlation  of  plane  strain  crack,  toughness  with 

strain  hardening  characteristics  of  a  low,  a  medium,  and  a  high 
strength  steel.  Appl.  Mater.  Res.  3  (1964)  88-101. 

2.  J.  M.  Krafft,  and  G.  R.  Irwin,  Symposium  on  Fracture  Toughness 

Testing.  ASTM,  STP  No.  3S1  (1965)  114-129;  Philadelphia,  PA. 

3.  J.  R.  Rice,  and  N.  Levy,  In  Physics  of  Strength  and  Plasticity. 

Ed.  A.  S.  Argon  (1969)  227-293.  MIT  Press,  Cambridge,  MA. 

4.  D.  S.  Dugdale ,  J,  Hech .  Phys .  Solids  8  (1960)  100. 

5.  K.  B.  Broberg,  Arkiv  for  Fvsik  18  (1960)  159. 

6.  A.  S.  Korhonen  and  H.  J.  Kleemola,  Met.  Trans.  9 A  (1978)  979. 

7.  J.  D.  Bryant,  D.  D.  Makel  and  H.  G.  F.  Wilsdorf,  Mater.  Sci.  Eng. 

7J7  (1986)  85-93. 

8.  J.  D.  Bryant,  D.  D.  Makel  and  H.  G.  F.  Wilsdorf,  In  Metal lurg . 

Appl.  of  Shock-Wave  and  High-Strain-Rate  Phenomena.  Eds.  L.  E. 
Murr,  K.  P.  Staudhammer  and  M.  A.  Meyers,  pp.  723-739.  M.  Dekker, 
New  York  (19S6) . 

9.  H.  G.  F.  Wilsdorf.  Acta  Metal  1.  30  (1982)  1247-1258. 

10.  A.  M.  Frcudenthal  and  J.  H.  Weiner,  J.  Appl.  Phys.  27  (1956)  44. 


108 


Appendix  I 
Appendix  II 
Presentations 
Publications 
Degrees  Awarded 


Materials  Science  and  Engineering,  77  (1986)  85-93 


85 


Observations  on  the  Effect  of  Temperature  Rise  at  Fracture  in  Two  Titanium 
Alloys 

J.  D.  BRYANT,  D.  D.  MAKEL  and  H.  G.  F.  WILSDORF 

Department  of  Materials  Science,  University  of  Virginia.  Charlottesville,  VA  22901  (U.S.A.) 

( Received  February  27,  1985;  in  revised  form  Aprii  26,  1985) 


ABSTRACT 

Two  titanium  alloys ,  the  cv-(3  alloy  Ti-8Mn 
and  the  mctastable  alloy  Ti  -  10V -2Fe  -3AI 
l  where  the  compositions  are  given  in  approxi¬ 
mate  weight  percent),  were  strained  to 
fracture  at  strain  rates  of  up  to  28  s~J  at 
room  temperature.  Investigation  of  their 
fracture  surfaces  by  scanning  electron  micros¬ 
copy  revealed  droplet-like  features  at  the  rim 
of  dimples  and  raised  features  in  the  regions 
between  and  inside  dimples.  These  were 
interpreted  as  molten  remnants  of  micro¬ 
ligaments.  A  calculation  of  the  affected 
volume,  taken  from  stereo-photogrammetry 
data,  confirmed  this  view.  It  was  clear  from 
experimental  evidence  that  the  melting  point 
was  exceeded  in  micro-ligaments  only  at  the 
moment  of  the  final  separation,  and  a  model 
is  proposed  to  explain  the  temperature  rise 
from  contributions  of  three  mechanisms 
which  can  be  distinguished  phenomeno¬ 
logically.  Numerical  calculations  have  shown 
that  the  melting  point  can  be  exceeded 
under  adiabatic  conditions  on  account  of 
the  alloys'  high  yield  stress,  low  heat  con¬ 
ductivity  and  high  stress  concentration 
factor. 


1.  INTRODUCTION 

Metals  which  fail  by  tensile  forces  in  the 
ductile  fracture  mode  initially  undergo 
uniaxial  deformation  until  necking  occurs, 
followed  by  cup-and-cone  fracture  or 
rupture.  At  least  90©  of  the  plastic  defor¬ 
mation  is  transformed  into  heat  which,  in 
standard  specimens  and  at  normal  strain 
rates,  hardly  affects  the  test.  However,  at 
strain  rates  t  higher  than  10  1  s'1  the  test 
may  be  influenced  by  a  temperature  rise 

•  sumo 


[1]  if  the  flow  stress  of  the  material  is  tem¬ 
perature  sensitive. 

For  fracture  tests  at  high  strain  rates  the 
temperature  rise  in  the  specimen  due  to 
uniform  straining  almost  coincides  with 
the  development  of  heat  at  the  crack  tip. 

Local  heating  at  the  crack  tip  has  been 
investigated  theoretically  by  Priest  [2]  and 
by  Rice  and  Levy  [3].  The  relationship 
between  stress  intensity  factor  K  and  strain¬ 
hardening  exponent  n  at  high  crack  propa- 
tion  rates  was  treated  by  Krafft  and  Irwin 
[4].  Reliable  measurements  of  the  tempera¬ 
ture  at  a  moving  crack  tip  in  metals  are 
unavailable  and  this  lack  of  experimental 
data  seems  to  have  prevented  progress  in  the 
understanding  of  this  important  problem. 

The  objective  of  this  investigation  is  to 
elucidate  the  energy  dissipation  at  the  crack 
tip  in  two  titanium  alloys  by  concentrating  on 
the  temperature  rise  in  a  rapidly  moving  crack. 
A  study  of  fracture  surfaces  using  scanning 
electron  microscopy  (SEM)  has  been  carried 
out  and  the  novel  features  seen  have  been 
interpreted.  It  appears  that  for  the  first  time 
the  temperature  rise  in  a  rapidly  moving  crack, 
final  separation  in  alloys  has  been  obtained 
from  experimental  data  within  a  narrow  range. 
A  model  is  being  advanced  to  explain  the 
temperature  rise  at  the  crack  tip  for  ductile 
metals. 


2.  EXPERIMENTAL  PROCEDURE 

Tensile  specimens  were  prepared  from  the 
metastahle  (3  titanium  alloy  Ti-10V-2Fe  3A1 
and  the  a- (3  alloy  Ti-8Mn  (where  the  com¬ 
positions  are  given  in  approximate  weight  per 
cent).  The  Ti-  10V-2Fe-3Al  specimens  were 
fabricated  with  a  19.05  mm  x  4.4  mm  gauge 
section,  with  the  diameter  increasing  to 
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9.53  mm  in  the  threaded  section.  Flat  bar 
specimens  were  made  from  Ti-8Mn,  with  a 
gauge  length  of  8.9  mm,  a  thickness  of 
0.68  mm  and  a  width  of  2.4  mm.  The  heat 
treatments  were  chosen  to  maximize  the 
ductility  of  the  specimens.  To  accomplish 
this,  a  solution  treatment  temperature  of 
780  °C  was  used  to  transform  the  as-received 
Ti-10 V-2Fe-3Al  materials  into  a  two-phase 
alloy  with  roughly  10%  of  the  alloy  left  in 
the  a  phase  in  the  form  of  large  (10-50  pm) 
plates.  The  formation  of  smaller  strengthening 
precipitates  was  accomplished  during  a  500  °C 
aging.  For  the  Ti-8Mn  specimens  a  solution 
temperature  of  800  °C  was  used,  followed  by 
a  furnace  cooling  to  600  °C.  This  alloy 
developed  large  acicular  a  precipitates. 
Specimens  with  identical  thermomechanical 
histories  were  subjected  to  different  strain 
rates,  the  slowest  of  which  was  0.0031  s'1 
up  to  7.8  s'1.  The  tests  were  performed  on 
an  MTS  hydraulic  tester.  As  stress-strain  data 
at  the  higher  strain  rates  could  not  be  recorded 
on  the  standard  x-y  plotter,  a  retention  tube 
oscilloscope  was  used  which  read  the  signals 
of  the  load  and  stroke  outputs  directly.  Cali¬ 
bration  was  provided  through  the  x-y  plotter 
by  simultaneous  recording  during  a  low  strain 
rate  test.  Typical  oscilloscope  traces  are  given 
in  Fig.  1. 

Fracture  surfaces  were  examined  in  a  scan¬ 
ning  electron  microscope  and  stereo-micro- 
graphs  were  taken  of  different  regions.  In 
addition,  the  deformation  underlying  the 


Kill.  1.  A  retention  tube  oscilloscope  was  placed  in 
the  .v  i  s  y  mode  and  Used  to  display  the  load  cell 
output  acainst  the  stroke  output  directly.  Stress 
strain  data  could  thus  be  recorded  up  to  the  blithest 
strain  rate  achievable  on  Ihe  MIS  (cross  head  speed. 
<)  -  m  s  1  l  See  Table  g 


fracture  surface  was  examined  using  a 
profile  technique.  The  upper  portion  of  the 
gauge  length  containing  the  fracture  surface 
was  plated  with  nickel  to  preserve  details  of 
the  surface.  The  sample  was  then  mounted 
in  Bakelite  and  ground  down  to  expose  the 
central  core  of  the  tensile  specimen.  Through 
the  action  of  two  different  etchants,  the 
nickel  plating  could  be  brought  down  to 
the  level  of  the  titanium  alloy,  allowing  an 
unobstructed  view  of  the  microstructure 
and  plastic  region  directly  below  the  fracture 
surface. 


3.  RESULTS 

The  features  of  the  fracture  surface  were 
seen  to  vary  with  the  strain  rate  employed. 
While  both  the  rapidly  strained  and  the 
slowly  strained  specimens  failed  primarily 
by  transgranuiar  fracture,  the  features 
within  the  transgranuiar  regions  were  strik¬ 
ingly  different.  While  the  slowly  strained 
specimens  exhibited  an  undulating  terraced 
surface  (there  was  evidence  of  decohesion 
along  a  phase  precipitates),  at  higher  strain 
rates  the  fracture  path  appeared  unperturbed 
by  the  precipitate  morphology.  In  these 
regions  the  fracture  surface  is  more  uniform, 
in  both  the  size  and  the  shape  of  the  dimples 
on  the  surface  as  well  as  in  the  smoother 
overall  topography.  Figure  2  compares  the 
transgranuiar  fracture  surfaces  of  Ti-lOV- 
2Fe-3Al  specimens  subjected  to  high  and 
low  strain  rates.  Complete  tensile  data  are 
given  in  Table  1. 

Larger  shear  lips  were  found  on  more 
rapidly  deformed  specimens,  as  well  as  an 
increase  in  the  reduction  in  area.  Failure 
through  the  void  sheet  mechanism  was  more 
prevalent  in  the  specimens  strained  at  higher 
rates.  Several  of  the  high  speed  tehts  resulted 
in  fracture  by  the  cup-and-eone  mode  usually 
associated  with  more  ductile  metals. 

Some  anomalous  features  within  the  shear 
lips  of  the  fracture  surface  were  found.  In 
Figs.  3  and  I  the  typical  dimple  structure  and 
surface  texture  from  a  transgranuiar  region 
can  ho  compared  with  regions  within  the 
shear  lip.  The  surface  in  tills  region,  although 
dimpled,  is  rougher  in  appearance.  Higher 
magnification  shows  the  surface  to  he  popu¬ 
lated  with  small  hemispherical  bumps  I  Figs.  5 
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and  6)  in  the  regions  between  and  inside  the 
dimples.  In  addition,  the  perimeter  of  the 
dimples  in  this  region  is  decorated  with  small 
spherical  globules.  The  smoothness  of  these 
features  suggested  that  they  represent  the 
globularizeu  remains  of  micro-ligaments. 
They  were  found  not  only  on  the  specimens 
subjected  to  the  higher  strain  rates  but  also 
occasionally,  on  the  fracture  surfaces  of  the 


Fig.  2.  The  variation  in  fracture  mode  for  Ti -10V- 
2 Fe  3AI  can  be  seen  by  comparing  these  two  micro 
graphs.  in  )  The  slowly  pulled  it  -  U.0053  s'  I  speci¬ 
men  shows  a  perturbed  undulating  fracture  surface, 
compared  with  the  i  b  i  smooth  uniformly  dimpled 
surface  of  the  rapidly  pulled  specimen  1C  -  13  s"1 ). 


samples  strained  at  lower  strain  rates.  In  both 
cases,  these  features  were  always  found  only 
near  the  outer  edge  of  the  specimen,  within 
the  shear  lip. 

A  reasonable  postulate  as  to  the  origin  of 
the  rough  textured  surfaces  within  and 


Fig.  3.  A  higher  magnification  shows  the  deep 
irregular  dimpled  surface  typical  of  a  transgranula r 
fracture  surface  for  Ti-lOV -2Fe-3Al. 


Fig.  4.  The  dimple  structure  in  the  shear  lip  is 
markedly  different.  This  micrograph  shows  the 
transition  of  normal  transgranular  fracture  surface 
topography  to  the  shallow  roughened  surface  near 
the  perimeter  of  the  tensile  Ti  hMn  specimen. 
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Fig.  5.  The  globularized  features  in  the  shear  lip  are 
evidenced  in  this  micrograph  of  the  fracture  surface 
of  Ti  »Mn  deformed  at  a  strain  rate  of  28  s-1. 


Fig.  6.  The  remnants  of  micro-ligaments  in  Ti  10V 
2 Fe  3A1  on  the  perimeter  of  the  dimples  are  seen 
to  be  globularized.  The  area  between  and  inside  the 
dimples  is  populated  with  small  hemispherical 
irregularities. 


between  the  dimples  is  that  the  fracture 
surface  in  these  regions  ha.--  been  coated  with 
tite  molten  remains  of  micro-ligaments.  The 
hemispherical  features  would  result  from  the 
action  of  surface  tension,  as  m  splat  cooling. 
To  test  tile  feasibility  of  tins  interpretation, 
tile  volume  of  material  :n  the  mn  ro-hgament - 
was  estimated  and  compared  with  the  mini¬ 
mum  volume  of  material  ret  paired  ?u  yield 
the  heni'.'piieneai  Airfare  features  (in  the 
;t"Ump!ion  of  a  triangular  ligament  i'rn- 
-ertion  wi’ ii  a  height  anti  a  iu-r  of  1 1  pm. 
the  volume  of  material  witir.n  tins  iig.inn  :, 
i  , f  . '  :s  imagined  to  extend  t ■  •  tile  perimeter 
of  tiie  i-ntire  dimple  i  won  Id  be  -u  1 1  :■  cent  t  ■  ■ 
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Fig.  7.  The  depth  of  the  plastic  zone  can  be  quantified 
by  noting  the  deformation  of  the  a  primary  precipi¬ 
tates  in  the  Ti-lOV  2Fe-3Al  alloy  directly  below 
the  fracture  surface. 


The  profile  technique  allowed  for  quanti¬ 
tative  measures  of  the  plastic  zone  beneath 
the  fracture  surface  by  noting  the  distortion 
of  the  underlying  precipitates  (Fig.  7).  As 
with  the  features  on  the  fracture  surface 
itself,  the  size  of  the  plastic  zone  was  seen 
to  be  a  function  of  the  strain  rate.  Measure¬ 
ments  showed  the  plastic  zone  to  vary  widely 
even  within  a  single  specimen.  The  higher 
strain  rates,  however,  tended  to  produce 
shallower  plastic  zones. 

To  obtain  a  better  view  of  the  underlying 
microstructure  in  the  regions  characterized 
by  the  globularized  features,  a  technique  was 
developed  wherein  the  nickel  plating  was 
selectively  etched  away  from  the  fracture 
surface  such  that  the  direction  of  the  shear 
which  operated  in  that  region  could  be  deter¬ 
mined  I  Fig.  8).  It  was  seen  that  in  these 
regions  the  precipitates  immediately  beneath 
the  fracture  surface  were  heavily  deformed, 
with  shear  strains  as  high  as  5. 


1  DISCISSION 

I  he  experimental  evidence  from  SUM 
stereo-micrographs  suggested  convincingly 
that  m  some  areas  of  the  fracture  .surface  the 
alloy’-,  melting  point  had  been  exceeded  at 
the  moment  of  final  .separation.  However,  it 
o  clear  from  an  analysis  of  fracture  surfaces 
studied  m  tins  investigation  that  the  overall 
fracture  surface  must  have  been  well  below 
tin-  inciting  temperature  of  the  alloy.  <  >nly 
small  patches  from  a  few  micrometers  to 


Fig.  S.  By  overetching  the  protective  nickel  plating 
of  the  fracture  surface,  the  direction  and  magnitude 
to  the  strain  in  the  immediate  region  below  the 
fracture  surface  of  Ti  b\ln  can  be  estimated. 


about  20  gm  in  diameter  showed  droplet¬ 
like  features,  and  these  observations  tvere 
made  on  the  shear  lips  only.  Therefore,  the 
calculation  of  the  temperature  rise  AT  will 
be  made  by  estimating  and  summing  the 
temperatures  produced  by  (i)  the  overall 
plastic  deformation  before  the  crack  develops, 
(ii)  the  plastic  deformation  at  the  crack  tip 
and  ( iii )  the  rupture  of  micro-ligaments. 

It  will  be  assumed  that  the  temperature 
increase  due  to  the  plastic  deformation 
preceding  fracture  is  adiabatic.  This  seems 
to  be  justified  since  (i)  the  alloys  have  a 
limited  ductility,  (ii)  the  time  from  the  start 
of  the  test  to  fracture  is  for  most  specimens 


TABLK 2 

Phv.Mcal  data  for  lira  titanium  alloys 


of  the  order  of  a  fraction  of  a  second  and 
(iii)  the  heat  conductivity  of  titanium  alloys 
is  very  low.  The  adiabatic  temperature  rise 
for  a  homogeneous  solid  and  uniform  defor¬ 
mation  is  usually  given  by 


r1=~  |o 

pc  J 


where  p  is  the  density,  c  is  the  specific  heat, 
a  is  the  stress  and  e  is  the  strain.  If  work 
hardening  is  assumed  to  occur  according  to 
the  usual  power  law,  AT,  can  be  calculated 
as  follows: 


AT,  = 


pc  n  +  1 


with  n  being  the  strain-hardening  exponent. 
When  the  values  listed  in  Table  2  were  used 
and  eqn.  (2)  was  employed,  a  temperature 
increase  AT,  of  21  °C  for  Ti-8Mn  was  found. 
The  value  for  e,  is  of  course  estimated,  but 
it  appears  to  be  reasonable  in  view  of  the 
relatively  small  strain  to  fracture.  Since  AT, 
is  small  compared  with  the  melting  point  of 
the  alloys,  a  variation  in  e,  is  of  little 
consequence  to  the  final  result.  From  an 
equation  that  has  been  employed  in  a  more 
general  approach  to  the  temperature  rise  in 
tensile  specimens  [1],  i.e. 

ae 

AT,  =  (3! 
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Work  hardening  exponent 
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X'alues  for  the  following  alloys 

Ti  10V  2Fe  3AI 

1 160  MPa 

2  pfi  ■  10"8  cm 

4. bCi  Mb  m  3 

4.9  ■  102  J  kg’1  'C“l 

1.09  ■  10  W  m  1  C1 

0.32 

TO  MP.i  m1  2  I  at  20  Cl 

I  10  C.l'a 

II  7  Cjl’.t 
o.  1 


Ti  SUn 

1370  MPa 
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a  value  of  33  °C  can  be  determined;  this 
appears  to  be  high.  Although  unimportant 
for  the  titanium  alloys  investigated  here, 
due  to  uniform  straining  may  affect 
the  temperature  rise  at  crack  tips  for  highly 
temperature-sensitive  metals  and  alloys. 

With  the  initiation  of  a  crack  in  the  center 
of  the  specimen,  the  plastic  deformation 
becomes  localized  in  front  of  the  crack  tip. 
Rice  and  Levy  [3]  have  calculated  this 
temperature  increase.  They  assumed  a 
non-hardening  plastic  solid,  a  small  plane 
strain  plastic  zone  at  the  tip  and  the  absence 
of  thermomechanical  coupling  effects.  This 
model  considers  the  case  of  a  moving  heat 
source  ahead  of  the  crack  tip  (at  this  point 
the  model  is  converted  to  a  moving  coordinate 
system).  For  the  stationary  crack  the  plastic 
work  rate  f(x,  y,  t )  is  determined  as 


fix.  t) 


4(1  -  v)Oq2 

- - x 

irG 


where  u;  is  the  maximum  dimension  of  the 
plastic  zone.  For  the  moving  coordinate 
system  the  plastic  work  rate  is 


(fix,  y> 


2l  1  —  iMo02 
n  G 


|1  (1  -  .x  'cji1,2i 


In  this  equation,  5D(y  I  is  the  Dirac  delta 
function  employed  to  account  for  the  dis¬ 
continuity  in  displacement  at  the  plastic 
zone  and  V  is  the  crack  speed.  For  the 
running  crack  with  a  growing  plastic  zone, 
the  work  rate  would  he  the  sum  of  g(x.  y) 
and  fix.  y.  t).  Thus  thegt.v.  y )  term  will 
dominate  under  the  condition  that  V.  the 
speed  of  the  crack  tip,  greatly  exceeds  the 
plastic  zone  growth  rate  dco/dt.  When  the 
equations  for  the  configuration  of  localized 
heat  generation  are  combined,  the  tempera¬ 
ture  rise  ATir,  0.  r  I  at  the  stationary  crack 


The  equivalent  equation  for  the  temperature 
rise  near  the  head  of  a  moving  crack  tip  is 
more  sensitive  to  the  distribution  of  the 
work  rate  (sites  of  heat  generation)  and  is 
described  by 


AT2 


tt1/2  (1  -v2)K2  SV\112 
2  E(pck)112  VtoJ 


0.886 


(1  -v2)K 
E(pck )1/2 


By  substitution  of 


7T  K2 


(7) 


it  is  found  that 


A  T2  =  1.414 


( 1  ~  v2  )Ko0 
E(pck)1'2 


vm 


(8) 


In  this  model,  it  should  be  noted  that  the 
temperature  increase  is  greatest  at  the  crack 
tip  itself  and  drops  to  a  negligible  level  at 
the  plastic  zone  boundary. 

The  calculations  were  made  with  the  room 
temperature  values  of  the  stress  intensity 
factor  K  and  the  stress  o0  which  do  not  take 
into  account  the  heating  in  front  of  the  crack 
tip  during  the  actual  crack  propagation.  It  is 
known  that,  in  high  strength  alloys,  K  is 
increasing  while  o0  will  be  lower  with  increas¬ 
ing  temperature  15-10],  Since  no  specific 
measurements  were  found  in  the  literature, 
the  balance  between  K  and  o0  with  rising 
temperatures  could  not  be  assessed.  In  order 
to  gain  an  impression  of  the  effect  of  an 
increasing  stress  intensity  factor  on  AT2,  the 
dependence  of  AT2  on  V  was  calculated  for 
the  two  alloys  in  accordance  with  eqn.  (8). 
Three  values  for  K  were  assumed  and  AT, 
was  plotted  against  the  speed  of  the  crack 
in  Fig.  9.  It  is  seen  that  at  V  —  20  m  s'1  the 
temperature  rise  is  about  900  °C  when 
K  --  60  MPa  m1/2  for  Ti-10V-2Fe-3Al. 

At  sui  h  a  steep  temperature  rise  at  a  moder¬ 
ate  crack  speed  it  is  reasonable  to  expect 
that  K  would  approach  140  MPa  m1  2  at  the 
time  of  separation.  The  calculation  shows 
that  with  K  -  100  MPa  m1'2  a  temperature 
of  about  1500  V  would  have  been  reached. 

In  Fig.  10.  AT2  has  been  plotted  against  K 
for  three  crack  speeds.  When  a  reasonable 
value  for  V  of  50  m  s  1  for  IV  10V-2Fe-3Al 
is  taken,  it  is  found  that  AT  equals  the 
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Fig.  9.  Temperature  rise  at  the  crack  tip  dependent 
the  crack  speed  l'  for  three  stress  concentration 
factors  A :  — ,  Ti  10  V  -  2Fe-3AI ; - ,  Ti-8Mn. 


Fit;.  10.  Dependence  of  the  temperature  rise  at  the 
crack  tip  on  the  stress  concentration  factor  K  for 
three  crack  speeds  1':  -  ■  Ti  ■  10V~2Fe-3Al; 

Ti  fMn 


melting  point  with  the  room  temperature 
K  value  of  70  MPa  m1  2.  Despite  the 
uncertainty  of  the  changes  that  will  affect 
K  and  a„  during  crack  propagation,  the 
calculations  show  convincingly  that  the  tem¬ 
perature  rise  at  th1  crack  tip  is  substantial 
and  can  approach  the  melting  temperature. 

The  final  contribution  to  the  temperature 
increase  is  caused  by  deformation,  leading  to 
file  rupture  of  many  micro-ligaments  that 
'■ompletes  the  fracture.  After  the  growth  of 
void*  to  a  maximum  >i/.e  tile  remaining 
mii-m-ligaments  were  reduced  m  thickness 
fo  about  1  M'"  with  lengths  of  a  few  micro- 
ni'  NT'  a>  judged  from  earlier  studies  on  a 
variety  of  metals  and  alloys  [11.  1  2  | .  The 
deformation  of  these  micro-ligaments  can 
be  estimated  to  be  an  order  of  magnitude 
larger  t  ban  during  macro-necking  [  1  d  |.  On 


account  of  the  small  size  of  micro-ligaments, 
it  is  appropriate  to  use  a  dislocation  model 
for  calculating  the  temperature  rise  at  the 
point  of  their  rupture. 

A  number  of  attempts  have  been  made  to 
calculate  the  temperature  that  can  be  obtained 
from  the  energy  release  of  moving  dislocations 
on  glide  planes  [14,15].  These  models  assume 
a  dislocation  to  be  a  moving  source  of  heat. 
Freudenthal  and  Weiner  [14]  define  the 
heat  flux  q  for  n  dislocations  producing  a 
total  slip  distance  6  =  nb  as 

t6i> 


with  r  being  the  shear  stress;  u  is  the  disloca¬ 
tion  velocity  and  mb  is  the  distance  between 
dislocations.  After  addressing  the  heat  con¬ 
duction  from  a  heat  source  of  length  A, 
Freudenthal  and  Weiner  arrive  at  the  follow¬ 
ing  equation: 


AT, 


r  //cAiA 
mKt  V  7 r  / 


1/2 


(10) 


where  Kt  is  the  thermal  conductivity  and  k 
is  the  thermal  diffusivity  (k  =  Kt/cp).  An 
equivalent  equation  for  AT3  was  derived 
also  by  Eshelby  and  Pratt  [15]. 

Actual  calculations  show  that  AT3  rises 
by  a  few  degrees  Celsius  if  one  glide  plane 
with,  say,  500  dislocations  moving  at  a  speed 
of  10°  cm  s'1  is  involved.  However,  Eshelby 
and  Pratt  pointed  out  that  the  temperature 
would  rise  to  higher  levels  if  more  slip  planes 
are  active;  clearly,  AT3  increases  with  the 
amount  of  shear  and  smaller  distances  between 
active  planes,  leaving  r  unchanged.  Recently, 
the  above  model  has  been  extended  by 
Armstrong  et  al.  [16]  who  considered  the 
temperature  rise  when  piled-up  dislocations 
on  one  slip  plane  were  suddenly  released  and, 
in  addition,  the  relai.onship  between  shear 
band  density  and  AT3.  However,  since  the 
titanium  alloys  investigated  here  do  not 
contain  hard  particles  behind  which  large 
pile-ups  could  form,  their  model  appears  to 
he  not  applicable  for  the  deformation  of 
micro-ligaments  at  the  point  of  separation. 

Attempts  to  adopt  this  model  without  or 
with  the  proposed  modifications  were  not 
successful.  The  problem  is  to  accommodate 
a  very  complex  geometrical  shape,  i.c.  the 
rim  of  dimples  down  to  a  given  depth,  to 


a  first  approximation  the  following  equation 
is  obtained : 


discrete  glide  planes,  varying  glide  distances, 
varying  distances  between  glide  planes  and 
varying  stresses  to  the  framework  of  the 
model.  It  appeared  to  be  more  appropriate 
in  the  light  of  the  geometry  of  the  rapidly 
deforming  micro-ligaments  to  apply  a  form 
of  the  Freudenthal-Weiner  model  in  which 
glide  takes  place  simultaneously  on  the  most 
active  parallel  glide  planes  in  a  given  volume. 
In  that  case  a  group  of  dislocations  can  be 
considered  to  move  on  one  plane  between 
“insulating"  planes  (Fig.  11).  Obviously, 
this  model  provides  for  the  upper  limit  of 
A T3  and  a  AT3  calculation  according  to 
eqn.  (10)  would  then  represent  the  lower 
limit.  The  temperature  rise  for  active  planes 
near  the  "insulating”  planes  would  have  an 
intermediate  value.  When  the  heat  conduction 
in  the*  direction  is  neglected,  the  temperature 
distribution  given  by  Carslaw  and  Jaeger 
[  17]  is 
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l  is  the  distance  between  the  active  glide 
plane  and  the  "insulated"  boundary  plane 
and  should  have  a  very  small  value;  then  the 
second  term  of  eqn.  (Ill  can  be  neglected. 
When  q  according  to  eqn.  (9)  is  inserted,  to 


Fur  1  1.  MorM  of’ dislocation  motion  on  glide  plant* 
f  which  lakes  place  in  .1  crystal  volume  with  many 
active  glide  planes.  Planes  I  have  been  called 
“insulating  planes  by  Fn  udenthal  and  Weiner  (1-1) 
to  signify  that  virtually  no  heal  loss  to  the  surround¬ 
ing  volume  takes  place  with  reference  to  plane  (' 

■  IbMwirg  according  to  ret.  1  1  1 


t8 


(12) 


In  calculating  AT3,  for  the  rupture  of 
ligaments,  values  for  r,  p  and  c  have  to  be 
taken  that  are  adjusted  for  operating  at  the 
increased  temperature  at  the  crack  tip 
which  with  ATj  30  °C  and  AT2  1500°C 
is  estimated  to  be  1530°C.  Since  measure¬ 
ments  for  the  three  parameters  above  are 
not  available  for  the  alloys  investigated,  the 
following  estimates  were  based  on  the  trend 
of  similar  titanium  alloys  at  higher  tempera¬ 
tures.  t  is  assumed  to  be  350  MPa,  p  to  equal 
4.3  g  cm-3  and  c  to  equal  0.8  X  107  erg  g'1 
°C~1.  If  it  is  assumed  that  the  maximum 
dislocation  density  is  involved  in  the  rupture 
process  (i.e.  about  1012  dislocations  cm'2) 
and  further  a  local  strain  for  micro-ligaments 
is  calculated  (7m  ^5),  then  6  —  150  nm  and 
l  =  35  nm.  According  to  eqn.  (12)  the  tem¬ 
perature  rise  AT3  equals  308  °C.  This  result 
is  not  claimed  to  be  precise;  however,  there 
can  be  little  doubt  that  AT3  is  in  the  range 
of  several  hundred  degrees  Celsius.  Com¬ 
parable  calculations  using  the  equation  of 
Eshelby  and  Pratt  [15]  for  a  number  of 
closely  spaced  active  glide  planes  gives  results 
for  the  temperature  rise  of  the  same  order 
of  magnitude. 

ThesumA  T,  F  AT2  +  AT3  is  -  .utl840°C 
and  thus  exceeds  the  melting  point  of  the 
alloys  by  about  240  °C.  If  the  extremely  short 
time  for  the  rupture  of  micro-ligaments  is 
considered,  it  is  reasonable  to  say  that,  at  the 
moment  when  the  melting  point  is  exceeded, 
the  voids  contain  a  vacuum  and  a  deposition 
of  molten  or  vaporized  material  can  take 
place  on  the  dimple  surfaces  which  remain 
in  their  original  shape.  This  would  explain 
the  minute  hillocks  and  the  frozen  droplet 
configurations  observed  (Figs.  5  and  0)  while 
in  border  cases,  where  only  the  utmost  tips 
have  reached  the  melting  point,  no  such 
deposits  are  visible  (Fig.  4). 

The  theoretical  treatment  is  based  on  the 
contention  that  the  processes  leading  to  the 
temperature  rise  are  adiabatic.  The  low  heat 
conductivity  of  the  titanium  alloys,  which 
is  40  times  lower  than  that  of  copper,  was 
given  earlier  as  an  important  reason.  Adiabatic 
effects,  however,  should  also  be  discernible 
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in  the  behavior  of  the  plastic  zone.  Figures  7 
and  8  were  evaluated  and  in  places  show  a 
shear  of  about  5  which  is  a  clear  indication 
that,  within  the  given  time  interval,  adiabatic 
conditions  prevailed  during  crack  propagation. 
A  case  for  adiabatic  behavior  in  titanium 
during  dynamic  plastic  deformation  has  been 
made  earlier  by  Culver  [18]. 


5.  CONCLUSIONS 

Experimental  evidence  pertaining  to  the 
temperature  at  a  crack  tip  was  obtained  by 
SEM  stereo-photogrammetry  of  fracture 
surfaces.  The  novel  features  observed  were 
droplet-like  irregularities  at  the  rim  of  dimples 
and  numerous  small  bumps  at  the  bottom  of 
the  corresponding  dimples.  An  estimate  of 
the  affected  volume  of  micro-ligaments 
supported  the  view  that  in  certain  limited 
areas  a  very  small  volume  was  heated  above 
the  melting  point  at  the  moment  of  final 
separation.  The  experimental  conditions, 
i.e.  deformation  to  fracture  at  a  high  strain 
rate,  and  certain  properties  of  the  alloys 
such  as  high  yield  stress,  low  thermal  con¬ 
ductivity  and  high  stress  concentration  factor 
made  a  temperature  rise  above  the  melting 
point  a  distinct  possibility. 

It  was  thought  that  three  mechanisms 
could  be  distinguished  which  were  associated 
with  the  temperature  rise  in  the  specimen: 

(i|  uniform  straining  to  necking;  (ii)  energy 
dissipation  at  the  crack  tip;  (iii)  deformation 
to  fracture  of  micro-ligaments.  Calculations 
of  temperature  increases  for  the  three  mech¬ 
anisms  yielded  approximately  30  °C,  1500°C 
and  300  °C  respectively  and  thus  supported 
the  conclusions  above.  The  proposed  model 
is  a  first  attempt  to  rationalize  a  melting 
phenomenon  at  the  fracture  of  two  titanium 
alloys  and  more  detailed  computer  calculations 
are  in  progress. 


This  research  was  sponsored  by  the  Office 
of  Naval  Research  under  Contract  N00014- 
82-0309,  NR  031-840.  We  wish  to  thank 
Dr.  Bruce  MacDonald  for  his  interest  and 
valuable  discussions. 


REFERENCES 

1  H.  P.  Stuwe,  Z.  Metallkd. ,  61  (1970)  180. 

2  A.  H.  Priest,  Proc.  Int.  Conf.  on  Dynamic 
Fracture  Toughness ,  Welding  Institute,  American 
Society  for  Metals,  Metals  Park,  OH,  1977,  p.  95. 

3  <1.  R.  Rice  and  N.  Levy,  in  A.  S.  Argon  (ed.), 
Physics  of  Strength  and  Plasticity,  Massachusetts 
Institute  of  Technology  Press,  Cambridge,  MA, 
1969,  p.  177. 

4  J.  M.  Krafft  and  G.  R.  Irwin,  Proc.  Symp.  on 
Fracture  Toughness  Testing  and  its  Applications, 
ASTM  Spec.  Tech.  Publ.  381,  1965,  p.  114. 

5  H.  G.  Baron,  J.  Iron  Steel  Inst.,  London,  182 
(1956)  354. 

6  K.  J.  Marsh  and  J.  D.  Campbell,  J.  Mech.  Phys. 
Solids,  11  (1963)  49. 

7  U.  S.  Lindholm,  A.  Nagy,  G.  R.  Johnson  and 
J.  M.  Hoegfeldt,  J.  Eng.  Mater.  Technol.,  102 
(1980)  376. 

8  A.  K.  Sechdev  and  J.  E.  Hunter,  Metall.  Trans. 

A,  13  (1982)  1063. 

9  S.  L.  Semiatin,  G.  D.  Lahoti  and  S.  I.  Oh,  in 

J.  Mescall  and  V.  Weiss  (eds. ),  Material  Behavior 
Under  High  Stress  and  Ultrahigh  Loading  Rates, 
Plenum,  New  York,  1983,  p.  119. 

10  A.  Kobayashi  and  N.  Ohtani,  in  K.  Kawata  and 
J.  Shioiri  (eds.),  High  Velocity  Deformation  of 
Solids,  Springer,  Berlin,  1978,  p.  208. 

11  H.  G.  F.  Wilsdorf,  Mater.  Sci.  Eng.,  59  ( 1983)  1. 

12  H.  G.  F.  Wilsdorf,  Scr.  Metall..  17  ( 1983)  1209. 

13  H.  G.  F.  Wilsdorf,  Z.  Metallkd..  75  (1984)  154. 

14  A.  M.  Freudenthal  and  J.  H.  Weiner.  J.  Appl. 
Phys.,  27  (1956)  44. 

15  J.  D.  Eshelby  and  P.  L.  Pratt,  Acta  Metall.,  4 
( 1956)  560. 

16  R.  W.  Armstrong,  C.  S.  Coffey  and  W.  L.  Elban, 
Acfa  Metall.,  3G  (1982)  2111. 

17  H.  S.  Carslaw  and  Jaeger,  Conduction  of  Heat 
in  Solids.  Oxford  University  Press.  New  York, 
1918,  224. 

18  R.  S.  Culver,  in  R.  VV.  Rohde,  B.  M.  Butcher. 

J.  R.  Holland  and  C  H.  Karnes  (eds.  I.  Metal¬ 
lurgical  Effects  at  High  Strain  Rates.  Plenum. 
New  York.  197  3,  p.  519. 


Materials  Science  and  Engineering,  81  (1986)  273-292 


273 


Low  Energy  Dislocation  Structures  Associated  with  Cracks  in  Ductile  Fracture 

K.  JAGANNADHAM 

Materials  Engineering  Department,  North  Carolina  State  University,  Raleigh,  NC  27650  ( U.S.A .) 

H.  G.  F.  WILSDORF 

Department  of  Materials  Science,  University  of  Virginia,  Charlottesville,  VA  22901  (U.S.A.) 

(Received  January  31, 1986) 


ABSTRACT 

Ductile  fracture  in  metals  and  alloys  is  pre¬ 
ceded  by  work  hardening  and  necking.  Micro- 
structural  processes,  previously  explored  by 
in  situ  transmission  electron  microscopy,  have 
been  re-examined  in  the  context  of  interac¬ 
tions  between  glide  dislocations  and  crack 
tips.  In  this,  the  parameters  which  signifi¬ 
cantly  influence  experimental  findings  have 
been  critically  reviewed  and  theoretical  as¬ 
pects  of  ductile  fracture  have  been  examined 
including  the  principle  of  energy  minimiza¬ 
tion.  Consideration  was  given  to  the  displace¬ 
ment  of  the  stress  singularity  ahead  of  the 
crack  tip  on  account  of  the  plastic  zone  and 
the  effects  of  thermal  instabilities  during 
crack  growth. 

Crack  propagation  occurs  in  ductile  ma¬ 
terials  mostly  by  void  coalescence.  This  re¬ 
quires  the  nucleation  of  microcracks  ahead 
of  the  crack  tip  as  was  observed  experimen- 
>>,"■  1  agreement  with  theoretical  expecta- 
Microcrack  initiation  was  seen  to  take 

;■■■('  at  dislocation  cell  walls  which  represent 

e  ’  d  frequently  formed  low  energy  dis- 
loca;  j n  configurations  in  work-hardened 
metals,  As  a  consequence,  too,  the  formation 
of  a  plastic  crack  and  its  growth  in  a  virgin 
crystal  are  distinctly  different  from  the  nu¬ 
cleation  and  growth  of  a  crack  in  a  crystal 
with  a  well-defined  cell  structure  that  was 
formed  through  prior  deformation.  Disloca¬ 
tion  multiplication  is  considered  to  take  place 
by  the  supercritical  bowing  of  the  longest 
free  link  lengths  present  within  the  cell  walls, 
as  also  observed  experimentally.  The  micro- 
structural  changes  that  take  place  in  the  plas¬ 
tic  zone  ahead  of  the  crack  tip  are  found  to 
be  more  important  in  fracture  than  the  micro¬ 
scopic  processes  immediately  at  the  tip  itself. 


1.  INTRODUCTION 

Ductile  fracture  of  metals  and  alloys  is 
more  complex  than  other  types  of  failure. 

This  is  so  mainly  because  a  variety  of  micro- 
structural  changes  take  place  in  the  material 
long  before  final  separation  on  account  of  the 
extensive  plastic  deformation  that  occurs 
before  fracture.  This  causes  substructural 
changes  associated  with  dislocation  configura¬ 
tions  which  depend  on  several  variables  char¬ 
acteristic  of  the  specific  materials  tested,  and 
similarly  on  the  imposed  conditions  of  defor¬ 
mation  including  strain  rate  and  temperature. 

The  work  hardening  of  a  material  during 
plastic  deformation  consists  of  an  increase  in 
the  stress  required  for  dislocation  multiplica¬ 
tion  and  movement.  It  is  controlled  by  speci¬ 
fic  materials  parameters  including  crystallo¬ 
graphy,  precipitates  and  stacking  fault  energy, 
for  example.  These  parameters  and  the  result¬ 
ing  microstructural  changes  should  be  incor¬ 
porated  into  theoretical  models  of  crack 
initiation  and  propagation,  since  an  under¬ 
standing  of  ductile  fracture  of  a  material 
obtained  on  the  basis  of  the  work -hardening 
curve  alone  necessarily  remains  rather  qualita¬ 
tive.  However,  using  various  empirical  rela¬ 
tions,  it  becomes  possible  to  correlate  the 
changes  in  microstructural  features  with  the 
ductile  behavior  preceding  fracture.  To  a 
considerable  extent,  these  are  derived  from 
experiments  wherein  foils  of  metals  and 
alloys  are  strained  to  fracture  by  in  situ 
techniques  in  the  transmission  electron  micro¬ 
scope  with  simultaneous  observation  of 
changing  dislocation  patterns.  These  have 
become  very  powerful  tools  in  the  drive  to¬ 
wards  our  understanding  of  ductile  fracture 
[1],  In  this  it  has  become  clear  that  a  general 
principle  of  energy  minimization  of  disloca- 
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tion  configurations  is  most  useful,  namely 
experience  has  shown  that  low  energy  dislo¬ 
cation  structures  are  developed  about  cracks, 
so  that  energy  minimization  introduces  a 
natural  organization  principle  and  simplifica¬ 
tion  into  the  complex  phenomena  [2-4],  As  a 
result,  it  becomes  possible  to  model  the  crack 
behavior  in  terms  of  work  hardening  and  low 
energy  dislocation  configurations. 


2.  GENERAL  CONSIDERATIONS 

Fracture  of  a  material  takes  place  when  a 
crack  can  propagate  on  decreasing  energy  of 
the  system.  The  total  energy  of  the  crack 
configuration  consists  of  different  energy 
terms  that  contribute  to  the  system  [4],  Of 
several  stress-strain  singularities  that  can 
exist  in  solids,  cracks  are  among  the  most 
important.  The  elastic  crack  tip  is  associated 
with  an  inverse  square  root  type  of  singularity 
of  the  stress  field  [5].  For  plastic  crack  tips, 
the  nature  of  the  singularity  was  analytically 
derived  and  it  was  shown  that  it  is  of  the 
same  type  but  shifted  into  the  plastic  zone 
developed  in  front  of  the  crack  tip  [6].  The 
spatial  variation  in  the  crack  tip  stress  field  is 
centrally  important  in  determining  the  plastic 
deformation  associated  with  the  tip.  Thus, 
generation  and  movement  of  dislocations  in 
the  region  ahead  of  the  tip  take  place  under 
the  action  of  the  crack  tip  stress  field. 

The  classification  of  ductile  and  brittle 
fracture  is  derived  from  the  nature  of  the 
crack  tip  [7,  8].  If  it  is  energetically  favorable 
to  nucleate  dislocations  in  front  of  the  crack 
tip  and  simultaneously  to  develop  ledge  steps, 
the  fracture  behavior  is  classified  as  ductile. 

A  propagating  crack  which  intermittently 
generates  lattice  dislocations  and  ledge  steps 
on  the  crack  surface  represents  a  moving 
.source  of  dislocations.  Ductile  fracture  is 
conceived  in  the  general  literature  as  follows. 
In  highly  ductile  materials,  extensive  disloca¬ 
tion  generation  and  movement  take  place  at 
the  crack  tip  so  that  the  stress  field  associated 
with  the  crack  is  screened  by  the  plastic  zone. 
In  this,  the  sharp  elastic  crack  tip  is  replaced 
by  a  rounded  and  blunted  tip  surrounded  by 
the  plastic  zone.  While  this  view  is  correct,  it 
is  incomplete  in  not  taking  into  account  the 
microstruetural  changes  that  actually  take 
place  ahead  of  the  crack  in  the  plastic  zone. 


even  though  their  influence  on  the  crack  tip 
is  most  important.  Specifically,  quite  com¬ 
monly  the  dislocations  in  the  plastic  zone 
ahead  of  the  crack  tip  rearrange  so  as  to  form 
cell  walls  (the  classical  low  energy  substruc¬ 
ture)  and,  if  so,  plentiful  new  crack  nuclea- 
tion  sites  are  thereby  provided.  Therefore,  the 
importance  of  the  principle  of  minimization 
of  energy  in  the  generation  of  the  dislocation 
substructure  associated  with  the  crack  tip 
cannot  be  overemphasized. 

The  formation  of  the  dislocation  configura¬ 
tions  ahead  of  the  crack  tip  is  influenced  by 
several  parameters  dependent  on  the  nature 
of  the  metal  or  alloy.  For  example,  in  close- 
packed  f.c.c.  and  h.c.p.  metals  and  alloys,  the 
formation  or  presence  of  stacking  faults  affects 
the  mobility  of  the  dislocations  [9],  and 
thence  the  resulting  dislocation  structures, 
in  accordance  with  the  stacking  fault  energy. 
In  particular,  in  low  stacking  fault  energy 
materials  such  as  stainless  steels,  slip  is  con¬ 
fined  to  narrow  slip  zones  so  that  dislocation 
cell  structures  cannot  readily  form  [1],  In¬ 
stead,  other  low  energy  configurations  such  as 
dipolar  mats  may  be  generated  [10],  As  a 
result,  the  mechanism  of  fracture  becomes 
different  in  low  compared  with  high  stacking 
fault  energy  f.c.c.  metals.  Additionally,  point 
defects  generated  by  moving  dislocations  in 
the  region  close  to  the  crack  tip  and  their  con¬ 
densation  to  stacking  fault  tetrahedra  provide 
another  important  microstruetural  feature. 

Its  major  effect  is  that  these  interfere  with 
dislocation  movement. 

The  above  certainly  does  not  exhaust  the 
list  of  important  defect  effects  in  fracture. 

The  interactions  of  dislocations  with  several 
other  defect  configurations  have  been  docu¬ 
mented  especially  in  close-packed  f.c.c.  and 
h  .c  .p .  structures  [11].  Radiation-induced 
voids  and  void  clusters  interact  directly  with 
a  growing  crack  in  that  the  crack  front  cuts 
through  spherical  voids  and  soft  coherent 
second-phase  particles  [12, 13],  Thus  the 
resistance  of  the  material  to  crack  propaga¬ 
tion  is  significantly  altered  also  by  these. 
However,  hard  non-coherent  second-phase 
particles  distributed  in  a  soft  matrix  act  as 
sources  of  dislocations  which  in  turn  react 
with  the  substructure  developed  in  the  plastic 
zone.  The  stress  field  associated  with  the 
crack  nucleates  voids  along  the  non-coherent 
interface  of  the  hard  second-phase  particles 


because  of  the  incompatibility  of  deformation 
[14].  Void  nucleation  along  dislocation  cell 
walls  [15]  and  other  high  energy  dislocation 
configurations  such  as  two-phase  interfaces 
[14]  are  accompanied  by  a  decrease  in  energy 
of  the  configuration.  The  nucleation  of  voids 
and  their  coalescence  constitute  a  major 
mechanism  of  ductile  fracture.  Other  micro- 
structural  features  such  as  grain  boundaries 
and  twin  boundaries  also  are  sources  of  dis¬ 
locations  [11] .  At  the  same  time.  fhe  exist¬ 
ence  of  mese  features  provides  obstacles  to 
the  movement  of  the  crack  tip  and,  at  times, 
significantly  alters  its  path.  Other  important 
lattice  and  crystal  defects  affecting  crack 
propagation  are  illustrated  in  Fig.  1. 

The  generation  and  movement  of  dislo¬ 
cations  in  metals  and  alloys  cause  plastic 
deformation,  and  the  density  of  dislocations 
increases  with  increasing  plastic  strain.  As  a 
result  the  free  link  length  of  dislocations  that 
can  act  as  sources  for  further  dislocations 
decreases.  Eventually,  low  energy  configura¬ 
tions  such  as  dislocation  cell  structures  are 
formed  provided  that  the  mobility  of  disloca¬ 
tions  is  adequate  and  at  least  two  different 
Burgers  vectors  are  available.  In  the  absence 
of  a  crack  tip  and  if  a  uniform  stress  is  applied, 
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F iu .  1.  Crack  initiation  and  propagation  in  ductile 
fracture  are  dependent  on  crystal  defects.  The  most 
important  of  these  defects  are  depicted  in  this  sche¬ 
matic  drawing  of  the  microstructural  regime  of 
ductile  fracture  [  1  j. 


the  stress-strain  relation  associated  with  the 
plastic  deformation  of  metals  and  alloys  has 
been  derived  from  the  various  stages  of  the 
development  of  dislocation  substructures 
[16].  The  presence  of  a  crack  (and  the  stress 
field  associated  with  it)  introduces  additional 
dislocations  into  the  region  of  the  plastic 
zone.  Thus,  if  pappj  is  the  density  of  disloca¬ 
tions  due  to  a  uniform  applied  stress  and 
P crack  the  density  of  dislocations  due  to  the 
crack  tip  stress  field,  assuming  a  simple  super¬ 
position,  the  total  density  of  dislocations  is 
given  by  ptotal  =  pappl  +  pcrack.  On  the  basis 
of  much  experimental  evidence  and  the  theory 
of  work  hardening  the  flow  stress  t  is  related 
to  the  total  dislocation  density  through  a 
simple  relation 

r  =  T0  +  aG1bp^  (1) 

where  r0  is  the  frictional  stress  against  disloca¬ 
tion  movement,  G1  the  shear  modulus,  b  the 
Burgers  vector  of  the  lattice  dislocations  in 
the  material  and  a  a  numerical  factor  nearly 
equal  to  0.5.  An  important  application  of  the 
principle  of  minimization  of  energy  is  that  the 
dislocations  rearrange  in  the  form  of  a  cell 
structure.  The  minimization  of  strain  energy 
associated  with  the  cell  structure  compared 
with  a  random  arrangement  of  dislocations 
[3]  accounts  for  a  similar  simple  relation  for 
the  cell  size  L  which  has  been  known  empiri¬ 
cally  for  a  long  time  and  is  given  by 


T~T0 

where  K  is  a  numerical  factor  equal  to  about 
10.  Thus  the  formation  of  dislocation  cell 
structures  is  considered  as  inevitable  provided 
that  the  conditions  mentioned  earlier  exist. 
Other  low  energy  dislocation  arrangements 
are  also  considered  possible  if  the  cell  struc¬ 
ture  cannot  be  formed  under  restricted  mo¬ 
bility  of  dislocations.  The  dislocations  gener¬ 
ated  ahead  of  the  crack  tip  follow  a  different 
pattern  from  those  in  uniaxial  deformation 
because  of  the  specific  nature  of  spatial  stress 
field  variations  associated  with  the  crack  tip. 
The  interaction  of  dislocations  generated  be¬ 
cause  of  the  crack  tip  stress  field  with  those 
present  as  a  result  of  the  uniform  applied 
stress  also  typically  causes  dislocation  cells 
and  in  any  event  forms  an  important  step  in 
the  growth  behavior  of  the  crack  [17], 


In  this  paper,  experimental  findings  are 
presented  on  in  situ  fracture  carried  out  in 
the  transmission  electron  microscope.  The 
behavior  of  cracks  in  ductile  metals  and  alloys 
in  general,  and  the  dislocation  substructure 
developed  in  front  of  the  crack  tip  in  particu¬ 
lar,  will  be  given  special  attention.  Further, 
the  influence  of  several  microstructural  fea¬ 
tures  on  the  ductile  fracture  of  the  material 
is  described.  The  experiments  on  in  situ  frac¬ 
ture  carried  out  in  the  electron  microscope 
r-nn  provide  large  amounts  of  information 
depending  on  several  important  factors.  For 
example,  cracks  introduced  in  well-annealed 
foils  when  subjected  to  a  small  deformation 
at  very  low  strain  rates  can  illustrate  the 
mechanism  of  generation  and  movement 
of  dislocations  on  the  operating  slip  planes 
under  these  rarified  conditions.  Also  the 
formation  of  ledge  steps  that  are  responsible 
for  the  opening  of  the  crack  surfaces  and 
rounding  of  the  crack  tip  can  be  illustrated 
[4],  However,  if  in  situ  experiments  are  car¬ 
ried  out  on  foils  with  no  perforation  and 
initially  not  transparent  or  semitransparent 
to  an  electron  beam  then,  as  deformation 
progresses  and  the  foil  becomes  thinner,  the 
nucleated  cracks  and  voids  grow  in  the  pres¬ 
ence  of  a  well-defined  cell  structure.  As  a 
result  the  foil  becomes  transparent  to  the 
electron  beam  because  of  thinning,  enabling 
direct  characterization  of  the  dislocation  sub¬ 
structure  and  its  influence  on  the  growth  of 
the  cracks  formed  during  deformation  (com¬ 
pare  Section  5).  The  latter  type  of  experiment 
is  relatively  more  important  in  actually  under¬ 
standing  ductile  fracture  in  macroscopic  speci¬ 
mens.  Similarly,  other  in  situ  experiments  are 
useful  in  understanding  the  influence  of  stack¬ 
ing  fault  energy,  the  number  of  available  slip 
systems  at  the  crack  tip  and  the  lattice  fric¬ 
tional  stress  associated  with  the  material. 

Thus,  it  becomes  possible  to  understand  in 
detail  the  different  mechanisms  of  fracture 
as  influenced  by  different  parameters  in  dif¬ 
ferent  materials.  At  the  same  time,  certain 
common  features  of  crack  growth  behavior 
in  ductile  fracture  can  he  obtained  through 
these  studies. 

Thus,  as  outlined,  in  situ  fracture  carried 
out  in  the  electron  microscope  with  a  proper 
choice  of  the  experimental  parameters  is  im¬ 
portant  to  die  understanding  of  ductile  frac¬ 
ture.  The  thickness  of  the  foil,  the  strain  rate 


and  temperature  at  which  the  foils  are  de¬ 
formed  and  the  operating  voltage  of  the  trans¬ 
mission  electron  microscope  are  some  of  the 
conditions  that  can  be  controlled.  The  micro- 
structural  mechanisms  of  fracture  can  be  cor¬ 
related  with  the  material  parameters.  Direct 
observation  and  simultaneous  dynamic  re¬ 
cording  of  the  substructural  changes  enable 
this  type  of  experiment  to  be  very  versatile. 
Modern  analytical  facilities,  namely  scanning 
and  transmission  attachments,  microdiffrac¬ 
tion  capabilities  and  high  voltage  operating 
conditions  for  examining  thicker  foils,  make 
these  experiments  very  advantageous  in  gain¬ 
ing  a  better  understanding  of  the  mechanisms 
of  ductile  fracture  in  macroscopic  specimens. 


3.  EXPERIMENTAL  OBSERVATIONS  OF  CRACK 
GROWTH  IN  ANNEALED  THIN  FOILS 


In  situ  deformation  and  fracture  experi¬ 
ments  can  be  classified  broadly  into  two  cate¬ 
gories,  namely  the  first  category  wherein 
a  perforation  is  present  as  a  result  of  jet 
thinning  of  the  foil  to  produce  an  electron- 
transparent  area,  and  the  second  wherein  a 
non-transparent  or  semitransparent  area  is 
present  without  a  perforation.  In  the  first 
category,  the  foil  area  close  to  the  perforation 
is  fully  transparent  to  the  electron  beam  and 
is  free  from  dislocations  if  the  starting  ma¬ 
terial  was  well  annealed.  The  cracks  present 
at  the  edge  of  the  perforation  move  into  the 
foil  under  the  action  of  the  applied  load. 
Usually,  the  electron-transparent  foils  are 
so  thin  that  plane  stress  conditions  exist  in 
the  foil  on  applying  the  load  [1].  In  addition, 
the  free  surface  exerts  forces  which  affect  the 
mobility  of  the  dislocations.  As  a  result,  no 
more  than  one  or  two  slip  systems  become 
operative  at  the  crack  tip  in  such  thin  foils 
while  the  formation  of  low  energy  dislocation 
substructures  is  hindered  because  of  the  con¬ 
straints  exerted  by  the  free  surface.  The  image 
force  present  in  these  foils  also  changes  the 
orientation  of  dislocations  with  respect  to 
the  free  surface  such  that  on  its  slip  plane 
(whether  or  not  this  is  the  one  usually  ex¬ 
pected)  the  dislocation  is  normal  to  the 
surface.  Macroscopic  fracture  is  far  more 
complex  in  such  foils  in  that  dislocations 
with  several  different  Burgers  vectors  are 
produced  because  of  slip  on  several  systems. 
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In  contrast,  the  simplicity  of  plastic  defor¬ 
mation  ahead  of  the  crack  tip  in  very  thin  foils 
illustrates  some  of  the  microscopic  processes 
at  the  tip,  namely  dislocation  generation  and 
movement  and  simultaneous  formation  of 
ledge  steps  on  the  crack  surface.  Electron  mi¬ 
croscopy  observations  of  a  foil  deformed  in 
situ  in  the  tensile  stage  shown  in  Fig.  2(a) 
illustrates  the  features  mentioned.  The  two 
slip  systems  inclined  to  the  crack  and  leading 
to  ledges  on  the  crack  surfaces  as  revealed  by 
the  slip  traces,  and  the  resulting  opening  of  the 
crack  faces,  are  a  significant  feature  of  this 
micrograph.  The  dislocation  substructure  in 
the  foil,  which  consists  of  dislocation  loops 
and  loop  debris,  is  similar  to  the  characteristic 
deformation  of  the  matrix  in  b.c.c.  metals 
irrespective  of  the  presence  of  the  crack. 


(b) 


Fm.  2.  i  a )  Crack  tip  in  N'b  1  Oat.'VV  produced  by 
transmission  electron  microscopy  in  situ  straining. 

It  should  !>•■  noted  that  two  slip  systems  are  active 
in  iront  of  the  tip.  (b)  Schematic  illustration  for  the 
discrete  dislocation  model  of  a  continuous  plastic 

crack  ,  lattice  dislocations; - ,  crack  disloca- 

!  ions,  i  Reprinted  by  permission  of  Internationaler 
Bin'h Versand  (l.m.b.H.  from  ref.  IS.) 


However,  the  dense  dislocation  substructure 
present  in  line  with  the  direction  of  the  crack 
in  the  foil  is  associated  with  the  crack  tip 
singularity  [19].  A  schematic  illustration  of 
this  type  of  crack  growth  behavior  wherein 
the  crack  tip  generates  dislocations  intermit¬ 
tently  and  moves  forward  in  the  manner 
originally  proposed  by  Orowan  [20]  is  shown 
in  Fig.  2(b). 

The  following  points  have  to  be  taken  into 
account  when  evaluating  results  obtained 
from  in  situ  straining  of  very  thin  foils  in  the 
transmission  electron  microscope  to  fracture. 
The  limited  number  of  slip  systems  by  which 
plastic  deformation  takes  place  prevents  the 
formation  of  low  energy  dislocation  cell  struc¬ 
tures  which  are  found  in  thicker  foils.  Also, 
because  of  the  lack  of  dislocations  which  cut 
each  other  and  thus  do  not  produce  jogs, 
point  defect  production  is  absent.  As  a  result, 
experiments  carried  out  with  very  thin  foils 
are  not  representative  of  macroscopic  ductile 
fracture  behavior  of  a  material.  Ductile  frac¬ 
ture  involves  extensive  plastic  deformation, 
the  generation  and  movement  of  dislocations 
with  several  different  Burgers  vectors,  and  the 
production  and  interaction  of  point  defects. 


4.  EXPERIMENTAL  OBSERVATIONS  OF  CRACK 
NUCLEATION  AND  GROWTH  IN  DEFORMED 
MATERIALS 

In  this  section,  experiments  in  the  second 
category  wherein  a  foil  is  thinned  to  a  thickness 
without  producing  a  perforation  and  thus  more 
nearly  resembles  a  bulk  sample  are  treated. 

A  small  area  of  the  foil  may  be  now  semi¬ 
transparent  to  a  high  voltage  electron  beam 
or  the  whole  foil  may  be  left  intentionally 
non-transparent.  Now  deformation  in  the  in 
situ  stage  gives  rise  to  extensive  deformation 
of  the  whole  foil,  provided  that  suitable  foil 
orientations  are  chosen.  If  so,  plastic  defor¬ 
mation  on  multiple  slip  systems  generating 
dislocations  with  different  Burgers  vectors 
representative  of  deformation  in  macroscopic 
specimens  takes  place  with  simultaneous  foil 
thinning  to  make  portions  of  the  foil  electron 
transparent.  Work  hardening  of  the  material 
due  to  plastic  deformation  and  dislocation 
interaction  is  observed  in  such  experiments 
[1]  In  metals  and  alloys  with  a  high  mobility 
of  dislocations,  rearrangement  into  low  energy 
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dislocation  cell  structures  with  alternating 
orientation  of  the  lattice  within  the  cell  walls 
in  the  opposite  direction  has  been  docu¬ 
mented  [21],  as  shown  in  Fig.  3(a).  A  brief 
description  of  cell  walls  and  their  influence 
on  the  fracture  behavior  of  the  material  is 
in  order  at  this  point. 

In  a  uniformly  strained  material,  the  cell 
size  L  given  by  eqn.  (2)  is  obtained,  which  is 
commensurate  with  the  average  applied  stress. 
However,  in  the  region  near  the  crack  tip,  a 


smaller  cell  size  of  dislocations  is  expected 
because  of  further  subdivisions  of  cells.  The 
cell  walls  are,  in  some  way,  similar  to  grain 
boundaries  wherein  two  lattices  of  different 
orientations  are  joined  together.  When  glide 
dislocations  with  different  Burgers  vectors 
rearrange  and  form  cell  walls  of  width  w,  the 
net  Burgers  vector  of  the  dislocations  in  the 
cell  wall  is  responsible  for  the  misorientation 
and  the  strain  energy  associated  with  it.  The 
larger  the  dislocation  density  and  the  width 


Fig.  3  Direct  observation  of  initiation  of  microcracks  in  beryllium  at  dislocation  cell  walls  during  high  voltage 
electron  microscopy  in  situ  straining  of  a  foil  approximately  4  /jm  thick:  (a)  cell  walls  formed  during  straining: 
1 1)  I  t  he  largest  microcrack  should  be  compared  with  region  A  in  (a). 


of  the  wall,  the  larger  should  be  the  strain 
energy  associated  with  these  cell  walls.  Mis- 
orientations  between  20°  and  30°  and  cell 
sizes  of  0.1-0. 4 /am  have  been  recorded  where¬ 
as  the  dislocation  densities  within  cell  walls 
average  about  5  X  1011  cm'2  [1],  If  disloca¬ 
tions  within  cell  walls  are  to  act  as  sources, 
free  dislocations  of  link  length  of  the  order 
of  10  nm  have  to  be  activated  within  cell 
walls  of  misorientation  greater  than  10°.  The 
necessary  stress  (about  5  X  109  Pa  for  beryl¬ 
lium)  suggests  that  nucleation  of  microcracks 
is  an  alternative  to  further  glide  dislocation 
production  which  becomes  energetically 
favorable,  a  result  substantiated  by  observa¬ 
tions  during  in  situ  fracture  and  documented 
in  Fig.  3(b)  [15,  21], 

A  schematic  illustration  of  a  cell  wall  con¬ 
taining  a  cross-grid  of  screw  dislocations  is 
shown  in  Fig.  4(a)  with  the  nucleation  of  an 
elastic  crack  along  the  interface.  The  strain 
energy  associated  with  the  misorientation 
across  a  boundary  between  the  cell  wall  and 
the  region  free  from  dislocations  is  respon¬ 
sible  for  the  nucleation  of  an  elastic  crack. 

The  stress-free  boundary  conditions  that  exist 
at  the  free  surfaces  of  an  elastic  microcrack 
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Fig.  4.  iaj  Schematic  illustration  showing  microcrack 
initiation  at  a  cell  wall,  (b)  dislocation  model  of 
microcrack  initiation  at  a  cell  wall.  Surface  disloca¬ 
tions  have  been  drawn  in  a  smaller  size. 


can  be  understood  as  the  presence  of  surface 
dislocations  of  the  opposite  sign  called  anti¬ 
crack  dislocations  [4].  The  stress  field  due  to 
the  anticrack  dislocations  cancels  with  the 
stress  associated  with  the  lattice  dislocations 
in  the  cell  walls.  A  schematic  illustration  of 
the  dislocation  configuration  associated  with 
the  cell  wall  and  that  of  the  elastic  crack  with 
anticrack  dislocations,  represented  as  edge 
dislocations  for  convenience,  is  shown  in 
Fig.  4(b).  The  attractive  force  between  the 
dislocations  of  opposite  signs  tends  to  move 
the  lattice  dislocations  towards  the  elastic 
crack  and  thus  to  annihilate  with  the  anti¬ 
crack  dislocations.  As  a  result,  a  ledge  step 
is  produced  on  the  surface  which  in  turn 
removes  the  misorientation  that  is  initially 
introduced  because  of  the  presence  of  the 
dislocations  in  the  cell  walls.  Therefore, 
nucleation  of  a  microcrack  and  further 
annihilation  of  dislocations  within  the  cell 
wall  give  rise  to  formation  of  a  microcrack 
with  free  surfaces  and  reduction  in  strain 
energy  associated  with  the  misorientation. 

The  total  energy  change  may  be  written  in 
the  simple  form 

AE  =  Ee,  ~Eei  +£1;  +  Et  +  EW  (3) 

where  £ei  is  the  strain  energy  associated  with 
the  cell  wall  initially  before  the  nucleation 
of  the  elastic  crack,  Eet  is  the  elastic  strain 
energy  of  the  total  configuration  after  par¬ 
tial  dislocation  annihilation,  Eyj  is  the  free 
surface  energy  associated  with  the  crack  sur¬ 
faces,  Et  is  the  frictional  energy  spent  in  the 
movement  of  lattice  dislocations  and  Ew  is  the 
work  done  by  the  applied  stress.  If  complete 
annihilation  of  the  dislocations  in  the  cell 
walls  takes  place,  Ee  =0.  The  term  £«._  in 
eqn.  (3)  is  the  energy  associated  with  the 
boundaries  of  the  cell  walls.  Further,  the 
above  equation  suggests  that,  if  the  frictional 
stress  is  small,  the  strain  energy  associated 
with  the  cell  walls  can  be  sufficiently  larger 
than  the  surface  energy  before  a  crack  nucle¬ 
ates  that  there  is  a  net  energy  decrease.  There¬ 
fore,  a  critical  misorientation  across  the  cell 
walls  is  needed  above  which  nucleation  of 
microcracks  along  cell  walls  becomes  thermo¬ 
dynamically  favorable. 

At  this  point,  it  is  useful  to  note  that  the 
cell  wall  has  a  boundary  separating  the  dislo¬ 
cation-free  area  from  that  within  which  dis¬ 
locations  are  present.  Across  this  interface, 


there  is  no  material  discontinuity  other  than 
the  orientation  difference.  Therefore  the 
energy  £e  associated  with  the  interface  is 
only  due  to  the  misorientation  and  is  given 
by  the  strain  energy  of  the  dislocations.  When 
the  foil  area  adjacent  to  the  elastic  crack  be¬ 
comes  free  from  any  defects,  the  dislocations 
originating  from  adjacent  cell  walls  and  the 
crack  tip  move  into  the  free  area  of  the  foil 
and  this  will  give  rise  to  further  thinning  and 
relaxation  of  the  stress  field  associated  with 
cell  walls. 


5.  NUCLEATION  AND  PROPAGATION  OF 
CRACKS  UNDER  NORMAL  RATES  OF  DEFOR¬ 
MATION 

In  this  station  the  formation  of  dislocatior 
substructures  in  thicker  foils  in  front  of  the 
crack  tip  will  be  presented  to  illustrate  the 


ture  is  an  important  feature  of  the  micro¬ 
graph.  The  gradual  thinning  of  the  region 
ahead  of  the  crack  tip  made  the  foil  electron 
transparent  to  the  beam.  Simultaneously,  the 
cell  structure  becomes  more  regular,  with 
regions  in  between  the  cell  walls  now  becom¬ 
ing  free  from  dislocations.  A  comparison  of 
the  dislocation  substructure  in  the  thicker 
areas  of  the  foil  adjacent  to  the  crack  flanks 
with  those  observed  in  the  thinned  areas 
illustrates  the  rearrangement  of  dislocations 
into  a  more  regular  arrangement  of  cell  walls 
with  lower  energy.  The  dense  dislocation  sub¬ 
structure  within  cell  walls,  separating  the 
areas  free  from  dislocations,  can  be  recog¬ 
nized  only  within  the  thinned  and  more  trans¬ 
parent  area  ahead  of  the  crack  tip. 

The  gradual  thinning  of  the  foil  ahead  of 
the  tip  requires  the  operation  of  several  slip 
systems,  with  a  minimum  of  two  systems. 

This  will  be  illustrated  for  the  thinning  of 


ductile  fracture  behavior  at  normal  strain 
rates.  Figure  5  shows  the  electron  microscopy 
observations  of  a  region  in  an  Nb-10  at .%  V 
foil  ahead  of  the  crack  tip.  The  operation  of 
several  slip  systems  which  generate  disloca¬ 
tions  with  multiple  Burgers  vectors  and  further 
rearrangement  in  the  form  of  the  cell  struc- 


Fit>.  V  Thinning  in  front  of  crack  tip  in  Nb  10at.%V 
:oil  i  transmission  electron  microscopy  in  situ  strain- 
m t!  i.  Dislocation  cells  are  in  evidence  in  thinned  and 

■nieke,-  '-eg ions. 


gold  foils  during  the  process  of  ductile  frac¬ 
ture.  First,  a  comparison  will  be  made  be¬ 
tween  a  crack  tip  moving  at  a  slow  rate  of 
straining  (less  than  10~4  s"1)  (Fig.  2(a))  and 
a  crack  “tip”  moving  during  strain  rates 
larger  than  10~4  s'1  (Fig.  6).  The  latter  con¬ 
sists  of  a  series  of  microcracks  which  have 
edges  corresponding  to  crystallographic  direc¬ 
tions.  The  spacing  of  the  ensuing  holes  is 
approximately  equal  to  the  dislocation  cell 
diameter  and  nucleation  occurred  at  cell 
walls  as  described  above.  While  the  whole  tip 
configuration  at  first  seems  to  be  similar  to 
the  Dugdale  model  in  Fig.  7,  it  is  clear  (i) 
that  two  glide  systems  are  needed  to  effect 
the  growth  of  each  microcrack  and  (ii)  that 
an  additional  system  of  translation  is  required 
in  order  to  link  the  microcracks  together. 

It  was  found  [22]  that  four  equally  stressed 
glide  systems  of  maximum  resolved  shear 
stress,  namely  ( 1 1 1 ) [  101  ] ,  (lll)[0li], 
(lll)[101]  and  (11 1  )£01 1],  operated  to 
enlarge  each  individual  microcrack,  thereby 
thinning  the  foil  parallel  to  [110]  as  shown 
fora  foil  orientation  of  (110)  in  Fig.  8. 

The  stress  field  associated  with  the  major 
crack  exhibits  a  singularity  ahead  of  the  tip 
in  the  presence  of  a  plastic  zone  [6],  The 
energy  changes  at  the  nucleation  of  a  micro¬ 
crack  in  the  form  of  a  slit  ahead  of  the  major 
crack  have  been  determined  and  are  presented 
in  Table  1.  A  schematic  illustration  of  the 
mechanism  of  nucleation  of  a  microcrack  is 


Fig.  6.  Crack  tip  consisting  of  numerous  microcracks  in  gold  foil  produced  by  in  situ  straining  in  a  transmission 
electron  microscope  at  e  51  10  3s~£  The  edges  are  parallel  to  (110)  and  (112). 
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g.  ",  i  .i  i  Dugdale's  representation  of  a  plastic  tensile  crack.  The  crack  region  is  represented  by  crack  dislocations 
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i n  I  he  void  tip  is  sharp  at  the  opposite  end.  (c)  Void  growth  by  emitting  lattice  dislocations  at  the  tip. 


-h'lwn  in  Fig.  7,  including  the  definition  of 
some  of  the  parameters  presented  in  Table  1. 
In  addition,  the  various  energy  terms  contri¬ 
buting  to  the  total  energy  of  the  configura- 
tion  and  the  constitutive  energy  equations 
i  onneeting  the  various  energy  terms  are  listed 
n  Table  1 .  while  the  details  of  the  analysis 

■  an  be  obtained  from  ref.  23. 

file  total  energy  of  the  crack  void  config- 
iiration  given  by  ,  in  Table  1,  last  column, 
mow-  tlic  decrease  in  energy  associated  with 

■  m  mn  leation  of  the  void  or  microcrack.  The 
.mpor'an!  parameters  that  determine  the 

;•  mu  ieation  of  a  mieroi  rack  are  the 


yield  stress  or  the  frictional  stress  in  the  lat¬ 
tice,  the  initial  distance  between  the  void  and 
the  crack,  and  the  crack  size.  It  is  clear  from 
Table  1  that  a  higher  lattice  frictional  stress 
resulting  in  higher  frictional  energy  is  not 
favorable  for  crack  nucleation.  Similarly,  if 
the  frictional  stress  is  very  low,  plastic  defor¬ 
mation  becomes  very  extensive  so  that  nuclea¬ 
tion  in  the  form  of  a  slit  is  again  not  favored. 
Instead,  nucleation  of  a  void  along  a  cell  wall 
should  be  more  favorable  energetically,  as 
described  earlier  [15]. 

The  development  of  crystallographic  holes 
ahead  of  the  crack  tip  at  high  strain  rates  of 
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deformation  needs  further  consideration.  The 
initial  slit  formation  along  [  110J  by  additional 
slip  on  ( 1 1 1 ) [  10 and  (111 )[01 1]  systems 
gives  rise  to  the  crystallographic  edges.  These 
zigzag  patterns  of  cracks  compare  very  well 
with  those  observed  in  bulk  specimens  de¬ 
formed  in  tensile  machines  at  normal  strain 
rates  [22],  A  schematic  illustration  of  the 
plastic  deformation  by  slip  responsible  for 
the  crystallographic  edges  parallel  to  (112) 
is  presented  in  Fig.  9(a),  clarifying  a  few  very 
important  features  of  deformation.  The  final 
step  in  the  growth  of  the  actual  crack  takes 
place  by  repetitive  operation  of  glide  in  a 
narrow  slip  zone  as  illustrated  schematically 
in  Fig.  9(a)  and  observed  experimentally  as 
shown  in  Fig.  9(b).  High  strain  rates  of  defor¬ 


mation  on  a  limited  number  of  slip  systems 
and  a  high  critical  strain  (greater  than  5  for 
gold)  can  only  be  accommodated  by  the 
movement  of  avalanches  of  substantial  num¬ 
bers  of  dislocations.  The  frictional  energy 
spent  in  the  movement  of  large  numbers  of 
dislocations  is  dissipated  in  the  form  of  heat. 
During  adiabatic  or  near-adiabatic  glide,  a 
temperature  rise  in  the  very  narrow  glide  zone 
will  take  place  which  offsets  work-hardening 
effects.  An  increase  in  the  temperature  of  the 
glide  plane  accompanied  by  softening  is  re¬ 
sponsible  for  continuation  of  slip  on  the 
same  slip  plane. 

An  estimate  of  the  critical  strain  for  work 
softening,  using  power  law  hardening  with  an 
exponent  n  of  about  0.4  [24, 25],  has  shown 


Kig.  9.  Coalescence  of  crystallographic  holes  by  shear  ong  a  very  narrow  zone,  (a)  schematic  representation; 
hi  shear  along  narrow  glide  zones  in  a  gold  foil  (e  ^  lU”®s~M:(c)  zigzag  crack  produced  by  in  situ  straining  in 
in  Nb  l(Ui*  :  V  foil.  Mai  By  courtesy  of  M.  J.  Marcinkowski. ) 


that  it  should  be  greater  than  4.4  with  dr/dT 
—  0.2  MPa  K  _1.  Similar  crack  growth  behavior 
exhibiting  zigzag  crystallographic  edges  has 
also  been  observed  in  b.c.c.  medals  and  alloys. 
Figure  9(c)  sh  'ws  as  an  example  an  Nb- 
10  at M  V  alloy  foil  deformed  in  situ  in  a  ten¬ 
sile  stage.  The  crack  flanks  which  were  thin¬ 
ned  during  the  crack  propagation  have  be.: 
examined  for  the  dislocation  substructure 
retained  in  the  foil.  Extensive  tilting  of  the 
foil  has  shown  the  absence  of  any  dislocation 
substructure  in  the  crack  flanks.  Therefore, 
it  is  concluded  that  the  foil  is  sufficiently 
thin  that  it  cannot  retain  dislocations.  Only 
very  few  segments  of  the  sources  responsible 
for  continuation  of  slip  on  the  crystallogra¬ 
phic  planes  and  those  left  in  the  foil  after 
extensive  thinning  can  be  observed.  These 
are  lined  up  parallel  to  the  crystallographic 
edges  of  the  crack.  In  contrast,  a  dense  dis¬ 
location  substructure  consisting  of  a  well- 
developed  cell  structure  similar  to  that  shown 
in  Fig.  5  has  been  observed  in  the  thicker 
areas  of  the  foil  adjacent  to  the  crack  flanks. 
Thus  it  is  concluded  that,  simultaneously 
with  the  zigzag  crack  propagation,  extensive 
thinning  of  the  foil  due  to  slip  on  several  slip 
systems  takes  place,  even  while  avalanches  of 
dislocations  on  a  few  crystallographic  planes, 
which  result  in  adiabatic  or  thermoplastic  slip 
and  glide  instability,  are  responsible  for  crys¬ 
tallographic  edges  and  zigzag  crack  propaga¬ 
tion  [24]. 

As  we  have  seen  so  far,  there  are  two  insta¬ 
bilities  associated  with  a  crack  tip:  the  first 
results  from  the  initiation  of  microcracks  at 
dislocation  cell  walls,  and  the  second  com¬ 
pletes  local  crack  propagation  by  adiabatic 
shear.  In  addition,  thermal  instability  occurs 
during  the  thinning  of  microligaments.  This 
can  be  deduced  from  the  presence  of  stacking 
fault  tetrahedra  found  in  ligaments  of  gold 
and  copper  [22]  as  shown  in  Fig.  10(a)  for 
gold.  Figure  10(b)  permits  accurate  measure¬ 
ments  to  be  made,  leading  to  the  calculation 
of  a  vacancy  density  (7V  in  gold  of  4  X  10~4 
[26|.  The  vacancies  produced  during  plastic 
deformation  do  not  form  stacking  fault  tetra- 
hedra  spontaneously  but  need  thermal  acti¬ 
vation  to  do  so.  A  minimum  temperature  of 
)o  C  was  estimated  for  the  nueleation  of 
'tacking  fault  tetrahedra  in  gold  [26].  Al¬ 
though  th:>  thermal  instability  is  of  minor 
miportaiM  e  in  cold  and  copper,  it  may  be¬ 


come  significant  in  ductile  materials  with  a 
low  thermal  conductivity. 

6.  CRACK  PROPAGATION  UNDER  LOW  STRAIN 
RATES 

Above,  in  situ  deformation  of  foils  at  strain 
rates  comparable  with  those  observed  in  mac¬ 
roscopic  specimens  and  the  resulting  disloca¬ 
tion  substructures  observed  at  the  crack  tip 
have  been  discussed.  This  section  is  devoted 
to  observations  at  the  other  extreme  of  strain 
rates,  namely  foils  deformed  at  very  slow 
strain  rates  [1].  Under  these  conditions,  glide 
dislocations  do  not  move  in  large  numbers  all 
at  the  same  time  but  rather  singly  so  that  each 
dislocation  comes  to  rest  in  or  near  an  equi¬ 
librium  position.  As  a  result,  thermal  insta¬ 
bility  and  softening  of  the  glide  planes  are  not 
expected.  Instead,  work  hardening  due  to  the 
development  of  dense  dislocation  substruc¬ 
tures  continues  on  numerous  slip  planes.  The 
dissipation  of  frictional  energy  causing  a  rise 
in  temperature  is  here  not  large  enough  to 
reduce  work  hardening.  As  a  result  of  increas¬ 
ing  work  hardening  through  deformation  on 
specific  slip  planes,  the  operation  of  second¬ 
ary  sources  with  different  Burgers  vectors 
begins  so  that  plastic  deformation  can  con¬ 
tinue.  The  final  outcome  is  a  more  uniform 
and  homogeneous  deformation  by  the  genera¬ 
tion  and  movement  of  dislocations  with 
several  different  Burgers  vectors,  as  indeed 
has  been  observed  experimentally  [1]. 

Dislocation  rearrangements  to  lower  the 
energy  of  the  configuration  naturally  gives 
rise  to  low  energy  cell  walls.  In  comparison 
with  the  cell  structure  observed  under  normal 
strain  rates,  the  cell  walls  formed  at  lower 
strain  rates  are  thinner  and  more  regular.  This 
regular  arrangement  of  the  dislocations  in  the 
cell  walls,  betokening  low  energy  configura¬ 
tions,  is  clearly  seen  in  Fig.  11.  In  addition, 
since  a  repetitive  operation  of  glide  on  a  few 
slip  planes  is  here  replaced  by  a  more  uniform 
deformation,  the  crack  tip  becomes  blunted 
and  thus  a  more  rounded  tip  is  observed.  At 
the  same  time,  thinning  of  the  foil  continues 
such  that  the  crack  is  linked  to  the  nearby 
holes  or  voids.  While  nueleation  of  micro¬ 
cracks  is  still  favored  along  the  cell  walls,  their 
growth  cannot  be  confined  any  longer  to 
specific  crystallographic  planes.  Instead,  a 
more  irregular  growth  of  microcracks  due  to 
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I-'ii£.  1  0.  Slacking  fault  tetrahedra  in  a  ligament  remnant  of  gold  (a)  after  in  situ  straining;  (b)  the  size  of  the 
stacking  fault  tetrahedra  is  between  1.5  and  4  nm. 


deformation  along  several  slip  planes  is  ob¬ 
served.  Substantial  rearrangement  of  disloca¬ 
tions  at  low  strain  rates  results  in  thinning  of 
the  foil  to  link  the  grown  microcracks  to¬ 
gether.  These  features  observed  in  a  foil  of 
gold  deformed  at  very  low  strain  rates  are 
illustrated  in  Fig.  12. 


(  HACKS  Nl  (  l.KATKD  AT  INTERSECTIONS 
Of  SI. II’  BANDS 

As  mentioned  earlier,  there  are  several 
parameters  that  influence  the  development 


of  dislocation  substructure  and  the  micro- 
structural  changes  in  ductile  materials.  Of 
these,  the  low  stacking  fault  energy  in  f.c.c. 
and  h.c.p.  metals  and  alloys  is  an  important 
parameter  that  restricts  the  mobility  of  dis¬ 
locations.  As  a  result,  plastic  deformation 
is  by  slip  in  the  planar  mode.  The  limited 
mobility  of  dislocations  bound  by  stacking 
faults  confines  their  movement  along  the 
slip  planes,  instead  of  developing  first  tangles 
and  later  on  cell  structures. 

This  type  of  plastic  deformation  is  charac¬ 
teristic  of  a-brass  and  austenitic  stainless  steel. 
It  is  illustrated  for  type  304  stainless  steel  in 
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Fig.  11.  Straining  at  e  *  10 ~4 s  1  leads  to  extensive  blunting  at  the  crack  tip-  (a)  microcrack  formation;  (b) 
blunted  crack  tip.  The  increase  in  low  energy  dislocation  wails  between  strain  increments  should  be  noted. 


the  electron  micrograph  shown  in  Fig.  13(a) 
[27],  A  schematic  illustration  of  the  con¬ 
tinuous  thinning  of  the  foil  due  to  plastic 
deformation  confined  to  two  intersecting 
narrow  glide  bands  is  shown  in  Fig.  13(b). 
When  thinning  occurs  along  each  band,  a  void 
is  nucleated  at  the  intersection  of  the  bands 
as  a  result  of  the  superposition  of  deforma¬ 
tion  associated  with  both  bands.  Thus,  several 
voids  nucleated  along  the  foils  at  points  of 
intersection  of  slip  hands  grow  as  a  result  of 


further  deformation.  The  growth  of  the  voids 
continues  only  in  the  direction  of  the  thin¬ 
ning  of  the  foil,  as  illustrated  in  Fig.  12(a). 
This  mechanism  of  void  formation  and  growth 
is  similar  but  different  from  that  of  voids 
when  incompatibility  of  plastic  deformation 
and  a  triaxial  state  of  stress  exist  ahead  of  a 
crack  tip  in  a  work-hardened  ductile  material 
[28].  In  the  present  situation  the  inhomogen¬ 
eity  of  plastic  deformation  is  responsible  for 
the  void  nucleation  and  growth. 
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Fig.  1 2.  Development  of  round  holes  at  a  crack  tip  during  the  straining  of  gold  at  e  <  10"4s  1. 


8.  DISCUSSION 

It  is  clear  from  the  results  presented  in  the 
preceding  sections  that  the  ductile  fracture 
behavior  of  metals  and  alloys  critically  de¬ 
pends  on  several  experimental  and  materials 
parameters.  Common  to  all  the  above  descrip¬ 
tions  of  crack  growth  is  the  importance  of 
the  dislocation  configurations  in  the  process 
zone.  In  particular,  low  energy  dislocation 
configurations  of  cell  structures  are  a  striking, 
and  perhaps  initially  unexpected,  microscopic 
observation.  Provided  that  the  dislocations  in 
the  region  ahead  of  the  crack  tip  are  mobile, 
rearrangement  into  ceil  wails  evidently  takes 
place  efficiently.  This  is  accounted  for  by  the 
principle  of  minimization  of  strain  energy 
associated  with  the  dislocation  configurations 
[4 1  and  may  be  ascribed  to  the  energy  changes 
associated  with  the  dislocation  rearrange¬ 
ments  13],  namely 

-  £,  -  Ee  +Et  +  EK  (4) 

where  £'e  is  the  strain  energy  of  the  final  cell 
structure  configuration,  £e  the  strain  energy 
of  the  initial  configuration,  Et  the  frictional 
energy  change  associated  with  the  dislocation 
rearrangement  and  Ew  the  work  done  by  the 
applied  stress.  First  of  all,  it  is  important  to 
note  that  the  term  due  to  the  applied  stress  is 
always  negative  while  the  term  due  to  the 
la  nice  frictional  stress  is  always  positive. 


Therefore,  whether  dislocation  rearrangement 
takes  place  or  not  depends  largely  on  the  fric¬ 
tional  stress.  If  the  frictional  stress  is  too  high, 
such  as  in  the  deformation  of  certain  b.c.c. 
metals  at  low  temperatures,  dislocation  re¬ 
arrangement  cannot  be  expected.  Similarly, 
in  f.c.c.  or  h.c.p.  metals  and  alloys  with  a  very 
low  stacking  fault  energy  the  mobility  of  dis¬ 
locations  is  very  low  and  cell  structures  can¬ 
not  form  at  moderate  dislocation  densities. 

On  the  assumption  that  the  frictional  stress 
is  very  low  and  Et  is  a  negligible  term,  the 
decrease  in  energy  of  the  formation  of  the 
final  configuration  depends  on  Ee  .  The  mini¬ 
mization  of  the  strain  energy  Ee;  in  terms  of 
the  cell  size  [3],  together  with  eqn.  (1),  dir¬ 
ectly  leads  to  the  simple  result  given  in  eqn. 
(2),  with  the  constant  K  of  the  order  of  10. 

These  results,  which  are  valid  for  a  uniform 
applied  stress,  should  now  be  applied  to  the 
region  in  front  of  the  crack  tip.  A  schematic 
illustration  of  the  cell  structure  envisaged  in 
the  region  ahead  of  the  crack  tip  is  shown  in 
Fig.  14.  It  is  important  to  distinguish  the 
dislocation  substructure  formed  in  front  of 
a  crack  tip  in  a  virgin  crystal  from  that 
present  in  a  crystal  wherein  a  crack  is  nucle¬ 
ated  after  considerable  plastic  deformation. 

In  the  former  situation  the  dislocation  con¬ 
figuration  is  formed  in  the  presence  of  both 
the  applied  stress  and  the  crack  tip  stress 
field,  and  the  latter  corresponds  to  the  prior 
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Fig.  13.  (a)  Hole  formation  at  intersection  of  crossing  glide  bands  in  a  type  304  stainless  steel  foil  0.6  gm  thick 
alter  in  situ  straining  in  the  high  voltage  electron  microscope;  (b)  schematic  drawing  depicting  hole  formation  as 
shown  by  (a),  ((a)  By  courtesy  of  Bauer  [27  ].) 


development  of  a  stable  cell  structure  before 
the  c rack  is  nucleated.  Therefore,  growth  of 
the  crack  into  a  region  where  there  is  micro¬ 
crack  nucleation  along  cell  walls  may  be  dis¬ 
tinguished  from  a  region  where  gradual  thin¬ 
ning  of  the  foil  with  a  more  regular  cell  wall 
arrangement  occurs.  Actually,  further  sub¬ 
division  of  a  pre-existing  cell  structure  is  also 
observed  in  the  crack  tip  stress  field.  This  can 
reduce  cell  diameters  to  one-tenth  or  less  of 
those  which  are  a  long  distance  from  the  tip. 

Measurements  of  cell  diameters  in  front  of 
the  crack  tip  in  Figs.  3(a)  and  3(b)  for  bervl- 
lium,  in  Fig.  5  for  Nb  10  at .%  V  and  in  Figs. 
Ilia),  1 1(b)  and  12  for  gold  correspond  to  an 
average  value  of  60  100  nm.  Using  this  value 


of  L  in  eqn.  (2),  the  shear  stress  r  in  the  Nb- 
10at.%  V  foil  is  determined  to  be  1.56  X  106 
Pa  or  E/90  where  E  is  Young’s  modulus.  The 
same  order  of  magnitude  of  shear  stress  was 
also  found  in  other  foils  examined  and  it  is 
considered  to  be  a  reasonable  value.  Here  the 
frictional  stress  r0  was  neglected  in  com¬ 
parison  with  the  resolved  shear  stress  t  in  the 
foil.  It  should  be  noted  that  these  results  are 
applicable  in  foils  with  prior  deformation  so 
that  a  well-developed  cell  structure  pre-exists 
ahead  of  the  crack  tip. 

When  a  crack  is  formed  in  a  well-annealed 
crystal,  the  resulting  dislocation  configuration 
must  correspond  to  the  crack  tip  stress  field 
over  a  wide  area.  In  particular,  we  consider 


Fig,  1  4,  Model  of  a  crack  tip  surrounded  by  disloca¬ 
tion  cells.  (By  courtesy  of  M  J.  Marcinkowski.) 


thu  plane  strain  situation  of  an  elastic  crack  in 
mode  I  present  in  a  thick  foil.  The  shear  stress 
component  responsible  for  the  generation  and 
movement  of  dislocations  is  given  by 


<tv  s  |  r.  ')  )  —  oa 


sin  0  sin  1 0 


(5) 


where  ( r,  >) )  coordinates  are  employed  [4], 
The  shear  stress  component  due  to  the  crack 
is  superimposed  on  the  component  due  to  the 
applied  stress.  Thus,  in  addition  to  the  dislo¬ 
cation  density  pappl  due  to  the  applied  stress, 
the  dislocation  density  pcrack  due  to  the  crack 
tip  stress  field  should  be  determined  and 
superimposed  to  obtain  the  total  dislocation 
density 


(6) 


Implicit  in  eqn.  (6)  is  the  assumption  that  the 
two  dislocation  density  terms  can  be  super¬ 
imposed.  The  resultant  flow  stress  may  be 
given  in  the  form 


and 

R  3 

—  sin  6  sin  ~6  —  r0 
2  r  20 


=  ^Pcrackln{(^-)  }  (8) 

6  L\b0  /crackl 

The  arrangement  of  lattice  dislocations  in  the 
form  of  a  low  energy  configuration  such  as  a 
cell  structure,  a  Taylor  lattice,  a  dipolar  mat, 
a  tilt  wall  or  a  dislocation  boundary  [10]  is 
to  be  expected  in  the  spatially  varying  stress 
field.  Using  eqn.  (2)  in  conjunction  with  eqns. 
(7)  and  (8),  the  cell  size  theoretically  ex¬ 
pected  at  the  crack  tip  can  be  calculated.  It 
is  found  that  the  resolved  shear  stress  com¬ 
ponent  due  to  the  crack  tip  when  calculated 
for  different  applied  stress  levels  [4]  remains 
almost  constant  over  a  limited  length  along 
certain  orientations  of  slip  planes.  Thus  a  cell 
structure  of  uniform  size  along  these  planes 
of  constant  shear  stress  results.  A  one-to-one 
correspondence  between  the  magnitude  of 
shear  stress  and  the  low  energy  cell  structure 
is  observed  [29], 

In  the  fracture  of  ductile  materials,  crack 
nucleation  and  propagation  into  a  region 
containing  stable  cell  structures  are  the  most 
common  case.  The  cell  structure  screens  the 
stress  field  associated  with  the  crack.  In  addi¬ 
tion,  the  crack  tip  in  ductile  materials  is  not 
sharp  because  of  the  formation  of  ledges 
through  multiple  slip  at  the  tip.  As  a  result, 
the  actual  stress  field  does  not  exhibit  a 
strong  singularity.  Therefore,  crack  growth 
necessarily  progresses  via  microstructural 
changes  ahead  of  the  tip,  such  as  foil  thinning 
and  the  microcrack  nucleation  described 
earlier.  The  coalescence  of  microcracks  or 
voids  leads  to  growth  of  the  crack  in  the 
macroscopic  sense. 

When  the  crack  encounters  the  stable  cell 
structure,  the  stress  associated  with  the  tip 
generates  additional  dislocations  from  sources 
present  in  the  cell  walls.  This  is  clearly  illus¬ 
trated  in  Fig.  15  by  the  nucleation  of  disloca- 
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Fig.  15.  Nucleation  of  numerous  dislocations  from  links  in  cell  walls  during  straining  of  an  Nb-10at.%V  foil  in 
the  transmission  electron  microscope,  thereby  continuing  the  thinning  process  at  the  crack  tip.  The  direction 
of  dislocation  motion  is  indicated  by  the  arrow. 


Fig.  16.  Ligament  remnants  in  a  fractured  Nb~10at.%V  foil.  It  should  be  noted  that  cells  are  absent  in  the  areas 
of  the  crack  flank;  incomplete  cells  are  visible  at  the  upper  crack  flank  near  the  non-transparent  foil. 
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tions  in  the  thin  area  of  a  propagating  crack 
tip  in  an  Nb  J0at.%V  foil.  Thus,  in  bulk 
samples,  the  increase  in  dislocation  density 
as  a  crack  advances  into  a  stable  cell  structure 
occurs  by  the  activation  of  sources  from  free 
dislocation  link  lengths  within  the  cell  walls. 
In  contrast,  in  foil  samples,  plastic  deforma¬ 


tion  ahead  of  the  tip  leading  to  thinning  of 
the  foil  such  that  the  foil  is  free  from  any 
dislocation  substructure  as  shown  in  Fig.  16 
is  very  important.  When  the  foil  is  free  from 
dislocations,  it  may  become  brittle  and  final 
separation  can  take  place.  The  same  process 
can  presumably  also  occur  in  the  bulk.  There- 
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fore,  in  the  fracture  of  ductile  materials  the 
nature  of  the  microstructural  features  formed 
ahead  of  the  crack  tip  is  very  important  as 
already  repeatedly  emphasized. 

Restrictions  on  the  mobility  of  disloca¬ 
tions,  such  as  through  the  presence  of  stack¬ 
ing  fault  ribbons  between  partials  in  extended 
dislocations,  have  a  significant  influence  on 
the  development  of  substructural  features. 
This  has  already  been  demonstrated  by  the 
example  of  the  confinement  of  sequences  of 
gliding  dislocations  into  narrow  slip  zones  in 
low  stacking  fault  energy  alloys  and  the  result¬ 
ing  nucleation  of  voids  due  to  foil  thinning 
at  points  of  intersection  of  slip  bands  as  dis¬ 
cussed  in  Section  7.  Conversely,  the  mobility 
of  dislocations  can  be  increased  by  raising  the 
temperature  of  deformation.  Then,  both  con¬ 
servative  and  non-conservative  motion  of  dis¬ 
locations  becomes  relatively  easier  so  that 
polygonization  and  low  energy  dislocation 
cell  structures  can  develop. 

In  this  connection  the  importance  of  strain 
rate  cannot  be  overemphasized.  Already, 
thermal  instabilities  giving  rise  to  crystallo¬ 
graphic  microcrack  growth  and  coalescence 
have  been  documented  in  Section  6  for  high 
strain  rates.  However,  irregular  microcrack 
configurations  are  a  natural  consequence  of 
low  strain  rates.  The  above  concept  of  the 
importance  of  local  temperature  increases  due 
to  straining  has  been  tested  in  recent  high 
strain  rate  experiments  on  macroscopic  speci¬ 
mens  in  two  titanium  alloys  [30],  These  have 
confirmed  that  the  heat  evolution  at  the  tip 
of  a  moving  crack  can  indeed  give  rise  to  large 
temperature  increases  in  these  alloys.  In  fact, 
a  considerable  number  of  dimples  have  been 
found  to  consist  of  solidified  regions  from 
molten  crack  ligaments,  constituting  irrefut¬ 
able  proof  of  the  occurrence  of  thermal  insta¬ 
bilities  as  discussed  above. 


9  SUMMARY  AND  CONCLUSIONS 

The  ductile  fracture  behavior  of  metals  and 
alloys  inferred  from  in  situ  experiments  in  the 
transmission  electron  microscope  has  been 
reviewed.  Direct  observations  of  dislocation 
substructures  developed  during  in  situ  frac¬ 
ture  of  foils  of  ductile  materials  have  been 
presented.  The  influence  of  several  parameters 


(including  some  associated  with  the  experi¬ 
mental  conditions  and  others  dictated  by  the 
nature  of  the  metal  or  alloy)  on  the  formation 
microstructural  features  ahead  of  the  crack 
tip  have  been  clarified.  Whereas  some  differ¬ 
ences  in  the  observations  can  be  ascribed  to 
the  nature  of  the  metal  or  alloy,  a  common 
feature  is  that  all  the  dislocation  configura¬ 
tions  discussed  exhibit  evidence  for  the  opera¬ 
tion  of  the  principle  of  energy  minimization. 
Specifically,  provided  that  the  mobility  of  the 
dislocations  is  fairly  high  and  more  than  one 
Burgers  vector  is  represented,  low  energy  cell 
structures  form . 

Microstructural  features,  such  as  crystallo¬ 
graphic  microcracks,  foil  thinning  due  to 
multiple  slip,  microcrack  nucleation  along 
cell  walls  in  thick  foils,  and  irregular  micro¬ 
cracks  due  to  the  formation  of  thinner  cell 
walls  during  very  slow  straining,  have  all  been 
illustrated  to  emphasize  the  conditions  of 
crack  formation  and  their  effect  on  the 
growth  of  cracks.  In  this,  the  influence  of 
restricted  mobility  of  dislocations  on  the 
confinement  of  slip  to  narrow  slip  zones  and 
the  resultant  thinning  to  form  voids  or  micro¬ 
cracks  has  been  clarified.  Also,  the  occasional 
gradual  thinning  of  the  foil  such  that  no  dis¬ 
locations  are  present  when  it  becomes  very 
thin,  together  with  the  resulting  brittle  be¬ 
havior,  has  been  noted. 

Two  types  of  crack  growth  behavior  were 
discussed,  namely  the  crack  growth  into  a 
region  free  from  any  dislocation  substructure 
and  the  movement  of  a  crack  into  a  region 
with  a  well-developed  stable  cell  structure. 

In  this  connection  the  importance  of  work 
hardening  in  the  ductile  fracture  behavior  of 
metals  and  alloys  has  been  emphasized.  The 
superposition  of  two  dislocation  density 
terms,  namely  the  dislocation  density  due  to 
the  applied  stress  and  the  dislocation  density 
due  to  the  crack  tip  stress  field,  has  been 
inferred.  Most  notably,  and  contrary  to  in¬ 
tuition,  it  has  been  found  that  there  should 
not  be  a  continuously  spatially  varying  cell 
size  in  the  region  of  the  crack  tip,  since  the 
theoretically  deduced  resolved  shear  stress 
component,  which  is  responsible  for  the 
movement  of  the  most  prevalent  dislocations, 
is  constant  over  a  length  of  specific  planes. 
Thus,  it  is  therefore  proposed  (and  appears  to 
be  confirmed  in  Fig.  5,  for  example)  that  a 
uniform  cell  size  will  form  along  the  planes 


of  those  orientations  in  the  region  of  constant 
shear  stress. 

In  situ  fracture  experiments  are  useful  for 
undemanding  macroscopic  ductile  fracture 
provided  that  a  proper  choice  of  experimental 
parameters  is  made.  The  direct  observation  of 
microstructural  features  in  the  transmission 
electron  microscope  to  gain  a  better  under¬ 
standing  of  ductile  fracture  is  an  invaluable 
advantage  associated  with  in  situ  fracture  ex¬ 
periments.  This  is  so,  especially  since  changes 
in  dislocation  substructures  in  the  region 
ahead  of  the  crack  tip,  in  the  plastic  zone, 
lead  to  microstructural  features  which  turn 
out  to  be  of  utmost  importance  for  the  under¬ 
standing  of  ductile  fracture  in  general. 
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